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Preface 


This monograph provides a largely self-contained and broadly accessible 
exposition of two cosmological applications of algebraic quantum field theory 
(QFT) in curved spacetime: a fundamental analysis of the cosmological evolu¬ 
tion according to the Standard Model of Cosmology and a fundamental study 
of the perturbations in Inflation. The two central sections of the book dealing 
with these applications are preceded by sections containing a pedagogical intro¬ 
duction to the subject as well as introductory material on the construction of 
linear QFTs on general curved spacetimes with and without gauge symmetry in 
the algebraic approach, physically meaningful quantum states on general curved 
spacetimes, and the backreaction of quantum fields in curved spacetimes via the 
semiclassical Einstein equation. The target reader should have a basic under¬ 
standing of General Relativity and QFT on Minkowski spacetime, but does not 
need to have a background in QFT on curved spacetimes or the algebraic ap¬ 
proach to QFT. In particular, I took a great deal of care to provide a thorough 
motivation for all concepts of algebraic QFT touched upon in this monograph, 
as they partly may seem rather abstract at first glance. Thus, it is my hope 
that this work can help non-experts to make ‘first contact’ with the algebraic 
approach to QFT. 
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Chapter 1 


Introduction 


Abstract. In this chapter we give a pedagogical introduction to algebraic quan¬ 
tum field theory and explain the concepts and the relevance of this framework 
in quantum field theory on curved spacetimes. This introduction should serve 
as a guide for the next chapter, where many concepts of and constructions in 
algebraic quantum field theory on curved spacetimes are reviewed in detail. 
Afterwards, we give a non-technical overview of the cosmological applications 
discussed in the final chapter of this monograph. 

1.1 A Pedagogical Introduction to Algebraic Quan¬ 
tum Field Theory on Curved Spacetimes 

Algebraic Quantum Field Theory (AQFT) |12j is a framework which focusses 
on the local and algebraic properties of QFT and thus aims for understanding 
structural properties of relativistic quantum field theory from first principles 
and in a model-independent fashion. In standard textbook treatments of QFT 
in Minkowski spacetime, the formalism of QFT is developed by constructing 
operators and deriving relations based on the vacuum state and the associated 
Hilbert space. However, this approach is not directly generalisable to curved 
spacetimes as we shall explain now. 

To this avail, we consider a quantized Hermitean scalar field 4>{x) and assume 
that it can be decomposed in two different ways as 



( 1 . 1 ) 


where al,ai and bl,bi are two sets of creation and annihilation operators with 
corresponding modes Ai(x), Bi{x) and vacua |Ha), |Hb) 


aj|Ho) = 0, 6i|Hf,) = 0. 
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The two sets of creating and annihilation operators are related by a Bogoliubov 
transformation 

h = atUi + f3ia\ . 

Mathematically the two possible decompositions of (j>{x) seem to be equivalent, 
whereas physically we would ask the question; Which of the two decompositions 
of 4){x) is better, or, alternatively, is there a preferred way to decompose 4'(x)l 
In Minkowski spacetime, we have Poincare symmetry at our disposal, in par¬ 
ticular Minkowski spacetime is time-translation invariant and we can construct 
a Hamilton operator H and obtain a related notion of ‘energy’. If H\Ila) = 0, 
but H\0,h) ^ 0, we would call Ifla) the ground state or vacuum state and choose 
to work with the decomposition of (j>{x) in terms of aj,ai. In other words, we 
would consider - i.e. represent - if){x) as an operator in the Fock space of the 
vacuum state. In curved spacetimes, these ideas fail in general because generic 
curved spacetimes are not time-translation invariant; prominent examples of 
such backgrounds are cosmological spacetimes with a metric line element 

ds^ =—dt^ + a{tYdxP‘, ( 1 . 2 ) 

where the function a(t) is non-constant in time. In the absence of time- 
translation invariance, no meaningful notion of a Hamilton operator or ‘energy’ 
exists, and thus we have no means to select or define a vacuum state. Even 
under these circumstances, one might still think that various possible decom¬ 
positions of the form (HU are in a sense equivalent, so that it does not really 
matter which of these one chooses to work with. However, in general one is 
facing the additional problem that 

Nha = = CC , 

i 

i.e. the ‘a-vacuum’ contains infinitely many ‘6-particles’, which in mathematical 
terms implies that |Oq) and |H{,) can not lie in the same Hilbert space. An exam¬ 
ple of this situation can be constructed by considering a cosmological spacetime 
of the form with a(t) = tanh(ct) where c is a constant with dimension of 
inverse time. The asymptotic regions of this spacetime for t —>■ ±oo are time- 
translation invariant and one can define corresponding asymptotic vacua |H±). 
One may then compute that e.g. the ‘-I-’-vacuum contains infinitely many 
particles which can be physically interpreted by saying that the expansion of 
space encoded in the functional form of a{t) creates infinitely many particles. 
This occurs because the quantum field ipi^) has infinitely many degrees of free¬ 
dom, which are all excited by the expansion. Consequently, the second equality 
sign in HU is in general purely heuristic because the two decompositions of 
4>ix) listed there are in general not related in a physically and mathematically 
meaningful way. 

From the above discussion we can infer that in the context of curved space- 
times the very notion of ‘particle’ is strictly speaking meaningless, because it 
relies on a preferred choice of vacuum state, which in general does not exist. 
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Consequently, one should not regard particles as a fundamental concept in quan¬ 
tum field theory, but at most as a derived concept which is meaningful in an 
approximative sense if the time scales on which the spacetime background is 
changing - or more general, the spacetime curvature scales - are large com¬ 
pared to the scales relevant for the physical processes we would like to discuss. 
On more conceptual grounds, we see that a Hilbert space can not be a primary 
object in the general construction of QFT models. 

Consequently, in algebraic quantum field theory one aims for constructing a 
QFT in a purely algebraic fashion without recourse to a Hilbert space. In fact, 
one first constructs an algebra which contains all observables of the theory and 
encodes their algebraic relations such as commutation relations and equations 
of motion. Hilbert spaces then appear in a second step as the spaces on which 
the algebra of observables can be represented. Different physical situations, e.g. 
states with a different temperature on Minkowski spacetime, are mathemati¬ 
cally modelled by in general inequivalent representations. Thus, the algebraic 
approach to QFT has the advantage that it enables us to discuss the physical 
properties of a physical system and the physical properties of a state of this 
system separately. Moreover, it allows us to treat all states of a physical system 
described by a QFT at once and in a coherent fashion. 

To understand the essential concepts of algebraic quantum field theory in 
curved spacetimes, it is advisable to investigate the easiest field model, the 
free Hermitean scalar field. We thus briefly sketch how a Klein-Gordon field (p 
propagating on a curved spacetime is treated in the algebraic framework. All 
concepts and notions we touch upon in the following will be explained in detail 
in the next chapter of this monograph. 

The first fundamental algebraic property of p is the Klein-Gordon equation 

Pp = (—□ + iri^ + ^R)(j) = 0 , 

where □ = is the d’Alembert operator, m is the mass, R is the Ricci 

curvature scalar and ^ quantifies a non-minimal coupling of (p to the scalar 
curvature. P is a (normally) hyperbolic differential operator, which means that 
classical solutions of the Klein-Gordon equation exist if we prescribe initial data, 
i.e. the value of p{x) and its time-derivative dt<p{x) at a fixed time t. Moreover, 
the solutions for given initial data are unique and the value of a solution p ed a 
point X depends only on the initial data in the past (or future) lightcone of x. 
Gonsequently, physical systems described by hyperbolic equations propagate in 
a causal and predictive fashion. In fact, all these statements are not correct on 
general spacetimes, but they hold on globally hyperbolic spacetimes. These are 
spacetimes which are of the form {M,g), where the spacetime manifold M may 
be decomposed as M x E, with R corresponding to ‘time’ and the manifold E 
corresponding to ‘space’, and where the causal (i.e. lightcone) structure specified 
by the metric g is such that the worldline of any physical observer, i.e. any 
inextendible timelike curve, hits an ‘equal-time surface’ {to} x Y exactly once. 
Minkowski spacetime and e.g. cosmological spacetimes are globally hyperbolic. 

Being a hyperbolic operator, the Klein-Gordon operator on a globally hyper¬ 
bolic spacetime has unique advanced and retarded Green’s functions Efi{x,y), 
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EA{x,y) = ER{y,x) which satisfy PxER{x,y) = PxEA{x,y) = S{x,y) and are 
such that Efi{x,y) {EA{x,y)) vanishes if x is not in the forward (backward) 
lightcone of y. Given these Green’s functions, we may construct the antisym¬ 
metric causal propagator E{x,y) = Efi{x,y) — EA{x,y) which is sometimes also 
called commutator function, spectral function or Pauli-Jordan distribution. The 
second fundamental algebraic property of the quantum field (j> are the canonical 
commutation relations (GGR) 

[(j){x),(l){y)] = iE{x,y). 

Clearly, E{x, y) vanishes if x and y are not causally related and thus the GGR 
encode the physical requirement that causally unrelated observables commute. 
Given coordinates {t,x) £ M = MxS one can show that ^’i, ^ 2 , ^ 2 )|ti=t 2 = 

6 {xi,X 2 ) and E(ti,xi,t 2 , ^ 2 )|ti=t 2 = 0- Consequently, the above covariant GGR 
are equivalent to equal-time GGR 

[dtfit, xi), (j){t, X 2 )] = i 6 {xi,X 2 ), [4>{t, xi), 4 >{t, £ 2 )] = 0 . 

E{x, y) is a singular object - a distribution - and diverges for x and y which 
are lightlike related. Consequently, the quantum field (j){x) is a singular object 
as well, which is rooted in the fact that it encodes infinitely many degrees of 
freedom. For this reason one often considers in the algebraic approach to QFT 
‘smeared fields’ 

Hf) = if, = j dx^J\d^\ (t){x)f{x ), 

M 

where /, called a ‘test function’, is infinitely often differentiable and has compact 
support in spacetime. Physically, / has to be interpreted as a ‘weighting func¬ 
tion’ such that (j>{f) models a ‘weighted measurement’ of the observable (f>{x). 
The compact and thus bounded support of / in spacetime reflects the physically 
realistic situation that detectors have a finite spatial size and measurements are 
performed in a finite time interval. 

The basic algebra of observables A{M) of the Hermitean scalar field </> is 
constructed by considering sums of products of smeared fields <f>{f) where / 
ranges over all possible test functions. The Klein-Gordon equation is encoded 
by identifying 4){f) with 0 if / is of the form Ph with a test function h and 
the ‘smeared GGR’ read (/>(/ 2 )] = iE{fi,f 2 ), where E{fi, f 2 ) is the 

causal propagator integrated with the test functions fi and / 2 . To have a 
notion of ‘taking the adjoint’, one introduces a *-operation specified by e.g. 
{4){f)(j){g))* = 4>{g)*4>{f)* = (j){g)(j){f). A state w on A{M) is a linear func¬ 
tional uj : A(M) C, which is positive and normalised, namely, uj{A*A) > 0 
for all A £ A{M), and a;(l) = 1. In this context, uj{A) for A £ A(M) has 
the physical interpretation of being the expectation value of A. Given an alge¬ 
braic state oj, one obtains a canonical representation tt^ of A{M) on a Hilbert 
space with vacuum via the so-called GNS construction, such that e.g. 
uj{4>{f)) = {0^\Trui{(l>{f))\0.^). In particular, the algebraic positivity condition 
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oj{A*A) > 0 ensures that the Hilbert space vector has positive norm 

for all A and the normalisation condition w(l) = 1 ensures that has unit 
norm. Conversely, given a Hilbert space H with the Klein-Gordon field re¬ 
alised as an operator (valued distribution) on H, the algebra constituted by 
these operators together with a normalised Hilbert space state are naturally 
a field algebra and a state in the abstract sense. Finally, an algebraic state 
u) on A{M) is uniquely determined, once all its n-point correlation functions 
w„(xi, • • • , Xn) = uj{(j){xi) ■ ■ ■ 4 >{xn)) are known. 

The algebra A{M) may be interpreted as the canonical quantization of a 
space of classical observables with a Poisson bracket defined by the causal prop¬ 
agator E. It is convenient to discuss physical properties such as e.g. gauge- 
invariance on the level of this classical space before passing to the quantized 
algebra A{M), because - in the simple case of free field theories - many such 
properties of the classical theory automatically carry over to the quantum the¬ 
ory. 

The algebra A{M) contains only products of quantum fields at different 
points, in particular it does not contain quantized versions of the products of 
scalar fields at the same point such as (j){x)‘^. The naive definition of as 

lima;_>j, (j){x)(j){y) leads to divergences because the two-point correlation function 
uj 2 {x,y) = uj{4>{x)4>(]j)) of any quantum state w is singular for x = y; this 
essentially follows from the singularities of the causal propagator E{x, y) and 
the fact that for every state to, uj 2 {x,y) — ui 2 {y,x) = iE{x,y) must hold on 
account of the canonical commutation relations. The singularity of uj 2 {x,y) for 
X = y is nothing but the ‘tadpole singularity’ well-known from perturbative 
QFT in Minkowski spacetime. A way to cure this singularity in Minkowski 
spacetime is to decompose the quantum field 4>{x) represented on the Hilbert 
space corresponding to a state oj into creation and annihilation operators and 
to define a normal ordered : (j){xY '-ui by ‘normal ordering’ the creation and 
annihilation operators in the naive square (j){x)‘^. One the algebraic level, this 
is equivalent to defining 

■.(t>{xf-u = iim {(t){x)4){y) - uj 2 {x,y)). (1.3) 

x^y 

If OJ is the vacuum state on Minkowski spacetime, then one finds that products 
of : (f>{x)‘^ luj at different points are well-defined and may be computed by the 
Wick theorem which implies e.g. 

w (: (j){xy : (j){yf :^) = 2 uj 2 {x, yf . 

Consequently, normal-ordered quantities form an algebra themselves. This ob¬ 
servation relies heavily on the UV-regularity properties of the Minkowski vac¬ 
uum state, i.e. loosely speaking products of : (j){x)'^ at different points are 
well-defined because the two-point correlation function of the Minkowski vac¬ 
uum is ‘singular but not too singular’. Mathematically, one has that the expres¬ 
sion oj2{x,y)‘^ is a well-defined distribution such that integrating it with any 
pair of test functions /i, /2 gives a finite result. If we consider for example the 
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massless case, then the two-point correlation function of the Minkowski vacuum 
has the form 

In order to mimic the procedure of normal ordering in Minkowski spacetime also 
in curved spacetimes, we need to consider the proper generalisation of the UV- 
regularity properties of the Minkowski vacuum. In turns out that Hadamard 
states are precisely the class of states in QFT on curved spacetimes which have 
the property that products of their correlation functions such as uj 2 {x,y)‘^ are 
well-defined. Hadamard states are characterised by having two-point functions 
of the form 




u{x,y) 

<j{x,y) 


+ V log(cr(a;, y)) + w{x, y)] = H {x, y) + 


'Ax,y) 

87r2 


where a{x, y) is one half the squared geodesic distance, and u, v, w are in¬ 
finitely often differentiable (smooth) functions. In this functional form, the UV 
divergences of uj 2 {x,y) are clearly visible, and they are completely contained 
in H{x,y). Moreover, we see that the massless Minkowski vacuum is in fact a 
Hadamard state with u = 1 and v = w = 0. It turns out that the Hadamard 
coefficients u, v, and, hence, H are completely specified by the parameters in 
the Klein-Gordon operator P and the local curvature in the neighbourhood of 
the points x and y. Hence, the singular part H is completely state-independent, 
and the two-point functions of two Hadamard states differ only in the regular 
part w. 

Following the above discussion, given any Hadamard state w, we can define 
meaningful normal ordered field expressions such as : 4){x)'^ '-ui by (II.3p . However, 
the paradigm in algebraic QFT is to define observables in a state-independent 
way and :(l){x)^ clearly fails to satisfy this property. Even worse, the regular 
part w in the correlation function uj 2 ix,y) of any Hadamard state is a highly 
non-local object because u} 2 ix,y) satisfies the Klein-Gordon equation in both 
arguments and thus is e.g. sensitive to the functional form of the metric g in the 
full past lightcone of x and y. Consequently : 4>{xY '-uj is a non-local observable, 
which is conceptually unsatisfactory. In order to cure this problem while still 
maintaining the property to have well-defined products at different points, we 
can define 

■(t>{xf-H= ^ix[i{(t>{x)(t>{y)-H{x,y)), 

x^V 

i.e. we subtract only the state-independent singular part and obtain a truly local 
observable which is independent of the geometry of spacetime far away from x. 
However, this definition of local normal ordering is not the only possibility. 
Imposing further algebraic properties such as canonical commutation relations 
with the linear field (p{x) and particular scaling and regularity properties with 
respect to the metric and the parameters in the Klein-Gordon equation, one 
finds that any expression which satisfies these conditions as well as locality 
must be of the form 


: (j){x)'^ : = : (j){x)'^ :h +{aR -f /3m^)I 
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where a and P are arbitrary dimensionless real constants. Consequently, we 
find that normal ordering - or equivalently the renormalisation of tadpoles - is 
inherently ambiguous on curved spacetimes. 

In this section we have only sketched many concepts of algebraic QFT on 
curved spacetimes and have explained them mostly at the basis of examples. 
Even though we will provide more details in the next chapter, we would like to 
stress already at this point that the algebraic construction of perturbative inter¬ 
acting models on curved spacetimes with and without local gauge symmetries is 
by now well-understood in conceptual terms [a s [n na na [m [13 [n [g n [n] . 


1.2 Outline of the Cosmological Applications 

We briefly outline the two cosmological applications of algebraic quantum held 
theory discussed in detail in the main and hnal chapter of this monograph. 

1.2.1 The Cosmological Expansion in QFT on Cnrved Space- 
times 

According to the Standard Model of Cosmology - the ACDM-model - our 
universe contains matter, radiation, and Dark Energy, whose combined energy 
density determines the expansion of the universe. In the ACDM-model, these 
three kinds of matter-energy are modelled macroscopically as a perfect huid 
and are thus completely determined by an energy density p and a pressure 
p, with different equations of state p = p{g) = wg, w = 0, 1 for matter, 

radiation and Dark Energy (assuming that the latter is just due to a cosmological 
constant) respectively. 

However, at least the contributions to the macroscopic matter and radiation 
energy densities which are in principle well-understood originate microscopically 
from particle physics. Hence, it should be possible to derive these contributions 
from hrst principles within QFT on curved spacetimes. However, in the stan¬ 
dard literature usually a mixed classical/quantum analysis is performed on the 
basis of effective Boltzmann equations in which the collision terms are computed 
within QFT on flat spacetime whereas the expansion/curvature of spacetime is 
taken into account by means of redshift/dilution-terms, see e.g. [22]. After a 
sufficient amount of cosmological expansion, i.e. in the late universe, the colli¬ 
sion terms become negligible and the energy densities of matter and radiation 
just redshift as dictated by their equation of state. 

As a first application of AQFT on curved spacetimes to cosmology, we aim 
to improve on this situation and to demonstrate that it is indeed possible to de¬ 
rive the form of the energy density in the ACDM-model microscopically within 
quantum held theory on curved spacetime: we model matter and radiation by 
quantum helds propagating on a cosmological spacetime and we show that there 
exist states for these quantum helds in which the energy density has the form 
assumed in the ACDM-model up to small corrections. Indeed, we hnd that 
these small corrections are a possible explanation for the phenomenon of Dark 
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Radiation, which shows that a fundamental analysis of the ACDM-model is 
not only interesting from the conceptual point of view but also from the phe¬ 
nomenological one. 

Due to the complexity of the problem and for the sake of clarity we shall 
make a few simplifying assumptions. On the one hand, we shall model both 
matter and radiation by scalar and neutral quantum fields for the ease of pre¬ 
sentation, but all concepts and principal constructions we shall use have been 
developed for fields of higher spin and non-trivial charge as well and we shall 
mention the relevant literature whenever appropriate. Thus, a treatment taking 
into account these more realistic fields is straightforward. On the other hand, 
we shall consider only non-interacting quantum fields and thus the effects of 
the field interactions which presumably played an important role in the early 
universe will only appear indirectly as characteristics of the states of the free 
quantum fields in our description. Notwithstanding, all concepts necessary to 
extend our treatment to interacting fields are have already been developed as 
pointed out at the end of the previous section. Finally, in this work we are only 
interested in modelling the history of the universe from the time of Big Bang 
Nucleosynthesis (BBN) until today. This restriction also justifies our approxi¬ 
mation of considering non-interacting quantum fields, as one usually assumes 
that field interactions can be neglected on cosmological scales after electron- 
positron annihilation, which happened roughly at the same time as BBN [22) . 

The quantum states which we will find to microscopically model the macro¬ 
scopic energy densities in the ACDM-model are generalised thermal excita¬ 
tions of so-called states of low energy, which are homogeneous and isotropic 
Hadamard states that minimise the energy density integrated against a weight¬ 
ing function / [25]. Whereas in e.g. Minkowski spacetime the vacuum is the 
only state of low energy in this sense, in general cosmological spacetimes states 
of low energy depend on the sampling function / and are thus non-unique as 
expected. Notwithstanding, we shall compute that for sufficiently large width of 
the sampling function /, the energy density in states of low energy on cosmolog¬ 
ical spacetimes of ACDM-type is negligible in comparison to the macroscopic 
energy density in the ACDM-model. This generalises the results found in 
for the special case of de Sitter spacetime. Consequently, states of low energy 
with sufficiently large width of their characteristic sampling function all deserve 
to be considered as ‘generalised vacuum states’, i.e. as a good approximation 
to the concept of ‘vacuum’ in cosmological spacetimes. The generalised thermal 
excitations of the states of low energy we shall consider are phenomenologically 
well-motivated. For the case of the massless, conformally coupled scalar field 
modelling radiation, they are just conformal transformations of thermal states 
in Minkowski spacetime (conformal KMS states) which phenomenologically re¬ 
flect the observation that the Cosmic Microwave Background (CMB) radiation 
is thermal. On the other hand. Dark Matter, which constitutes the major part 
of the cosmological matter density in the ACDM-model, is in many models 
considered to be of thermal origin and the quantum states we shall consider for 
the massive conformally coupled scalar field are thought to model this thermal 
origin in simple terms. 
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The analysis in Section 15^ is essentially identical to [13] and constitutes the 
first full publication of the results reported there. 

1.2.2 A Birds—Eye View of Perturbations in Inflation 

The inflationary paradigm is by now an important cornerstone of modern cos¬ 
mology. In the simplest models of Inflation, one assumes that a classical real 
Klein-Gordon field with a suitable potential V{(j)), coupled to spacetime met¬ 
ric via the Einstein equations, drives a phase of exponential expansion in the 
early universe. After this phase, the universe respectively its matter-energy 
content is thought to be almost completely homogenised, whereby the quan¬ 
tized perturbations of the scalar field and the metric are believed to constitute 
the seeds for the small-scale inhomogeneities in the universe that we observe 
today. 

Mathematically, this idea is usually implemented by considering the coupled 
Einstein-Klein-Gordon system on a Friedmann-Lemaitre-Robertson-Walker 
(FLRW) spacetime. Given a suitable potential V{(j)), this coupled system will 
have solutions which display the wanted exponential behaviour. In order to 
analyse the perturbations in Inflation, the Einstein-Klein-Gordon system is 
linearised and the resulting linear field theory is quantized on the background 
solution in the framework of quantum field theory in curved spacetimes. The 
theory of perturbations in Inflation thus constitutes one of the major applica¬ 
tions of this framework. 

However, a direct quantization of the linearised Einstein-Klein-Gordon sys¬ 
tem is potentially obstructed by the fact that this system has gauge symmetries. 
Thus the usual approach to the quantization of perturbations in Inflation, see 
e.g. the reviews IHImIEt] and the recent work [9], consists of first splitting 
the degrees of freedom of the perturbed metric into components which trans¬ 
form as scalars, vectors and tensors under the isometry group of the FLRW 
background, the Euclidean group. Subsequently, gauge-invariant linear combi¬ 
nations of these components as well as the perturbed scalar field are identified, 
which are then quantized in the standard manner. Thereby it turns out that 
the tensor components of the perturbed metric are manifestly gauge-invariant, 
whereas the vector components are essentially pure gauge and thus unphys¬ 
ical. The scalar perturbations instead are usually quantized in terms of the 
gauge-invariant Mukhanov-Sasaki variable, which is essentially a conformally 
coupled Klein-Gordon field with a time-dependent mass. In the recent work [9], 
this choice of dynamical variable has been shown to be uniquely fixed by cer¬ 
tain natural requirements. The relation between the quantized perturbations of 
the Einstein-Klein-Gordon system and the small-scale inhomogeneities in the 
present universe is usually established by relating the power-spectrum of the 
latter to the power spectrum of the former in several non-trivial steps, cf. e.g. 
HI EH EH- An approach which differs in the way this relation is made, and is 
closer to the spirit of stochastic gravity, may be found in the recent work [26j . 

The conceptual drawback of the standard approach to quantizing perturba¬ 
tions in Inflation is that this approach makes heavy use of the isometry group 
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and the related preferred coordinate system of FLRW spacetimes and is thus 
inherently non-covariant. In that sense, it is a bottom~up approach, which is 
of course well-motivated by the fact that it allows one to make explicit compu¬ 
tations. Notwithstanding, it seems advisable to check whether the same results 
can be obtained in a rather top-down approach, as this would provide a firm 
conceptual underpinning of the standard approach. Motivated by this, the 
quantum theory of the linearised Einstein-Klein-Gordon system on arbitrary 
on-shell backgrounds, and with arbitrary potential V{(f>) and non-minimal cou¬ 
pling to the scalar curvature has been developed in |14) . In order to deal with 
the gauge symmetries of this system, [14] follows ideas of [7], which deals with 
the gauge-invariant quantization of the vector potential on curved spacetimes. 
This approach was later used in m for quantizing linearised pure gravity on 
cosmological vacuum spacetimes and generalised in m in order to encompass 
arbitrary (Bosonic and Fermionic) linear gauge theories on curved spacetimes. 
In contrast to the BRST/BV approach to quantum gauge theories, see e.g. 
mm and |4] for an application to perturbative pure quantum gravity on 
curved spacetimes, the formalism used in m works without the introduction of 
auxiliary fields, at the expense of being applicable only to linear field theories. 
We shall review this general formalism to quantize linear gauge theories in Sec¬ 
tion 12.2.21 and the quantization of the linearised Einstein-Klein-Gordon system 
on arbitrary on-shell backgrounds in the first part of Section [3.31 

In the second part of that section, we consider the special case of FLRW 
backgrounds and review the results of m on comparing the quantum theory 
obtained from the general quantization of the linearised Einstein-Klein-Gordon 
system on on-shell backgrounds with the standard approach to the quantiza¬ 
tion of perturbations in Inflation. Thereby it turns out that the set of quantum 
observables in the standard approach, which is spanned by local observables of 
scalar and tensor type, is contained in the set of observables obtained in the 
general construction, but strictly smaller. However, one further finds that this 
discrepancy seems to be alleviated if one restricts to configurations of the lin¬ 
earised Einstein-Klein-Gordon system which vanish at spatial infinity, which 
apparently is a general assumption in the standard approach, see e.g. [23], 
because these configurations are considered to be ‘small’ and thus truly pertur¬ 
bative. Namely, it is argued in [14] that local observables of scalar and tensor 
type are sufficient for measuring this subset of configurations. 
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Chapter 2 

Algebraic Quantum Field 
Theory on Curved 
Spacetimes 


Abstract. In this chapter, we review all background material on algebraic 
quantum field theory on curved spacetimes which is necessary for understand¬ 
ing the cosmological applications discussed in the next chapter. Starting with a 
brief account of globally hyperbolic curved spacetimes and related geometric no¬ 
tions, we then explain how the algebras of observables generated by products of 
linear quantum fields at different points are obtained by canonical quantization 
of spaces of classical observables. This discussion will be model-independent 
and will cover both Bosonic and Fermionic models with and without local gauge 
symmetries. Afterwards, we review the concept of Hadamard states which en¬ 
compass all physically reasonable quantum states on curved spacetimes. The 
modern paradigm in QFT on curved spacetimes is that observables and their 
algebras should be constructed in a local and covariant way. We briefly review 
the theoretical formulation of this concept and explain how it is implemented 
in the construction of an extended algebra of observables of the free scalar field 
which also contains products of quantum fields at coinciding points. Finally, 
we discuss the quantum stress-energy tensor as a particular example of such an 
observable as well as the related semiclassical Einstein equation. 


2.1 Globally Hyperbolic Spacetimes and Related 
Geometric Notions 

The philosophy of algebraic quantum field theory in curved spacetimes is to set 
up a framework which is valid on all physically reasonable curved Lorentzian 
spacetimes and independent of their particular properties. Given this frame¬ 
work, one may then exploit particular properties of a given spacetime such as 
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symmetries in order obtain specific results or to perform explicit calculations. 
A class of spacetimes which encompasses most cases which are of physical in¬ 
terest are globally hyperbolic spacetimes. These include Friedmann-Lemaitre- 
Robertson-Walker spacetimes - in particular Minkowski spacetime - as well 
as Black Hole spacetimes such as Schwarzschild- and Kerr-spacetime, whereas 
prominent examples of spacetimes which are not globally hyperbolic are Anti 
de Sitter-spacetime (see e.g. [3 Chapter 3.5]) and a portion of Minkowski 
spacetime obtained by restricting one of the spatial coordinates to a finite in¬ 
terval such as the spacetimes relevant for discussing the Casimir effect. The 
constructions we shall review in the following are well-defined on all globally 
hyperbolic spacetimes. The physically relevant spacetime examples which are 
not globally hyperbolic are usually such that sufficiently small portions still have 
this property. Consequently, the algebraic constructions on globally hyperbolic 
spacetimes can be extended to these cases by patching together local construc¬ 
tions, see for instance [85l[33]. In this section we shall review the definition of 
globally hyperbolic spacetimes and a few related differential geometric notions 
which we shall use throughout this monograph. 

To this end, in this work a spacetime (M, g) is meant to be a Hausdorff, con¬ 
nected, smooth manifold M, endowed with a Lorentzian metric g, the invariant 
volume measure of which shall be denoted by dgX = y^\ det We will mostly 
consider four-dimensional spacetimes. However, most notions and results can be 
formulated and obtained for Lorentzian spacetimes with a dimension d differing 
from four and we will try to point out how the spacetime dimension affects them 
whenever it seems interesting and possible. We will follow the monograph by 
Wald [122] regarding most conventions and notations and, hence, work with the 
metric signature (—, -f, -f, -f). It is often required that a spacetime be second 
countable, or, equivalently, paracompact, i.e. that its topology has a countable 
basis. Though, as proven by Geroch in m, paracompactness already follows 
from the properties of (M, g) listed above. In addition to the attributes already 
required, we demand that the spacetime under consideration is orientable and 
time-orientable and that an orientation has been chosen in both respects. We 
will often omit the spacetime metric g and denote a spacetime by M in brief. 

For a point x G M, T^M denotes the tangent space of M at x and T*M 
denotes the respective cotangent space; the tangent and cotangent bundles of 
M shall be denoted by TM and T*M, respectively. If x : Afi — M 2 is a diffeo- 
morphism, we denote by x* the pull-back of x and by x* the push-forward of x- 
X* and X* map tensors on M 2 to tensors on Mi and tensors on Mi to tensors on 
M2, respectively; they furthermore satisfy x* = jl^^l Appendix Cj. In 

case gi and 52 are the chosen Lorentzian metrics on Mi and M 2 and x*5i = 92 ^ 
we call X an isometry; if x*9i = ^^52 with a strictly positive smooth function H, 
X shall be called a conformal isometry and Ll^g a conformal transformation of 
g. Note that this definition differs from the one often used in the case of highly 
symmetric or flat spacetimes since one does not rescale coordinates, but the met¬ 
ric. A conformal transformation according to our definition is sometimes called 
Weyl transformation in the literature. If x is an embedding x : Mi ^ M 2 , i.e. 
x(Mi) is a submanifold of M 2 and x a diffeomorphism between Mi and x(Mi), 
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it is understood that a push-forward x* of x is only defined on x(-^i) C M 2 . 

In case an embedding x : -^1 ^ M 2 between the manifolds of two spacetimes 
(Ml, 51 ) and (M 2 , 52 ) is an isometry between (Mi, 51 ) and (x(A^i), 52 |x(Mi)), 
we call X an isometric embedding, whereas an embedding which is a conformal 
isometry between (Mi,gi) and (x(Mi),g 2 lx(Mi)) shall be called a conformal 
embedding. 

Some works make extensive use of the abstract index notation, i.e. they use 
Latin indices to denote tensorial identities which hold in any basis to distinguish 
them from identities which hold only in specific bases. As this distinction will 
not be necessary in the present work, we will not use abstract index notation, 
but shall use Greek indices to denote general tensor components in a coordinate 
basis {i 9 ii}/i=o ,...,3 and shall reserve Latin indices for other uses. We employ the 
Einstein summation convention, e.g. Af = Sn=o shall lower Greek 

indices by means of = g(i9ii, d,^) and raise them by g^'^ = {g~^)^u- 

Every smooth Lorentzian manifold admits a unique metric-compatible and 
torsion-free linear connection, the Levi-Civita connection, and we shall denote 
the associated covariant derivative along a vector field v, i.e. a smooth section of 
TM, by V„. We will abbreviate by and furthermore use the shorthand 
notation for covariant derivatives of a tensor field T. 

Our definitions for the Riemann tensor Rap-yS, the Ricci tensor Rap, and the 
Ricci scalar R are 

‘^a-.P'Y '^OL’yP — Ra Py'^X 7 RaP — Ra /3A 7 ^ a 7 (^A) 

where Va are the components of an arbitrary covector. The Riemann tensor 
possesses the symmetries 

RaPyS — RpayS — RySaP : RaPyS P RaSPy “f RaySp — 0 (^■^) 

and fulfils the Bianchi identity 


(2.3) 



Moreover, its trace-free part, the Weyl tensor, is defined as 


CaPyS — RaPyS g {gaSgpy gaygps) R 

2 gpyRaS gaSRpy P gayRps') 


where the appearing coefficients differ in spacetimes with d pi A. In addition 
to the covariant derivative, we can define the notion of a Lie derivative along a 
vector field v. the integral curves c(s) of v with respect to a curve parameter 
s define, in general only for small s and on an open neighbourhood of c( 0 ), a 
one-parameter group of diffeomorphisms Xs [1221 Chapter 2.2]. Given a tensor 
field T of arbitrary rank, we can thus define the Lie derivative of T along v as 
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If Xg is a one-parameter group of isometries, we call v a Killing vector field, 
while in case of Xs being a one-parameter group of conformal isometries, we 
shall call v a conformal Killing vector field. It follows that a Killing vector field 
V fulfils Cyg — 0, while a conformal Killing vector field v fulfils C^g = fg with 
some smooth function / |1221 Appendix C.3]. 

In order to define what it means for a spacetime to be globally hyperbolic, we 
need a few additional standard notions related to Lorentzian spacetimes. To wit, 
following our sign convention, we call a vector Vx S TxM timelike if g{vx, Vx) < 0, 
spacelike if g{vx, Vx) > 0, lightlike or null if g{vx, Vx) = 0, and causal if it is either 
timelike or null. Extending this, we call a vector field v : M ^ TM spacelike, 
timelike, lightlike, or causal if it possesses this property at every point. Finally, 
we call a curve c : K. D / —;> M, with I an interval, spacelike, timelike, lightlike, 
or causal if its tangent vector field bears this property. Note that, according to 
our definition, a trivial curve c = a; is lightlike. As {M, g) is time orientable, we 
can split the lightcones in TM at all points in M into ‘future’ and ‘past’ in a 
consistent way and say that a causal curve is future directed if its tangent vector 
field at a point is always in the future lightcone at this point; past directed causal 
curves are defined analogously. 

For the definition of global hyperbolicity, we need the notion of inextendible 
causal curves; these are curves that ‘run off to infinity’ or ‘run into a singular 
point’. Hence, given a future directed curve c parametrised by s, we call x a 
future endpoint of c if, for every neighbourhood O of x, there is an sq such that 
c(s) S O for all s > sq. With this in mind, we say that a future directed causal 
curve is future inextendible if, for all possible parametrisations, it has no future 
endpoint and we define past inextendible past directed causal curves similarly. A 
related notion is the one of a complete geodesic. A geodesic c is called complete 
if, in its affine parametrisation defined by Vdc/ds^ = 0, the affine parameter s 
ranges over all R.. A manifold M is called geodesically complete if all geodesics 
on M are complete. 

In the following, we are going to define the generalisations of flat spacetime 
lightcones in curved spacetimes. By I~^(x,M) we denote the chronological fu¬ 
ture of a point x relative to M, i.e. all points in M which can be reached by 
a future directed timelike curve starting from x, while J^(x,M) denotes the 
causal future of a point x, viz. all points in M which can be reached by future 
directed causal curve starting from x. Notice that, generally, x S J^{x,M) 
and I~^{x,M) is an open subset of M while the situations x ^ I^{x,M) and 
J'^{x,M) being a closed subset of M are not generic, but for instance present 
in globally hyperbolic spacetimes [122] . In analogy to the preceding definitions, 
we define the chronological past I~{x, M) and causal past J~ {x, M) of a point x 
by employing past directed timelike and causal curves, respectively. We extend 
this definition to a general subset O C M by setting 

I^{0,M)= y I^{x,M) J^{0,M)= y J^{x,M); 

xeo xeo 

additionally, we define I{0, M) = I+{0, M)U/-(C>, M) and J(C>, M) = J+{0, M)U 
J~ {O, M). As the penultimate prerequisite for the definition of global hyper- 
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bolicity, we say that a subset O of M is achronal if I^{0, M) fl O is empty, i.e. 
an achronal set is such that every timelike curve meets it at most once. Given 
a closed achronal set O, we define its future domain of dependence D^{0,M) 
as the set containing all points x G M such that every past inextendible causal 
curve through x intersects O. By our definitions, D~^{0,M) C J^{0,M), but 
note that M) is in general considerably larger than D^{0, M). We define 

D-{0,M) analogously and set D{0,M) = D+{O, M) U D-(O, M). D{0,M) 
is sometimes also called the Cauchy development of O. With this, we are fi¬ 
nally in the position to state the definition of global hyperbolicity (valid for all 
spacetime dimensions). 

Definition 2.1.1. A Cauchy surface is a closed achronal set E C M with 
I?(E, M) = M. A spacetime {M,g) is called globally hyperbolic if it contains a 
Cauchy surface. 

Although the geometric intuition sourced by our knowledge of Minkowski 
spacetime can fail us in general Lorentzian spacetimes, it is essentially satisfac¬ 
tory in globally hyperbolic spacetimes. According to Definition 12.1.11 a Cauchy 
surface is a ‘non-timelike’ set on which every ‘physical signal’ or ‘worldline’ 
must register exactly once. This is reminiscent of a constant time surface in 
flat spacetime and one can indeed show that this is correct. In fact, Geroch has 
proved in m that globally hyperbolic spacetimes are topologically K x E and 
Bernal and Sanchez [nmniin] have been able to improve on this and to show 
that every globally hyperbolic spacetime has a smooth Cauchy surface E and is, 
hence, even diffeomorphic to ffi. x E. This implies in particular the existence of 
a (non-unique) smooth global time function < : M —>■ R, i.e. t is a smooth func¬ 
tion with a timelike and future directed gradient field Vt; t is, hence, strictly 
increasing along any future directed timelike curve. In the following, we shall 
always consider smooth Cauchy surfaces, even in the cases where we do not 
mention it explicitly. 

In the remainder of this chapter, we will gradually see that globally hyper¬ 
bolic curved spacetimes have many more nice properties well-known from flat 
spacetime and, hence, seem to constitute the perfect compromise between a 
spacetime which is generically curved and one which is physically sensible. Par¬ 
ticularly, it will turn out that second order, linear, hyperbolic partial differential 
equations have well-defined global solutions on a globally hyperbolic spacetime. 
Hence, whenever we speak of a spacetime in the following and do not explicitly 
demand it to be globally hyperbolic, this property shall be understood to be 
present implicitly. 

On globally hyperbolic spacetimes, there can be no closed timelike curves, 
otherwise we would have a contradiction to the existence of a smooth and strictly 
increasing time function. There is a causality condition related to this which 
can be shown to be weaker than global hyperbolicity, namely, strong causality. 
A spacetime is called strongly causal if it can not contain almost closed timelike 
curves, i.e. for every x G M and every neighbourhood Oi ^ x, there is a 
neighbourhood O 2 C Ox oi x such that no causal curve intersects O 2 more than 
once. One might wonder if this weaker condition can be filled up to obtain full 
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global hyperbolicity and indeed some references, e.g. [51 [5S]. define a spacetime 
(M, g) to be globally hyperbolic if it is strongly causal and J~^ {x) fl J~ (y) is 
compact for all x,y £ M. One can show that the latter definition is equivalent 
to Definition l2.1.1l [51 1122] which is, notwithstanding, the more intuitive one in 
our opinion. 

We close this section by introducing a few additional sets with special causal 
properties. To this avail, we denote by exp^, the exponential map at x € M. A 
set O C M is called geodesically starshaped with respect to a; S O if there is 
an open subset O' of T^M which is starshaped with respect to 0 G T^M such 
that exp,j, : O' ^ O is a diffeomorphism. We call a subset O C M geodesically 
convex if it is geodesically starshaped with respect to all its points. This entails 
in particular that each two points x, y in O are connected by a unique geodesic 
which is completely contained in O. A related notion are causal domains, these 
are subsets of geodesically convex sets which are in addition globally hyperbolic. 
Finally, we would like to introduce causally convex regions, a generalisation 
of geodesically convex sets. They are open, non-empty subsets O C M with 
the property that, for all x,y € O, all causal curves connecting x and y are 
entirely contained in O. One can prove that every point in a spacetime lies 
in a geodesically convex neighbourhood and in a causal domain m and one 
might wonder if the case of a globally hyperbolic spacetime which is geodesically 
convex is not quite generic. However, whereas Friedmann-Lemaitre-Robertson- 
Walker spacetimes with flat spatial sections are geodesically convex, even de 
Sitter spacetime, which is both globally hyperbolic and maximally symmetric 
and could, hence, be expected to share many properties of Minkowski spacetime, 
is not. 


2.2 Linear Classical Fields on Curved Spacetimes 

As outlined in Section [O] the ‘canonical’ route to quantize linear classical field 
theories on curved spacetimes in the algebraic language is to first construct the 
canonical covariant classical Poisson bracket (or a symmetric equivalent in the 
case of Fermionic theories) and then to quantize the model by enforcing canon¬ 
ical (anti)commutation relations defined by this bracket. In this section, we 
shall first review how this is done for free field theories without gauge symmetry 
before discussing the case where local gauge symmetries are present. 


2.2.1 Models without Gauge Symmetry 

We shall start our discussion of classical field theories without local gauge sym¬ 
metries by looking at the example of the free Klein-Gordon field, which is the 
‘harmonic oscillator’ of QFT on curved spacetimes. In discussing this example it 
will become clear what the basic ingredients determining a linear field theoretic 
model are and how they enter the definition and construction of this model in 
the algebraic framework. 
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The Free Neutral Klein Gordon Field 

In Physics, we are used to describe dynamics by (partial) differential equations 
and initial conditions. The relevant equation for the neutral scalar field is the 
free Klein-Gordon equation 


(-□ +/)<^ = + /)(/) = 0 


with the d’Alembert operator □ and some scalar function / of mass dimension 
2. The function / determines the ‘potential’ V{(j)) = ^f4>^ of the Klein-Gordon 
field and may be considered as a background field just like the metric g. Usually 
one considers the case / = i.e. 

P(l3 = {-n + ^R + m"^) (j) = 0 (2.4) 


such that / is entirely determined in terms of a constant mass m > 0 and the 
Ricci scalar R where the dimensionless constant ^ parametrises the strength of 
the coupling of (j) to R. In principle one could consider more non-trivial coupling 
terms with the correct mass dimension such as / = {R°‘^RajVf'/{ vtAR), however 
these may be ruled out by either invoking Occam’s razor or by demanding 
analytic dependence of / on the metric and m like in [701 Section 5.1]. 

The case (12.41) with ^ = 0 is usually called minimal coupling, whereas, in 
four dimensions, the case ^ = 1/6 is called conformal coupling. While the 
former name refers to the fact that the Klein-Gordon field is coupled to the 
background metric only via the covariant derivative, the reason for the latter 
is rooted in the behaviour of this derivative under conformal transformations. 
Namely, if we consider the conformally related metrics g and g = Qfg with 
a strictly positive smooth function 12, denote by V, □, and R the quantities 
associated to g and by V, □, and R the quantities associated to g, then the 
respective metric compatibility of the covariant derivatives V and V and their 
agreement on scalar functions imply [1221 Appendix D] 


1 


1 




(2.5) 


This entails that a function cj) solving (—□ + gi?)</ = 0 can be mapped to a 
solution (j) of (—□ + ^R)4> = 0 by multiplying it with the conformal factor U 
to the power of the conformal weight —1, i.e. (f> = We shall therefore 

call a scalar field f) with an equation of motion (—□ + ^R)(l> = 0 conformally 
invariant. In other spacetimes dimensions d ^ A, the conformal weight and the 
magnitude of the conformal coupling are different, see [1221 Appendix D]. 

Having a partial differential equation for a free scalar field at hand, one would 
expect that giving sufficient initial data would determine a unique solution on all 
M. However, this is, in case of the Klein-Gordon operator at hand, in general 
only true for globally hyperbolic spacetimes. To see a simple counterexample, 
let us consider Minkowski spacetime with a compactified time direction and 
the massless case, i.e. the equation (9/ — — dy — df)(j) = 0. Giving initial 
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conditions (l)\t=o = 0, dt<p\t=o = 1, a possible local solution is (p = t. But this 
can of course never be a global solution, since one would run into contradictions 
after a full revolution around the compactified time direction. 

In what follows, the fundamental solutions or Green’s functions of the Klein- 
Gordon equation shall play a distinguished role. Before stating their existence, 
as well as the existence of general solutions, let us define the function spaces we 
shall be working with in the following, as well as their topological duals, see e.g. 
[22l Chapter VI] for an introduction. 

Definition 2.2.1. By r(M) = C^{M,'R) we denote the smooth (infinitely 
often continuously differentiable), real-valued functions on M equipped with 
the usual locally convex topology, i.e. a sequence of functions /„ S T{M) is 
said to converge to f £ T{M) if all derivatives of /„ converge to the ones of / 
uniformly on all compact subsets of M. 

The space ro(M) = C'q“(M,R) is the subset of T{M) constituted by the 
smooth, real-valued functions with compact support. We equip ro(M) with 
the locally convex topology determined by saying that a sequence of functions 
fn £ To{M) converges to f £ ro(M) if there is a compact subset K C M such 
that all fn and / are supported in K and all derivatives of fn converge to the 
ones of / uniformly in K. 

By Tq (M) = Tq{M) C, T^{M) = T{M) (g)R C we denote the complexifi- 
cations of ro(M) and T{M) respectively. 

The spaces rsc(M) and rtc(M) denote the subspaces of T{M) consisting 
of functions with spacelike-compact and timelike-compact support respectively. 
I.e. supp / n E is compact for all Cauchy surfaces E of (M, g) and all / G 
rsc(M), whereas for all / £ rtc(M) there exist two Cauchy surfaces Ei, E 2 
with supp / C J“(Ei, M) n J+(E 2 , M). 

By rQ(M) we denote the space of distributions, i.e. the topological dual of 
Tq{M) provided by continuous, linear functionals ro(M) —>■ R, whereas r'(M) 
denotes the topological dual of r(M), i.e. the space of distributions with compact 
support. and rQ*'(M) denote the complexified versions of the real-valued 


spaces. 


For / G r(M) and u G Tq{M) D r(M) D To{M) with compact overlapping 
support, we shall denote the (symmetric and non-degenerate) dual pairing of / 
and u by 



M 


The physical relevance of the above spaces is that functions in ro(M), so- 
called test functions, should henceforth essentially be viewed as encoding the 
localisation of some observable in space and time, reflecting the fact that a de¬ 
tector is of finite spatial extent and a measurement is made in a finite time 
interval. From the point of view of dynamics, initial data for a partial differen¬ 
tial equation may be encoded by distributions or functions with both compact 
and non-compact support, whereas solutions of hyperbolic partial differential 
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equations like the Klein-Gordon one are typically distributions or smooth func¬ 
tions which do not have compact support on account of the causal propagation 
of initial data; having a solution with compact support in time would entail 
that data ‘is lost somewhere’. Moreover, fundamental solutions of differential 
equations will always be singular distributions, as can be expected from the fact 
that they are solutions with a singular i5-distribution as source. Finally, since 
(anti)commutation relations of quantum fields are usually formulated in terms 
of fundamental solutions, the quantum fields and their expectation values will 
also turn out to be singular distributions quite generically. Physically this stems 
from the fact that a quantum field has infinitely many degrees of freedom. 

Let us now state the theorem which guarantees us existence and properties 
of solutions and fundamental solutions (also termed Green’s functions or prop¬ 
agators) of the Klein-Gordon operator P. We refer to the monograph for 
the proofs. 

Theorem 2.2.1. Let P : r(M) —>■ T{M) be a normally hyperbolic operator on 
a globally hyperbolic spacetime {M,g), i.e. in each coordinate patch of M, P 
can be expressed as 

P = + A^^d^ + B 

with smooth functions A^, B and the metric principal symbol —g^’^dfj_d^. Then, 
the following results hold. 

1. Let f S Vo{M), let E be a smooth Cauchy surface of M, let {uo,ui) G 
ro(S) X ro(S), and let N be the future directed timelike unit normal vector 
field o/S. Then, the Cauchy problem 

Pu = f, u|e = Uo, VAru|E = Mi 

has a unique solution u G T{M). Moreover, 

supp u C J (supp / U supp uq U supp ui,M ). 

A unique solution to the Cauchy problem also exists if the assumptions on 
the compact support of f, uq and ui are dropped. 

2. There exist unique retarded Er and advanced Ea fundamental solutions 
(Green’s functions, propagators) of P. Namely, there are unique contin¬ 
uous maps Eji/a ■ Tq{M) —>■ T{M) satisfying P o Eji/a = E^ijA o P = 
idro(M) and supp En/Af C J^{supp f,AI) for all f G ro(M). 

3. Let f, q G ro(M). If P js formally selfadjoint, i.e. {f,Pg) = {Pf,q), then 
Efi and Ea are the formal adjoints of one another, namely, (/, E^^jAg') = 

{pA/nf, g)- 

f. The causal propagator (Pauli-Jordan function) of P defined as E = En — 
Ea is a continuous map Tq{VI) -a rsc(M) C T{M) satisfying: for all 
solutions u of Pu = 0 with compactly supported initial conditions on a 
Cauchy surface there is an f G To{M) such that u = Ef. Moreover, 
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for every f € ro(M) satisfying Ef = 0 there is a g € To{M) such that 
f = Pg. Finally if P is formally self-adjoint, then E is formally skew- 
adjoint, i.e. if, Eg) = -{Ef,g) . 

The Klein-Gordon operator P is manifestly normally hyperbolic. Moreover, 
one can check by partial integration that P is also formally self-adjoint. Hence, 
all above-mentioned results hold for P. 

By continuity and the fact that r(M) C rQ(M), the operators En/j^ and E 
define bi-distributions Eh/a, E G rQ(M^) which we denote by the same symbol 
via e.g. 

ER/A{f, 9 ) = {f,ER/Ag) = / dgx dgy ER/A{x,y)f{x)g{y). 

M2 

In terms of integral kernels of these distributions, some of the identities stated 
in Theorem 12.2.11 read 


PxER/A{x,y) = 6{x,y), EA{x,y) = ER{y,x), E{y,x) =-E{x,y). 

The support properties of E^/j^ entail that E{f,g) vanishes if the supports 
of / and g are spacelike separated. On the level of distribution kernels, this 
implies that E(x, y) vanishes for spacelike separated x and y. In anticipation 
of the quantization of the free Klein-Gordon field, this qualifies E(x, y) as a 
commutator function. In the classical theory instead, E(x,y) defines a Poisson 
bracket or symplectic form. To see this, we first need to specify the vector space 
on which this bracket should be evaluated. 

Definition 2.2.2. By Sol (Sobc) we denote the space of real (spacelike-compact) 
solutions of the Klein-Gordon equation 

Sol = {(I)G r(M) I= 0}, Soisc = Sol n rsc(M). 

By £ we denote the quotient space 

£: = ro(M)/p[ro(M)], 

which is the labelling space of linear on-shell observables of the free neutral 
Klein-Gordon field. 

The fact that £ is the labelling space of (classical) linear on-shell observables 
of the free neutral Klein-Gordon field follows from the observation that each 
equivalence class [/] G £ defines a linear functional on Sol by 

Sol 9 </)^%](<^) = (/,</)), 

where we note that, in the classical theory, Sol plays the role of the space of 
pure states of the model. As ^ is a solution of the Klein-Gordon equation 
(</)) does not depend on the representative / G [f] and is well-defined. The 
observable (/, (f) may be interpreted as the ‘smeared classical field’ <j){f) ~ (/, (/>). 
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The classical observable i.e. the observable that gives the value of a 

configuration </) at the point x, may be obtained by formally considering (l>{f) 
with f = 5x- 

We know that every (j) G Sol is in one-to-one correspondence with initial 
data given on an arbitrary but fixed Cauchy surface S of {M,g). Analogously 
the support of a representative f G [f] G £ can be chosen to lie in an arbitrarily 
small neighbourhood of an arbitrary Cauchy surface. 

Lemma 2.2.1. Let [f] G £ be arbitrary and let E be any Cauchy surface of 
(M,g). Then, for any bounded neighbourhood 0(E) of E, we can find a g G 
Tq{M) with supp g C 0(E) and g G [/]. 

Proof. Let us assume that 0(E) lies in the future of supp /, i.e. J~ (supp /, M)n 
0(E) = 0, the other cases can be treated analogously. Let us consider two 
auxiliary Cauchy surfaces Ei and E 2 which are both contained in 0(E) and 
which are chosen such that E 2 lies in the future of E whereas Ei lies in the past 
of E. Moreover, let us take a smooth function y G r(M) which is identically 
vanishing in the future of E 2 and fulfils y = 1 in the past of Ei and let us define 
g = f — PxEnf ■ By construction and on account of the properties of both a 
globally hyperbolic spacetime (M, g) and a retarded fundamental solution En 
on M, xEnf has compact support, hence g G [/]. Finally, supp g is contained 
in a compact subset of J+(supp /, M) fl 0(E). □ 

We now observe that the causal propagator E induces a meaningful Poisson 
bracket on £. 

Proposition 2.2.1. The tuple {£,t) with r : £ x £ —>■ R defined by 

t([/]: [ 5 ]) = if, Eg) 

is a symplectic space. In particular 

1. T is well-defined and independent of the chosen representatives, 

2. T is antisymmetric, 

3. T is (weakly) non-degenerate, i.e. 'r([/],[g]) = 0 for all [g] G £ implies 

[/] = [ 0 ]. 

Proof. T is independent of the chosen representatives because P o E = 0. r is 
antisymmetric because E is formally skew-adjoint, cf. the last item of Theorem 
12.2.11 and because (•,•) is symmetric. The non-degeneracy of r follows again 
from the last item of Theorem l2.2.1l and the fact that (•, •) is non-degenerate. □ 

In standard treatments on scalar field theory, one usually defines Poisson 
brackets at ‘equal times’, but as realised by Peierls in m. one can give a 
covariant version of the Poisson bracket which does not depend on a splitting 
of spacetime into space and time, and this is what we have given above. To 
relate the covariant form r to an equal-time version, we need the definition of 
a ‘future part’ of a function / G r(M). 
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Definition 2.2.3. We consider a temporal cutoff function x of the form dis¬ 
cussed in the proof of Lemma r2.2.11 i.e. a smooth function x which is identically 
vanishing in the future of some Cauchy surface E 2 and identically one in the 
past of some Cauchy surface Ei in the past of E 2 . Given such a y, we define 
for an arbitrary / G T{M) the future part f'^ and the past part f~ by 

/"^ = (i-x)/, f~ = xf- 

The relation of the covariant picture to the equal time-picture can be now 
shown in several steps. 

Theorem 2.2.2. Let (•,-)goi be defined on tuples of solutions with compact 
overlapping support by 

Sol X Sol 3 {< 1 ) 1 , ( 1 ) 2 ) (<^1, (?i'2)sol = > 2 ) • 

Moreover, let E be an arbitrary Cauchy surface of {M, g) with future-pointing 
unit normal vectorfield N and canonical measure dE induced by dgX . 

1. The causal propagator E : To{M) —>■ TsdM) descends to a bijective map 
E:£^ Sohc. 

2. (•,‘)goj is antisymmetric and well-defined on all tuples of solutions with 
compact overlapping support, in particular this bilinear form does not de¬ 
pend on the choice of cutoff y entering the definition of the future part. 

3. For all f G To{M) and all (f G Sol, (/,</)) = {Ef,(j))g^^. In particular, 
(•,-)soi is well-defined on all tuples of solutions with spacelike-compact 
overlapping support. 

4 . For all f,g G Tq{M), t([/], [ 5 ]) = {Ef,Eg)^^y, thus the causal propagator 
E : Tq(M) —>■ VsciM) descends to an isomorphism between the symplectic 
spaces {£,t) and (Sobc, (•, •)goi). 

5. For all (j)i,(j )2 S Sol with spacelike-compact overlapping support, 

</>2)gol ^ y N>^j^{(i)i,(j)2) , jfj.{(flA2) = (t)lV^,(j)2 - . 

E 

6. For all f G ro(E) it holds 

VNEf\^ = f, Ef\^=0. 

On the level of distribution kernels, this entails that 

VNE{x,y)\^xT. = S^{x,y), E{x,y)\^x^ = 0, 

where is the S-distribution with respect to the canonical measure on E. 
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Proof. We sketch the proof. The first statement follows from the last item of 
Theorem l2.2.1l The fact that (•, •)gQj is well-defined follows from the observation 
that two different definitions 4>'^, fC' of the future part differ by a compactly 
supported smooth function f = (j)f consequently the supposedly different 

definitions of the bilinear form differ by (/> 2 )goi — ^ 2)301 = {Pf, 4 ’ 2 ) = 

{f,P(j) 2 ) = 0. Note that this partial integration is only possible because / 
has compact support, in particular, (•,')gQj is non-vanishing in general. The 
antisymmetry of (•, •)gQj follows by similar arguments and P(f)'^ = P{(f) — (j)~) = 
—P(f)~. The third statement follows from the fact that Enf is a valid future 
part of Ef, thus {Ef, 0)goi = {PEr/, (jf) = (/, cf). The fourth statement follows 
immediately from the first and third one, whereas the fifth one follows from 
V^j^((/>i, (/> 2 ) = (t> 2 P(t>i — 4 ’iP4’2 by an application of Stokes theorem, see e.g. 
[38] , where also a proof of the last statement can be found. □ 

We now interpret the previous results. As argued above, elements [/] S £ 
label linear on-shell observables </<(/) — (/, </>), i.e. the classical field (p smeared 
with the test function /. The causal propagator E induces a non-degenerate 
symplectic form r on £, which we may interpret as {(/>(/), ^( 5 )} ~ [g]) = 

{f,Eg), or, formally, as {(j){x), (j){y)} = E{x,y). On the other hand, since (£,t) 
and (Solsc, (‘j ■)soi) symplectically isomorphic, we can equivalently label lin¬ 
ear on-shell observables by Sobc 3 u, i.e. by (u, the classical field ‘sym¬ 
plectically smeared’ with the test solution m, where this symplectic smearing con¬ 
sists of integrating a particular expression at equal times. The last result of the 
above theorem implies that the covariant Poisson bracket {0(x), (piy)} = E(x, y) 
has the well-known equal-time equivalent 

{Viv0(x)|E,(/)(?/)|E} = VjvA(x,y)|ExE = s^{x,y) , 
{(p{x)\^,(p{y)\^} = E{x, 2 /)|exe = 0 , 

which may be interpreted as equal-time Poisson brackets of the field 4>{x) and 
its ‘canonical momentum’ Vn(P{x). Further details on the relation between the 
equal-time and covariant picture can be found e.g. in [1231 Chapter 3]. 

General Models without Gauge Symmetry 

The previous discussion of the classical free neutral Klein-Gordon field revealed 
the essential ingredients defining this model. Following e.g. nniiie], this can 
be generalised to define an arbitrary linear field-theoretic model on a curved 
spacetime. 

Definition 2.2.4. A real Bosonic linear field-theoretic model without local 
gauge symmetries on a curved spacetime is defined by the data (M, V, P), where 

1. M ~ {M,g) is a globally hyperbolic spacetime, 

2 . V is a real vector bundle over M, the space of smooth sections r(V) of 
V is endowed with a symmetric and non-degenerate bilinear form (•, •)y 
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which is well-defined on sections with compact overlapping support and 
given by the integral of a fibrewise symmetric and non-degenerate bilinear 
form ((•,.))v :r(V) xr(V)^r(M), 

3. P : r(V) —>■ r(V) is a Green-hyperbolic partial differential operator, i.e. 
there exist unique advanced E^ and retarded E^ fundamental solutions 
of P which satisfy P o E^j^ = E^^j^ o P = id|r(,(v) and supp E^^^f C 
J=^(supp /, M) for all f €To{V); moreover P is formally self-adjoint with 
respect to (•, •)y. 

A real Fermionic linear field-theoretic model without local gauge symmetries 
on a curved spacetime is defined analogously with the only difference being that 
((•, ■))^ and (•, ■)^ are not symmetric but antisymmetric. Complex theories can 
be obtained from the real ones by complexification. 

The relevance of the given data is as follows. Classical configurations $ of 
the linear field model under consideration are smooth sections $ G r(V) of the 
vector bundle V. We recall that V is locally of the form AI x V with a real 
vector space V which implies that locally $ is a smooth function from M to V, 
see e.g. [SH [51] for background material on vector bundles. We shall denote by 
ro(V), rtc(V), rsc(V) the subspaces of ro(V) consisting of smooth sections of V 
with compact, timelike-compact and spacelike-compact support, respectively. 

The operator P specifies the equation of motion for the field model, the 
formal self-adjointness of P is motivated by the fact that equations of motion 
arising as Euler-Lagrange equations of a Lagrangean are generally given by a 
formally self-adjoint P. In fact the (formal) action S'($) = ^ (<I>,P<i>)y leads to 
the Euler-Lagrange equation Pd) = 0. 

In the Klein-Gordon case we are dealing with an operator which is normally 
hyperbolic, i.e. the leading order term is of the form As reviewed 

in Theorem 12.2.11 this operator has a well-defined Cauchy problem, i.e. it is 
Cauchy-hyperbolic^ and consequently unique advance and retarded fundamental 
solutions exist such that the operator is Green-hyperbolic. Example of partial 
differential operators which are Cauchy-hyperbolic, but not normally hyperbolic 
are the Dirac operator and the Proca operator which defines the equation of 
motion for a massive vector field, see e.g. [6]. On the other hand, the distinction 
between Cauchy-hyperbolic operators and Green-hyperbolic operators does not 
matter in most examples although one can construction operators which are 
Green-hyperbolic but not Cauchy-hyperbolic, cf. [ 6 ] for details. 

Based on the data given in Definition 12.2.41 a symplectic space (Bosonic 
case) or inner product space (Fermionic case) can be constructed in full analogy 
to the Klein-Gordon case, in particular, the following can be shown. 

Theorem 2.2.3. Under the assumptions of Definition \2. 2^ let E^ = E^ — E^ 
denote the causal propagator of P, and let Sol C r(V), Sobc C Fsc(V) denote 
the space of smooth (smooth and spacelike-compact) solutions of P^ = 0 . 

1. The tuple where 


£: = Fo(v)/p[ro(v)], 
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t:£x£^R, T{[f],[g]) = {f,EPg)^ , 

is a well-defined symplectic (Bosonic case) or inner product (Fermionic 
case) space. In particular t is well-defined and independent of the chosen 
representatives and moreover non-degenerate and antisymmetric (Bosonic 
case) or symmetric (Fermionic case). 

2. Let [/] £ £ be arbitrary and let S be any Cauchy surface of {M, g). Then, 
for any bounded neighbourhood 0{E) ofE, we can find a g G ro(V) with 
supp g C 0{E) and g G [/]. 

3. The causal propagator E^ : ro(V) —?> rsc(V) descends to a bijective map 
£ —>■ Solsc and for all f G ro(V) and all $ G Sol, 

where for all $i, $2 G Sol with spacelike-compact overlapping support, the 
bilinear form (■,')goi *5 defined as 

($l,$2)sol= • 

4 . (■,')goi may be computed as a suitable integral over an arbitrary but fixed 

Cauchy surface S of {M,g) with future-pointing normal vector field N and 
induced measure dE. If there exists a ‘current’ j : r(V) x r(V) T*M 
such that = {{^i,P^ 2 ))v - {{^ 2 ,P^i))v for all $i,<l >2 G 

r(V), then 

E 

5. The tuple (Solsc, (•, ■)goi) is a well-defined symplectic (Bosonic case) or 
inner product (Fermionic case) space which is isomorphic to {£,t). 

As with the Klein-Gordon field, the last statement implies in physical terms 
that the symplectic respectively inner product space can be constructed both 
in a covariant and in an equal-time fashion, and that the two constructions 
give equivalent results. In many cases, the equal-time point of view is better 
suited for practical computations and for proving particular further properties 
of the bilinear form r, cf. the following discussion of theories with local gauge 
invariance. 

2.2.2 Models with Gauge Symmetry 

The discussion of linear field theoretic models with local gauge symmetries on 
curved spacetimes is naturally more involved than the case where such symme¬ 
tries are absent. However, as in this monograph we will only be dealing with 
linear models and simple observables, it will not be necessary to introduce aux¬ 
iliary fields like in the BRST/BV formalism [^|^. Instead, we shall review 
an approach which has been developed in [39] for the Maxwell field, used for 
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linearised gravity in [42] and then further generalised to arbitrary linear gauge 
theories in . For linear models and simple observables this approach and the 
BRST/BV formalism give equivalent results, however, non-linear models and 
more general observables are not tractable in the way we shall review in the 
following. 


A Toy Model 


We outline the essential ideas of this approach at the example of a toy model. 
We consider as a gauge field $ = C r(V) a tuple of two scalar fields 

on a spacetime (M, g) satisfying the equation of motion 


= 





= 0 , 


where V is the (trivial) vector bundle V = M x The gauge transformations 
are given by the following translations on configuration space $ i—!> $+Ae, where 
the gauge transformation operator K : T{M) —>■ r(V) is the linear operator 
defined by Ke = (e, e)‘ for a smooth function e G T{M). One may check that 
P o K = 0 holds which is equivalent to the gauge-invariance of the action 
5'($) = with = JMdgX +4>2(t>2)- 

Clearly, the linear combination ■0 = — (/)2 is gauge-invariant and satisfies 

—00 = 0 , and it would be rather natural to quantize $ by directly quantiz¬ 
ing 0 as a massless, minimally coupled scalar field. This would be much in 
the spirit of the usual quantization of perturbations in Inflation, where gauge- 
invariant linear combinations of the gauge field components, e.g. the Bardeen- 
Potentials or the Mukhanov-Sasaki variable, are taken as the fundamental fields 
for quantization, see the last chapter of this monograph. However, in general it 
is rather difficult to directly identify a gauge-invariant fundamental field like 0 
whose classical and quantum theory is equivalent to the classical and quantum 
theory of the original gauge field. Notwithstanding, an indirect characterisa¬ 
tion of such a gauge-invariant linear combination of gauge-field components, 
which can serve as a fundamental field for quantization, is still possible. In 
the toy model under consideration we consider a tuple / = (/i,/ 2 ) G ro(V) 
of test functions fi G Tq{M). We ask that f = fi + f 2 — 0, where 
At : r(V) ^ r(M) is the adjoint of the gauge transformation operator K 
i.e. Jj^^dgX eK^f = {Ke,f)v. Clearly, any / satisfying this condition is of 
the form / = {h, —h)* for a test function h. We now observe that the pair¬ 
ing between a gauge field configuration $ and such an / is gauge-invariant, 
i.e. (4) -I- Ae,/)v = (/,‘h)v + Jj^dg^ Thus we can consider 

the ‘smeared field’ $(/) ~ (/, $)v, with / = (/i, —/i)‘ and arbitrary h, as a 
gauge-invariant linear combination of gauge-field components which is suitable 
for playing the role of a fundamental field for quantization. We can compute 
(/,^)v = iMdg X 0/i, and observe that, up to the ‘smearing’ with /i, this in¬ 
direct choice of gauge-invariant fundamental field is exactly the one discussed 
in the beginning. If one chooses h to be the delta distribution S{x,y) rather 
than a test function, one even finds (/, $) = 0 (a;), whereas for general h, (/, $) 
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can be interpreted as a weighted, gauge-invariant measurement of the held con- 
hguration $. Moreover, as already anticipated, in general gauge theories with 
more complicated gauge transformation operators K it it usually extremely dif- 
Hcult to classify all solutions of f = 0 , which would be equivalent to a direct 
characterisation of one or several fundamental gauge-invariant Helds such as ip, 
whereas working implicitly with the condition f = 0 is always possible. 

General Models 

From the previous discussion we can already infer most of the additional data 
which is needed in addition to the data mentioned in Dehnition 12.2.41 in order 
to specify a linear Held theoretic model with local gauge symmetries on curved 
spacetimes. 

Definition 2.2.5. A real Bosonic linear Held-theoretic model with local gauge 
symmetries on a curved spacetime is deHned by the data (M, V, W, P, K), where 

1. M ~ {M,g) is a globally hyperbolic spacetime, 

2 . V and W are real vector bundles over M, the spaces of smooth sec¬ 
tions V and W are endowed with symmetric and non-degenerate bilin¬ 
ear forms (•, •)y and (•, •)yy which are well-deHned on sections with com¬ 
pact overlapping support and given by the integral of Hbrewise symmetric 
and non-degenerate bilinear forms ((•,-))v : r(V) x r(V) —>■ r(M) and 
((•,-))T^:r(w)xr(w)^r(M), 

3. P : r(V) — >■ r(V) is a partial differential operator which is formally self- 
adjoint with respect to (•, •)y, 

4. K : r(W) —?> r(V) is a partial differential operator such that P o K = 0; 
moreover R = iff o K : r(W) —!> r(W) is Cauchy-hyperbolic and there 
exists an operator T : r(>V) —?> r(V) such that A) P = P + T o AT'I' is 
Green-hyperbolic and b) Q = o T is Cauchy-hyperbolic. 

A real Fermionic linear Held-theoretic model with local gauge symmetries on 
a curved spacetime is deHned analogously with the only difference being that 
((■) ■))v’ (■’ ■)v> ■))w (■’ ■)w symmetric but antisymmetric. Com¬ 

plex theories can be obtained from the real ones by complexiHcation. 

These data have the following meaning. Sections of V are conHgurations of 
the gauge Held $, whereas local gauge transformations are parametrised via the 
gauge transformation operator K by sections of W. The differential operator 
P deHnes the equation of motion for the gauge Held $ via = 0. The formal 
self-adjointness of P is motivated by P being the Euler-Lagrange operator 
of a local action S'($), e.g. whereas the gauge-invariance 

condition PoK = 0 implies gauge-invariance of the action iS'(4>). This condition 
implies (for K ^ 0) that P can not be Cauchy-hyperbolic, because any ‘pure 
gauge conHguration’ = Ke with e S FoCW) of compact support solves the 
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equation of motion = 0 with vanishing initial data in the distant past, 
whereas for Cauchy-hyperbolic P the unique solution with vanishing initial 
data is identically zero. 

The Cauchy-hyperbolicity ol R = o K implies that for every $ S r(V) 
there exists an e G r(>V) such that $' = $ + itTe satisfies the ‘canonical gauge¬ 
fixing condition’ = 0. The existence of the gauge-fixing operator T such 

that the gauge-fixed equation of motion operator P = P + T o is Green- 
hyperbolic implies that every solution of = 0 in fact satisfies = 0 
up to gauge-equivalence; consequently, the dynamics of the ‘physical degrees of 
freedom’ is ruled by a hyperbolic equation of motion even if P is not hyperbolic. 
Finally, the condition that Q = o T is Cauchy-hyperbolic implies that the 
gauge-fixing = 0 is compatible with the hyperbolic dynamics of P$ = 0. 
T is in general not canonical and the following constructions will not depend 
on the particular choice of T in case several T with the required properties 
exist, thus we do not consider the gauge-fixing operator T as part of the data 
specifying the model. 

Apart from the toy model discussed above, a simple example of a linear 
gauge theory which fits into Definition 12.2.51 is the Maxwell field (on a trivial 
principal P(1 (-bundle) which after all was the inspiration for the formulation 
of this definition. This model is specified by (in differential form notation) 

}V = MxR, V = W ®T*M = T*M , 

(ei,e2)w=y ei A *62 , = J $i A *$2 , 

M M 

P = dU, K = T = d, 

P = dU + dd), K'fK = K'^T = dU = U. 

We would like to construct a (pre-)symplectic or (pre-)inner product space 
corresponding to the data given in Definition 12.2.51 by following as much as 
possible the logic of the case without gauge symmetry. To this avail we need a 
few further definitions of section spaces. 

Definition 2.2.6. As before, we denote by Sol C r(V) (Sobc C rsc(V)) the 
spaces of smooth solutions of the equation P<1> = 0 (with spacelike-compact 
support). By Q and Qsc we denote the space of gauge configurations (with 
spacelike-compact support), by t/sc,o we denote the gauge configurations induced 
by spacelike-compact gauge transformation parameters 

g = K [r(w)], g,, = go r,c(w), Ssc.o = k [rsc(w)]. 

In general, t/sc,o S Qsc- By kero(Ar^) we denote the space of gauge-invariant 
test-sections and by £ the labelling space of linear gauge-invariant on-shell 
observables 


keTo{K^) = {fGToiV)\K^f = 0} 


f = kero(K^)/P[ro(V)] . 
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Our discussion of the toy model in the previous subsection already indicated 
why £ defined above is a good candidate for a labelling space of linear gauge- 
invariant on-shell observables. First of all we observe that £ is well-defined 
because P[ro(V)] C kero(it''^) owing to P o K = 0. Moreover, we have by 
construction for arbitrary $ G Sol, e G r(>V), / G kero(it''^) and g G ro(V) 


{f + Pg,^ + Ke)^ = {f,<I>)^ . 


Consequently, every element [/] of £ induces a well-defined linear functional on 
Sol/C/, i.e. on gauge-equivalence classes of on-shell configurations, by Sol/t/ 9 
[<i>] !->■ C>[/]([<&]) = (/, <i>)y. Being gauge-invariant, these functionals correspond 
to meaningful (physical) observables. On the level of classical observables, the 
fact that the physical degrees of freedom of the gauge field propagate in a causal 
fashion is reflected in the following generalisation of Lemma [2.2.1l which is proved 
in [63] . 

Lemma 2.2.2. Let [/] € £ be arbitrary and let E be any Cauchy surface of 
{M,g). Then, for any bounded neighbourhood 0(E) of E, we can find a g G 
kero(it'^) with supp g C C’(E) and g G [/]. 

In constructing the classical bracket for models without gauge symmetry the 
last statement of Theorem 12.2.11 which in fact holds for the causal propagator 
E^ of any Green-hyperbolic operator P on an arbitrary vector bundle V, has 
been crucial. In the following, we review results obtained in [631 Theorem 
3.12+Theorem 5.2] which essentially imply that, although P is not hyperbolic, 
the causal propagator E^ of the gauge-hxed equation of motion operator P = 
P + ToiGl is effectively a causal propagator for P up to gauge-equivalence. The 
crucial observation here is that P and P coincide on kero(iir^) which implies that 
E^ restricted to kero (if is independent of the particular form of the gauge 
fixing operator T. 

Theorem 2.2.4. The causal propagator E^ of P = P + T o iff satisfies the 
following relations. 

1. h G kero(if^) and E^h G ffsc,o if o,nd only if h G P[ro(V)], with Qsc,o 
defined in Detinition \2.2.b\ 

2. Every h G Solgc can be split as h = hi + h 2 with hi G E^ [kero(if^)] and 
^2 G Gsc,0- 

3. E^ descends to a bijective map £ —)> Sobc/I/sc.o- 

f. EP is formally skew-adjoint w.r.t. (•,-)y on kero(if^), i.e. 



for all hi,h 2 G kero (if ^). 
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5. LetT' : r(>V) —>■ r(V) he any dijferential operator satisfying the properties 
required for the operator T in Definition \2.2.5\ and let E^ be the causal 
propagator of P' = P + T' o . Then E^ satisfies the four properties 
above. 

Given these results, we can now construct a meaningful bracket on £ by 
generalising Proposition 12. 2 .11 

Proposition 2.2.2. The tuple {£,t) with t : £ x £ ^ M. defined by 

Ti[f]d9]) = 

is o pre-symplectic space (Bosonic case) or pre~inner product space (Fermionic 
case). In particular, 

1. T is well-defined and independent of the chosen representatives, 

2. T is antisymmetric (Bosonic case) or symmetric (Fermionic case). 

3. Let T' : r(>V) —>■ r(V) he any differential operator satisfying the properties 
required for the operator T in Definition \2. 2. 5\ and define t' in analogy to 
T but with the causal propagator E^ of P' = P + T' o instead of E^. 
Then t' = t. 

We stress that r is in general not weakly non-degenerate, cf. the last state¬ 
ment of Theorem 12.2.51 This is a particular feature of gauge theories, cf. e.g. 
fTTUllg] for a discussion of the physical interpretation of this non-degeneracy 
in the case of the Maxwell field. The last statement above indicates that r is 
independent of the gauge-fixing operator T and in this sense, gauge-invariant. 
Indeed, we shall see in what follows that r can be rewritten in a manifestly 
gauge-invariant form. The form of r given here can be derived directly from 
the action ^ (<i>,P$)y by Peierls’ method in analogy to the derivation 

for electromagnetism in |110| . see also [751110]) for a broader context. 

As in the case without local gauge symmetries it is interesting and useful 
to observe that the covariant pre-symplectic or pre-inner product space {£, r) 
can be understood equivalently in an equal-time fashion. In fact the following 
statements have been proved in [631 Proposition 5.I-|-Theorem 5.2] (or can be 
proved by slightly generalising the arguments used there). 

Theorem 2.2.5. Under the assumptions of Detinition \2.S7^ let (■,')gQj he the 
bilinear form on Sol defined for $i, $2 G Sol with spacelike-compact overlapping 
support by 

where denotes the future part o/$, see Definition \2. 2.3\ This bilinear form 
has the following properties. 

1. (‘I*i,$2)soi well-defined for all $i,<I >2 G Sol with spacelike-compact 
overlapping support. In particular, it is independent of the choice of future 
part entering its definition. 
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2. (•,-)soi is antisymmetric (Bosonic case) or symmetric (Fermionic case). 

3. (•,-)soi is gauge-invariant, i.e. 

($l ,$2 + ife)sol = ($l,$ 2 )sol 

for all $i ,$2 G Sol, e G r(>V) s.t. <I>i and e have spacelike-compact 
overlapping .support. 

4 . (•r)soi computed as a suitable integral over an arbitrary but fixed 

Cauchy surface S of (M,g) with future-pointing normal vector field N and 
induced measure dE. If there exists a ‘current’ j : r(V) x r(V) —>■ T*M 
such that V^j^($i,$ 2 ) = -P'& 2 ))v - ((‘^ 2 ,-P‘l’i))v for all $i,<i >2 G 

r(V), then 

($i,$ 2 )soi = / dE 7V'^j^(<i>i,$2). 

5. For all $ G Sol and all h G keTo{K^), 

6. E^ descends to an isomorphism of pre-symplectic (Bosonic case) or pre- 
inner product spaces (Fermionic case) E^ : {E,t) —>■ (Solsc/^/sc, 0 : (•) •)soi)- 

'7- If Gsc,o Q Gsc, then T is degenerate, i.e. there exists [0] ^ [h] € S s.t. 
T{[h\,'£) = 0. 

The fourth and fifth statement in the above theorem show that one can view 
the observable Sol/t/ 9 [$] i -9 (/i, <f>)y ~ ^{h), i.e. the ‘covariantly smeared 
classical field’, equivalently as an ‘equal-time smeared classical field’ Sol/t/ 9 
[$] H> {H, $)g^j with H = E^h G [H] G Sobc/^/sc.o- 

2.3 Linear Quantum Fields on Curved Space- 
times 

Given the (pre-)symplectic or (pre-)inner product spaces of classical linear ob¬ 
servables constructed in the previous section for Bosonic and Fermionic theories 
with or without local gauge symmetries, there are several ‘canonical’ ways to 
construct corresponding algebras of observables in the associated quantum the¬ 
ories; these constructions differ mainly in technical terms. 

In the Bosonic case, one can consider the Weyl algebra corresponding to the 
pre-symplectic space (Sjt), which essentially means to quantize exponentials 
exp(i (/, <I>)y) of the smeared classical field $(/) — (/, $)v rather than the 
smeared classical field itself. This mainly has the technical advantage that one 
is dealing with a C'*-algebra corresponding to bounded operators on a Hilbert 
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space, i.e. to operators with a bounded spectrum. However, sometimes it is 
also advisable in physical terms to consider exponential observables rather than 
linear ones as fundamental building blocks, e.g. in case one is dealing with 
finite gauge transformations, cf. [8]. The Weyl algebra is initially constructed 
under the assumption that the form r is non-degenerate such that the the 
pre-symplectic space (£,t) is in fact symplectic. However, the construction of 
the Weyl algebra is also well-defined in the degenerate case [12] and even in 
case £ is not a vector space, but only an Abelian group [8]. In the Fermionic 
case it is not necessary to consider exponential observables in order to access 
the advantages of a C'*-algebraic framework, as a C'*-algebra can already be 
constructed based on linear observables due to the anticommutation relations 
of Fermionic quantum fields, see e.g. M- 

As in this monograph we shall not make use of the C'*-algebraic framework, 
it will be sufficient to consider the algebra of quantum observables constructed 
by directly quantizing the smeared classical fields $(/) — (/, $)y themselves 
both in the Bosonic case and in the Fermionic case. We shall first construct the 
Borchers-Uhlmann algebra A{M) corresponding to the linear model defined by 
the data {M,V,yV,P,K) cf. Definition 12.2.51 where we can consider a model 
{M,V,P) without gauge invariance as the subclass (M, V,>V = V,P,K = 0). 
The algebra A{M) contains only the most simple observables, in particular it 
does not contain observables which correspond to pointwise powers of the field 
such as or the stress-energy tensor. We shall discuss a larger algebra 

which contains also these observables in Section 12.6.11 

In order to construct the Borchers-Uhlmann algebra, we recall that the la¬ 
belling space £ consists of equivalence classes of test sections [/] corresponding to 
the classical observables $(/) — (/, with $ G [d)] G Sol/t/ (= Sol in the ab¬ 
sence of gauge symmetries, i.e. for AT = 0) and / G [/] G f = kero(Ar^)/P[ro(V)] 
(= ro(V)/P[Fo(V)] for K = 0). We recall that omitting the equivalence classes 
in the notation $(/) ~ (/, d))y is meaningful because the latter expression is 
independent of the chosen representatives of the equivalence classes. For the 
quantum theory, we need complex expressions and therefore consider the com- 
plexification £c = £ C of the labelling space £. 

With this in mind, we represent the quantum product of two smeared fields 
<l>(/i)$(/ 2 ) by the tensor product [/i] 0 [/ 2 ]- On the long run, we would like 
to represent d>(/) as an operator on a Hilbert space, we therefore need an 
operation which encodes ‘taking the adjoint with respect to a Hilbert space 
inner product’ on the abstract algebraic level. We define such a ^-operation by 
setting [$(/)]* = $(/), corresponding to [fi]* = [/i], and [$(/i) • • • $(/n)]* = 
[$(/„)]* • • • [^(/i)]* ■ Observables would then be polynomials P of smeared fields 
(tensor polynomials of elements of £c) which fulfil P* = V (in the Fermionic 
case we need further conditions e.g. P has to be an even polynomial). To 
promote $(/) to a proper quantum field with the correct (anti)commutation 
relations, we define 

= HfMg) t ^gW) = *t([/], [g])l = * (/, 5 )^ I, ( 2 . 6 ) 
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where we recall tlmt E^ is the causal propagator of the gauge-fixed equation of 
motion operator P = P + T o (= P in the absence of gauge symmetries K = 
0), I is the identity and the — (+) sign applies for the Bosonic (Fermionic) case. 
Recall that t([/], [g]) vanishes if the supports of / and g are spacelike separated, 
the above canonical (anti)commutation relations (CCR/CAR) therefore assure 
that observables commute at spacelike separations. In the case without gauge 
symmetries, one can write the CCR/CAR formally as 

[$(a;),$(y)]^ = iE^{x,y)l. 

Recall that this is nothing but the covariant version of the well-known equal¬ 
time CCR/CAR. In the case where gauge symmetries are present, this expres¬ 
sion does not make sense because the equality holds only when smeared with 
‘gauge-invariant’ test sections / S kero(Ar^). Finally, we remark that dynam¬ 
ics is already encoded by the fact that £c consists of equivalence classes with 
[Pg] = [0] G £c for all g G Fg (V) and that it is convenient to have a topology on 
the algebra A{M) in order to be able to quantify to which extent two abstract 
observables are ‘close’, i.e. similar in physical terms. We subsume the above 
discussion in the following definition. 

Definition 2.3.1. Consider a linear Bosonic or Fermionic (gauge) field theory 
defined by (M, V, W, P, AT) (with W = V and AT = 0 in the absence of gauge 
symmetries), cf. Definition 12.2.41 and Definition 12.2.51 and let {£,t) be the 
corresponding (pre-)symplectic or (pre-)inner product space constructed as in 
Theorem 12. 2. 31 and Proposition 12.2.21 The Borchers-Uhlmann algebra A{M) of 
the model (M, V, W, P, K) is defined as 

A(M) = Ao(M)/I, 

where Ao(M) is the direct sum 


Ao(M) = 0f®" 

n—0 

(£’®° = C) equipped with a product defined by the linear extension of the 
tensor product of f®”, a ^-operation defined by the antilinear extension of 
([/i] ® ® [fn])* = [fn] (D • • • (D [/i], and it is required each element of Ao{M) 

is a linear combination of elements of f®" with n < Umax < oo. Additionally, 
we equip Ao{M) with the topology induces by the locally convex topology of 
ro(V). Moreover, I is the closed *-ideal generated by elements of the form 
—*t([/], [5])©([/]C)[ff]Tb]®[/]), where — (-I-) stands for the Bosonic (Fermionic) 
case, and A{M) is thought to be equipped with the product, ^-operation, and 
topology descending from Ao{M). If O is an open subset of M, A(0) denotes 
the algebra obtained by allowing only test sections with support in O. 

A{M), in contrast to Ao{M), depends explicitly on the metric 5 of a space- 
time (M, g) via the causal propagator and the equation of motion. However, 
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now and in the following we shall omit this dependence in favour of notational 
simplicity. 

We recall that, by Lemma [2.2.21 every equivalence class [/] G Sc is so large 
that it contains elements with support in an arbitrarily small neighbourhood of 
any Cauchy surface of M. This implies the following well-known result, which 
on physical grounds entails the predictability of observables. 

Lemma 2.3.1. The Borchers-Uhlmann algebra A{M) fulfils the time-slice ax¬ 
iom. Namely, let E be a Cauchy surface of (M, g) and let O be an arbitrary 
neighbourhood of E. Then A{0) =A{M). 

We now turn our attention to states. Let 21 be a topological, unital *- 
algebra, i.e. 21 is endowed with an operation * which fulfils (AB)* = B*A* 
and (A*)* = A for all elements A,B in 21. A state w on 21 is defined to be 
a continuous linear functional 21 —>■ C which is normalised, i.e. a;(I) = 1 and 
positive, namely, uj{A*A) > 0 must hold for any A G 21. Considering the special 
topological and unital ^-algebra A{M), a state on A{M) is determined by its 
n-point correlation functions 

w„(/l, • ■■ ,fn)=UJ i^ifl) ■ ■ ■ $(/n)) • 

which are distributions in T'^{M'^). Given a state oj on 21, one can represent 
21 on a Hilbert space Hui by the so-called GNS construction (after Gel’fand, 
Naimark, and Segal), see for instance [551 H]. By this construction, algebra 
elements A G 21 are represented as operators 7r^(A) on a common dense and 
invariant subspace of Hcj, while lu is represented as a vector of G Hui such 
that for all A G 21 

a;(A) = (H^^|7rn(A)|H^^). 

Conversely, every vector in a Hilbert space H gives rise to an algebraic state on 
the algebra of linear operators on TL. 

Among the possible states on A{M) there are several special classes, which 
we collect in the following definition. Some of the definitions are sensible for 
general *-algebras, as we point out explicitly. 

Definition 2.3.2. Let 21 denote a general ^-algebra and let A{M) denote the 
Borchers-Uhlmann algebra of a linear Bosonic or Fermionic (gauge) field theory. 

1. A state w on 21 is called mixed, if it is a convex linear combination of 
states, i.e. uj = Awi + (1 — X)uj 2 , where A < 1 and uJi uj are states on 21. 
A state is called pure if it is not mixed. 

2. A state w on A{M) is called even, if it is invariant under $(/) !->■ —<!>(/), 
i.e. it has vanishing n-point functions for all odd n. 

3. An even state on A{M) is called guasifree or Gaussian if, for all even n, 

n/2 

Wn(/lj‘‘‘ I/n) = ^ ' J_J_W2 (/7r„(2i—1) j/7r„(2i)) ■ 

i—1 
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Here, 5'^ denotes the set of ordered permutations of n elements, namely, 
the following two conditions are satisfied for 7r„ S 5'^: 

7r„(2f — 1) < 7r„(2i) for 1 < i < n/2 , 

7r„(2f — 1) < 7r„(2i + 1) for 1 < f < n/2 . 

4. Let at denote a one-parameter group of ^-automorphisms on 2t, i.e. for 
arbitrary elements A, B of 21, 

at {A*B) = {at{A))* at{B), at {as{A)) = at+s{A), ao{A)=A. 

A state w on 21 is called at-invariant if uj{at{A)) = uj{A) for all A G 21. 

5. An at-invariant state w on 21 is said to satisfy the KMS condition for an 
inverse temperature /3 = T~^ > 0 if, for arbitrary elements A, B of 21, the 
two functions 

FAsit) = uj {Bat{A)) , GAB{t) = w (at(A)B) 

extend to functions Fab{z) and Gab{z) on the complex plane which are 
analytic in the strips 0 < Im z < /3 and —/3 < Im z < 0 respectively, 
continuous on the boundaries Im z G {0,/3}, and fulfil 

FAB{t + iP) = GAB{t) . 

The KMS condition (after Kubo, Martin, and Schwinger) holds naturally for 
Gibbs states of finite systems in quantum statistical mechanics, i.e. for states 
that are given as W/ 3 (A) = TvpA with a density matrix p = exp (—/3i7)(Tr exp (—/3i7)) 
H the Hamiltonian operator of the system, and Tr denoting the trace over the 
respective Hilbert space. This follows by setting 

at{A) = , 

making use of the cyclicity of the trace, and considering that exp (—/3iJ) is 
bounded and has finite trace in the case of a finite system. In the thermody¬ 
namic limit, exp (—/3i7) does not possess these properties any more, but the 
authors of [Bl] have shown that the KMS condition is still a reasonable con¬ 
dition in this infinite-volume limit. Physically, KMS states are states which 
are in (thermal) equilibrium with respect to the time evolution encoded in the 
automorphism at. In general curved spacetimes, there is no ‘time evolution’ 
which acts as an automorphism on A{M). One could be tempted to introduce a 
time evolution by a canonical time-translation with respect to some time func¬ 
tion of a globally hyperbolic spacetime. However, the causal propagator E^ 
will in general not be invariant under this time translation if the latter does 
not correspond to an isometry of {M,g). Hence, such time-translation would 
not result in an automorphism of A{M). There have been various proposals to 
overcome this problem and to define generalised notions of thermal equilibrium 
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in curved spacetimes, see m for a review. Ground states (vacuum states) may 
be thought of as KMS states with inverse temperature (3 = T~^ = oo. 

Quasifree or Gaussian states which are in addition pure are closely related to 
the well-known Fock space picture in the sense that the the GNS representation 
of a pure quasifree state is an irreducible representation on Fock space, see 
e.g. m Section 3.2]. In this sense, a pure quasifree state is in one-to-one 
correspondence to a specific definition of a ‘particle’. 

For the remainder of this chapter and most of the remainder of this mono¬ 
graph, the only model we shall discuss is the free neutral Klein-Gordon field 
for simplicity. However, all concepts we shall review will be applicable to more 
general models. 


2.4 Hadamard States 


The power of the algebraic approach lies in its ability to separate the algebraic 
relations of quantum fields from the Hilbert space representations of these rela¬ 
tions and thus in some sense to treat all possible Hilbert space representations 
at once. However, the definition of an algebraic state reviewed in the previous 
subsection is too general and thus further conditions are necessary in order to 
select the physically meaningful states among all possible ones on A{M). 

To this avail it seems reasonable to look at the situation in Minkowski space- 
time. Physically interesting states there include the Fock vacuum state and 
associated multiparticle states as well as coherent states and states describing 
thermal equilibrium situations. All these states share the same ultraviolet (UV) 
properties, i.e. the same high-energy behaviour, namely they satisfy the so- 
called Hadamard condition, which we shall review in a few moments. A closer 
look at the formulation of quantum field theory in Minkowski spacetime reveals 
that the Hadamard condition is indeed essential for the mathematical consis¬ 
tency of QFT in Minkowski spacetime, as we would like to briefly explain now. 
In the following we will only discuss real scalar fields for simplicity. Analyses of 
Hadamard states for fields of higher spin can be found e.g. in |43l [67l 1107] . 

The Borchers-Uhlmann algebra A{M) of the free neutral Klein-Gordon field 
(j) contains only very basic observables, namely, linear combinations of products 
of free fields at separate points, e.g. However, if one wants to treat 

interacting fields in perturbation theory, or the backreaction of quantum fields 
on curved spacetimes via their stress-energy tensor, ones needs a notion of 
normal ordering, i.e. a way to define field monomials like at the same 

point. To see that this requires some work, let us consider the massless scalar 
field in Minkowski spacetime. Its two point function reads 


^2{x,y) 


uj{(j){x)(j){y)) = lim 


eto dTT^ (x - yY -I- fe{xo - 2/o) + ' 


(2.7) 


where (x — yY denotes the Minkowskian product induced by the Minkowski 
metric and the limit has to be understood as being performed after integrating 
UJ 2 with at least one test function (weak limit). u} 2 {x,y) is a smooth function if 
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X and y are spacelike or timelike separated. It is singular at (x — = 0, but 

the singularity is ‘good enough’ to give a finite result when smearing uj 2 {x,y) 
with two test functions. Hence, UJ 2 is a well-defined (tempered) distribution. 
Loosely speaking, this shows once more that the product of fields (f){x)(f){y) is 
‘well-defined’ at non-null related points. However, if we were to define 'p'^{x) 
by some ‘limit’ like 

(j)'^{x) = lim (j){x)(j){y ), 

x^y 

the expectation value of the resulting object would ‘blow up’ and would not 
be any meaningful object. The well-known solution to this apparent problem 
is to define field monomials by appropriate regularising subtractions. For the 
squared field, this is achieved by setting 

■.(j)^{x): = lim [(j){x)(j){y) -W 2 (x, y)I) , 

x^y 

where of course one would have to specify in which sense the limit should be 
taken. Omitting the details of this procedure, it seems still clear that the Wick 
square :4>^(x ): is a meaningful object, as it has a sensible expectation value, i.e. 
u)[:(j}^{x):) = 0. In the standard Fock space picture, one heuristically writes the 
field (operator) in terms of creation and annihilation operators in momentum 
space, i.e. 



and defines : 4‘‘^{x) : by writing the mode expansion of the product (j){x)(j){y), 
‘normal ordering’ the appearing products of creation and annihilation operators 
such that the creation operators are standing on the left hand side of the anni¬ 
hilation operators, and then finally taking the limit a; —>■ y. It is easy to see that 
this procedure is equivalent to the above defined subtraction of the vacuum ex¬ 
pectation value. However, having defined the Wick polynomials is not enough. 
We would also like to multiply them, i.e., we would like them to constitute an 
algebra. Using the mode-expansion picture, one can straightforwardly compute 

■(t)^{x)-.:(jp-[y): = :<l)'^{x)(l)‘^{y): +4 ■.(j){x)(j){y): uj 2 {x,y) + 2[uj2{x,y)f , 

which is a special case of the well-known Wick theorem, see for instance |76] . 
The right hand side of the above equation is a sensible object if the appearing 
square of of the two-point function ui 2 {x,y) is well-defined. In more detail, 
we know that uj 2 {x,y) has singularities, and that these are integrable with test 
functions. Obviously, (a; 2 (a:,y))^ has singularities as well, and the question is 
whether the singularities are still good enough to be integrable with test func¬ 
tions. In terms of a mode decomposition, one could equivalently wonder whether 
the momentum space integrals appearing in the definition of : 4>^{x)4i^{y ): via 
normal ordering creation and annihilation operators converge in a sensible way. 
The answer to these questions is ‘yes’ because of the energy positivity prop¬ 
erty of the Minkowskian vacuum state, and this is the reason why one usually 
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never worries about whether normal ordering is well-defined in quantum field 
theory on Minkowski spacetime. In more detail, the Fourier decomposition of 
the massless two-point function uj 2 reads 

uj 2 {x, y) = lim [ dk eiko)S{k^) , (2.8) 

eto (Zny J 

where 0(fco) denotes the Heaviside step function. We see that the Fourier trans¬ 
form of ui 2 has only support on the forward lightcone (or the positive mass shell 
in the massive case); this corresponds to the fact that we have associated the 
positive frequency modes to the creation operator in the above mode expansion 
of the quantum field. This insight allows to determine (or rather, define) the 
square of u} 2 {x,y) by a convolution in Fourier space 

{uj 2 {x,y)f = lim —[ dq [dp 0(go) (5(g^) 0(po) <5(p^) uO+P'lU-y) 
eto [Zirp J J 

= hm [ dk [ dq e{qo) 6{q^) 0(fco - 9o) <5((fc - qf) . 

eto [Zirp J J 

R4 Ri 

Without going too much into details here, let us observe that the above ex¬ 
pression can only give a sensible result (a distribution) if the integral over q 
converges, i.e. if the integrand is rapidly decreasing in q. To see that this is 
the case, note that for an arbitrary but fixed k and large q where here ‘large’ 
is meant in the Euclidean norm on K.^, the integrand is vanishing on account of 
(5(g^) and 0(fco — qo) as ko — qo < 0 for large q. Loosely speaking, we observe the 
following: by the form of a convolution, the Fourier transform of UJ 2 is multiplied 
by the same Fourier transform, but with negative momentum. Since the UJ 2 has 
only Fourier support in one ‘energy direction’, namely the positive one, the in¬ 
tersection of its Fourier support and the same support evaluated with negative 
momentum is compact, and the convolution therefore well-defined. Moreover, 
as this statement only relies on the large momentum behaviour of Fourier trans¬ 
forms, it holds equally in the case of massive fields, as the mass shell approaches 
the light cone for large momenta. 

The outcome of the above considerations is the insight that, if we want to 
define a sensible generalisation of normal ordering in curved spacetimes, we have 
to select states whose two-point functions are singular, but regular enough to 
allow for pointwise multiplication. Even though general curved spacetimes are 
not translationally invariant and therefore do not allow to define a global Fourier 
transform and a related global energy positivity condition, one could think that 
this task can be achieved by some kind of a ‘local Fourier transform’ and a 
related ‘local energy positivity condition’ because only the ‘large momentum 
behaviour’ is relevant. In fact, as showed in the pioneering work of Radzikowski 
|10.311104] , this heuristic idea can be made precise in terms of microlocal analysis, 
a modern branch of Mathematics. Microlocal analysis gives a rigorous way 
to define the ‘large momentum behaviour’ of a distribution in a coordinate- 
independent manner and in the aforementioned works |10311104] , it has been 











2.4. HADAMARD STATES 


41 


shown that so-called Hadamard states, which have already been known to allow 
for a sensible renormalisation of the stress-energy tensor |1201112111123j , indeed 
fulfil a local energy positivity condition in the sense that their two-point function 
has a specific wave front set. Based on this, Brunetti, Fredenhagen, Kohler, 
Hollands, and Wald [Ill[l9ll20l[lll[72] have been able to show that one can as a 
matter of fact define normal ordering and perturbative interacting quantum field 
theories based on Hadamard states essentially in the same way as on Minkowski 
spacetimes. Though, it turned out that there is a big conceptual difference to flat 
spacetime quantum field theories, namely, new regularisation freedoms in terms 
of curvature terms appear. Although these are finitely many, and therefore lead 
to the result that theories which are perturbatively renormalisable in Minkowski 
spacetime retain this property in curved spacetimes mED], the appearance of 
this additional renormalisation freedom may have a profound impact on the 
backreaction of quantum fields on curved spacetimes, as we will discuss in the 
last chapter of this monograph. 

As we have already seen at the example of the massless Minkowski vacuum, 
Hadamard states can be approached from two angles. One way to discuss them 
is to look at the concrete realisation of their two-point function in ‘position 
space’. This treatment has lead to the insight that Hadamard states are the 
sensible starting point for the definition of a regularised stress-energy tensor 
|1201ll2lTll23] . and it is well-suited for actual calculations in particular. On the 
other hand, the rather abstract study of Hadamard states based on microlocal 
analysis is useful in order to tackle and solve conceptual problems. Following 
our discussion of the obstructions in the definition of normal ordering, we shall 
start our treatment by considering the microlocal aspects of Hadamard states. 
A standard monograph on microlocal analysis is the book of Hormander na, 
who has also contributed a large part to this field of Mathematics [lU 140] . 
Introductory treatments can be found in |1151 [TOl 15511112] . 

Let us start be introducing the notion of a wave front set. To motivate 
it, let us recall that a smooth function on R™ with compact support has a 
rapidly decreasing Fourier transform. If we take an distribution u in Fq(IR.™) 
and multiply it by an / G ro(M”’') with /(xq) ^ 0, then uf is an element of 
r'(R’”), i.e., a distribution with compact support. If fu were smooth, then its 
Fourier transform fu would be smooth and rapidly decreasing. The failure of fu 
to be smooth in a neighbourhood of xq can therefore be quantitatively described 
by the set of directions in Fourier space where fu is not rapidly decreasing. Of 
course it could happen that we choose / badly and therefore ‘cut’ some of the 
singularities of u at xq. To see the full singularity structure of u at xq, we 
therefore need to consider all test functions which are non-vanishing at Xq. 
With this in mind, one first defines the wave front set of distributions on R”’' 
and then extends it to curved manifolds in a second step. 

Definition 2.4.1. A neighbourhood F of fcg G R™ is called conic if fc G F implies 
Afc G F for all A G (0, oo). Let u G Fq(R"*). A point (xq, /cq) S R™ x (R’" \ {0}) 
is called a regular directed point of u if there is an / G ro(R™) with /(xq) 0 
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such that, for every n S N, there is a constant Cn S K. fulfilling 

\Mk)\<C„{l + \k\)-'^ 

for all fc in a conic neighbourhood of fco- The wave front set WF(u) is the 
complement in R.™ x (R™ \ {0}) of the set of all regular directed points of u. 

We immediately state a few important properties of wave front sets, the 
proofs of which can be found in [75] (see also [112] '). 

Theorem 2.4.1. Let u G rQ(R™). 

a) If u is smooth, then VIE (u) is empty. 

b) Let P be an arbitrary partial differential operator. It holds 

WF(Pu) C WF(m) . 

c) Let U, V C M™, let u G rg(ld), and let x ■ U ^ V be a diffeomorphism. 

The pull-back X*(m) of u defined by = u{x*f) for all f G Eq(U) 

fulfils 

WF(x*u) = x*WF(u) = {ix-\x),x*k) \ (x, fc) G WF(w)} , 

where x*k denotes the push-forward of x iri the sense o/cotangent vectors. 
Hence, the wave front set transforms covariantly under diffeomorphisms 
as an element o/T*R”^, and we can extend its definition to distributions 
on general curved manifolds M by patching together wave front sets in 
different coordinate patches of M. As a result, for u G Fq(M), WF(u) C 
T*M \ {0}, where 0 denotes the zero section ofT*M. 

d) Let ui, U 2 G Fq(M) and let 

WF(mi)©WF('U2 ) = {(x, ki + ^ 2 ) I (x, ki) G WF(mi), (x, ^ 2 ) G WF(m2 )} • 

//WF(iii)©WF(tt 2 ) does not intersect the zero section, then one can define 
the product U 1 U 2 in such a way that it yields a well-defined distribution 
in Fq(M) and that it reduces to the standard pointwise product of smooth 
functions if ui and U 2 are smooth. Moreover, the wave front set of such 
product is bounded in the following way 

WF(miU2) C WF(mi) U WF(u2) U (WF(ui) © WF(u2)) . 

Note that the wave front set transforms as a subset of the cotangent bundle 
on account of the covector nature of k in exp(jfcx). The last of the above 
statements is exactly the criterion for pointwise multiplication of distributions 
we have been looking for. Namely, from (12.81) and (12.71) one can infer that the 
wave front set of the Minkowskian two-point function (for m > 0) is [115] 


WF(w 2 ) = {{x,y,k, -k) G T*M^ \ x^y, {x - yf = 0, fc||(x - y), fco > 0}(2.9) 
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U {{a;, X, k, —k) G T*M^ \ k^ = 0, ko > 0} , 

particularly, it is the condition fcp > 0 which encodes the energy positivity of 
the Minkowskian vacuum state. We can now rephrase our observation that 
the pointwise square of ui 2 {x,y) is a well-defined distribution by noting that 
WF(a; 2 ) 0 WF(a; 2 ) does not contain the zero section. In contrast, we know that 
the (^-distribution d{x) is singular at a; = 0 and that its Fourier transform is a 
constant. Hence, its wave front set reads 

WF(5) = {(0,fc)|fcGR\{0}}, 

and we see that the ^-distribution does not have a ‘one-sided’ wave front set 
and, hence, can not be squared. The same holds if we view 5 as a distribution 
S{x^y) on Fo(R.^). Then 

WF((5(a;, y)) = {{x, x, k, —k) | fc G R \ {0}} . 

The previous discussion suggests that a generalisation of (12.91) to curved 
spacetimes is the sensible requirement to select states which allow for the con¬ 
struction of Wick polynomials. We shall now define such a generalisation. 

Definition 2.4.2. Let cu be a state on A{M). We say that oj fulfils the 
Hadamard condition and is therefore a Hadamard state if its two-point func¬ 
tion (jj 2 fulfils 

WF(a; 2 ) = {{x,y,kx,-ky) G T*M‘^\{0} \ {x,kx) ^ iy,ky), fca:i>0} . 

Here, {x,kx) {isky) implies that there exists a null geodesic c connecting x 
to y such that kx is coparallel and cotangent to c at a; and ky is the parallel 
transport of kx from a; to y along c. Finally, t> 0 means that the covector kx 
is future-directed. 

Having discussed the rather abstract aspect of Hadamard states, let us now 
turn to their more concrete realisations. To this avail, let us consider a geodesi- 
cally convex set O in {M,g), see Section [2T] By definition, there are open 
subsets O'x C TxM such that the exponential map exp,,. : —>• O is well- 

defined for all X & O, i.e. we can introduce Riemannian normal coordinates 
on O. For any two points x, y G O, we can therefore define the half squared 
geodesic distance a{x,y) as 

(^{x,y) = i5(exp-^(y),exp-i(y)) . 

This entity is sometimes also called Synge’s world function and is both smooth 
and symmetric on O x O. Moreover, one can show that it fulfils the following 
identity 


= 2cr, 


( 2 . 10 ) 


where the covariant derivatives are taken with respect to x (even though this 
does not matter by the symmetry of ct), see for instance [571 [im]. Let us 
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introduce a couple of standard notations related to objects on O x O such as 
a. If V and W are vector bundles over M with typical fibers constituted by the 
vector spaces V and W respectively, then we denote by VMW the exterior tensor 
product of V and W. V W is defined as the vector bundle over M x M with 
typical fibre V(^W. The more familiar notion of the tensor product bundle V®W 
is obtained by considering the pull-back bundle of V Kl W with respect to the 
map M B X {x,x) G Typical exterior product bundles are for instance 
the tangent bundles of Cartesian products of M, e.g. T*M^T*M = T*M'^. A 
section of V Kl W is called a bitensor. We introduce the Synge bracket notation 
for the coinciding point limits of a bitensor. Namely, let B be a smooth section 
of V K W. We define 


[B\{x) = lim B{x,y). 

V^x 

With this definition, [B] is a section of V (8 W. In the following, we shall 
denote by unprimed indices tensorial quantities at cc, while primed indices denote 
tensorial quantities at y. As an example, let us state the well-known Synge rule., 
proved for instance in [23l 1102) . 

Lemma 2.4.1. Let B be an arbitrary smooth bitensor. Its covariant derivatives 
at X and y are related by Synge’s rule. Namely, 

= {B]^- [B,^] ■ 

Particularly, let V be a vector bundle, let fa be a local frame of V defined on 
O <Z M and let x, y G O. If B is symmetric, i.e. the coefficients Bab'(x, y) of 

B{x,y) =B°-^\x,y) fa{x)® fbfy) 


fulfil 


then 


B‘^^\x,y) = B^'-{y,x), 

[B-A = [B.A = \[B]-,^- 


The half squared geodesic distance is a prototype of a class of bitensors of 
which we shall encounter many in the following. Namely, a fulfils a partial 
differential equation (12.101) which relates its higher order derivatives to lower 
order ones. Hence, given the initial conditions 


[cr] = 0 , [a.^] = 0 , [a.f,,,] = g^u 

which follow from the very definition of cr, one can compute the coinciding 
point limits of its higher derivatives by means of an inductive procedure, see for 
instance [371 m [Ml ins • As an example, in the case of [cr;^i,p], one differentiates 
(12.101) three times and then takes the coinciding point limit. Together with the 
already known relations, one obtains 


{^■,vvp\ — 0 ■ 
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At a level of fourth derivative, the same procedure yields a linear combination of 
three coinciding fourth derivatives, though with different index orders. To relate 
those, one has to commute derivatives to rearrange the indices in the looked- 
for fashion, and this ultimately leads to the appearance of Riemann curvature 
tensors and therefore to 

— '^i.Ti^QitT T Rfirvo) • 

A different bitensor of the abovementioned kind we shall need in the following 
is the bitensor of parallel transport g^, {x, y). Namely, given a geodesically convex 
set O, X, y € O, and a vector v = d^' in TyM, the parallel transport of v 

from y to X along the unique geodesic in O connecting x and y is given by the 
vector V in T^M with components 


This definition of the bitensor of parallel transport entails 








g^„,a.P = -a.^ . 


In fact, the first two identities can be taken as the defining partial differential 
equation of g^, and its initial condition (one can even show that the mentioned 
partial differential equation is an ordinary one). Out of these, one can obtain 
by the inductive procedure outlined above 

bp';a]=0> [9p'-cfs] = \R‘^ua0- 


With these preparations at hand, let us now provide the explicit form of Hadamard 
states. 


Definition 2.4.3. Let uj 2 be the two-point function of a state on A{M), let t 
be a time function on {M,g), let 

c^eix, y) = a(x, y) + 2ie{t{x) - t{y)) + , 


and let A be an arbitrary length scale. We say that W2 if of local Hadamard form 
if, for every xq € M there exists a geodesically convex neighbourhood O of xq 
such that a;2(a;, y) on O x O is of the form 


= Hill — 

etO StT^ 


u{x,y) 


v{x,y)log 


(^e{x,y) 

{he{x,y) + w{x,y)) . 


(Xeix,y) 

A2 


■w{x,y) 


Here, the Hadamard coefficients u, v, and w are smooth, real-valued biscalars, 
where v is given by a series expansion in a as 

OO 

U = ^ U„cr” 

n—Q 
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with smooth biscalar coefficients The bidistribution shall be called 

Hadamard parametrix, indicating that it solves the Klein-Gordon equation up 
to smooth terms. 

Note that the above series expansion of v does not necessarily converge on 
general smooth spacetimes, however, it is known to converge on analytic space- 
times [SS]. One therefore often truncates the series at a finite order n and asks 
for the w coefficient to be only of regularity C", see m- Moreover, the local 
Hadamard form is special case of the global Hadamard form defined for the first 
time in m- The definition of the global Hadamard form in m assures that 
there are no (spacelike) singularities in addition to the lightlike ones visible in 
the local form and, moreover, that the whole concept is independent of the 
chosen time function t. However, as proven by Radzikowski in |104| employing 
the microlocal version of the Hadamard condition, the local Hadamard form 
already implies the global Hadamard form on account of the fact that u }2 must 
be positive, have the causal propagator E as its antisymmetric part, and ful¬ 
fil the Klein-Gordon equation in both arguments. It is exactly this last fact 
which serves to determine the Hadamard coefficients u, u, and re by a recursive 
procedure. 

To see this, let us omit the subscript e and the scale A in the following, 
since they do not influence the result of the following calculations, and let us 
denote by P^ the Klein-Gordon operator P = —mf acting on the 
a:-variable. Applying P^ to h, we obtain potentially singular terms proportional 
to cr“" for n = 1,2,3 and to log cr, as well as smooth terms proportional to 
positive powers of a. We know, however, that the total result is smooth because 
Px{h + w) =0 since uj 2 is a bisolution of the Klein-Gordon equation and w is 
smooth. Gonsequently, the potentially singular terms in Pxh have to vanish 
identically at each order in tr and log tr. This immediately implies 

PxV = Q. ( 2 . 11 ) 

and, further, the following recursion relations 


-PxU + 2vo.^a^ + {Ux(J-2)vo = Q , ( 2 - 12 ) 

2u.^af^ + {Oxa - 4)u = 0 , (2.13) 

-PxVn + 2.{n + l)vn+i-Ii(r.+ {n+T) {Ux(J+ 2n)vn+i =1^ , Vn > 0 . (2.14) 

To solve these recursive partial differential equations, let us now focus on 
()2.13|) . Since the only derivative appearing in this equation is the derivative 
along the geodesic connecting x and y, (12.1311 is in fact an ordinary differential 
equation with respect to the affine parameter of the mentioned geodesic, u is 
therefore uniquely determined once a suitable initial condition is given. Com¬ 
paring the Hadamard form with the Minkowskian two-point function (12.7L the 
initial condition is usually chosen as 
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which leads to the well-known result that u is given by the square root of 
the so-called Van Vleck-Morette determinant, see for instance [37l f23ll56l I102j . 
Similarly, given u, the differential equation (12.121) is again an ordinary one with 
respect to the geodesic afSne parameter, and it can be immediately integrated 
since taking the coinciding point limit of (12.121) and inserting the properties of 
cr yields the initial condition 

bo] = ^[Pxu]. 

It is clear how this procedure can be iterated to obtain solutions for all 
Particularly, one obtains the initial conditions 

= 2(n+l)(n + 2)[^"'^"] 

for all n > 0. Moreover, one finds that u depends only on the local geometry of 
the spacetime, while the Vn and, hence, v depend only on the local geometry and 
the parameters appearing in the Klein-Gordon operator P, namely, the mass m 
and the coupling to the scalar curvature These observations entail that the 
state dependence of UJ 2 is encoded in the smooth biscalar w, which furthermore 
has to be symmetric because it is bound to vanish in the difference of two-point 
functions yielding the antisymmetric causal propagator E, viz. 

^ 2 {x, y) - UJ 2 {X, y) = W 2 (a;, y) - uj 2 {y, x) = iE{x, y). 

More precisely, this observation ensues from the following important result ob¬ 
tained in [89l IIO] . 

Theorem 2.4.2. The Hadamard coefficients Vn are symmetric biscalars. 

This theorem proves the folklore knowledge that the causal propagator E is 
locally given by 

iE = lim — n(h^ — h-A . 
eto Stt^ 

Even though we can in principle obtain the as unique solutions of or¬ 
dinary differential equations, we shall only need their coinciding point limits 
and coinciding points limit of their derivatives in what follows. In this respect, 
the symmetry of the Vn will prove very valuable in combination with Lemma 
12.4.11 In fact, employing the Hadamard recursion relations, we find the following 

results m- 

Lemma 2.4.2. The following identities hold for the Hadamard parametrix 
h{x,y) 




[Pxh] = [Pyh] = -6[ni], 
[{Pxh)-^i_i'] 


[{Pxh)-,n] ~ {iPv^);n'] ~ 

[(Py/i);^] = —2[ni];^ . 
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It is remarkable that these rather simple computations will be essentially 
sufficient for the construction of a conserved stress-energy tensor of a free scalar 
quantum field [91]. Particularly, the knowledge of the explicit form of, say, [ui] 
is not necessary to accomplish such a task. However, if one is interested in 
computing the actual backreaction of a scalar field on curved spacetimes, one 
needs the explicit form of [uij. One can compute this straightforwardly by the 
inductive procedure already mentioned at several occasions and the result is 
|9ll|36] 


{6^-l)m'^R {6^ - 1) R'^ (1 - 

24 288 120 

720 720 

m'^ , {6^-l)m^R , {6^ - 1) R^ , {I - 5^)aR 
“ “8~ 24 288 120 

f 

720 


(2.15) 


The Hadamard coefficients are related to the so-called DeWitt-Schwinger 
coefficients, see for instance [ssiisniins], which stem from an a priori completely 
different expansion of two-point functions. The latter have been computed for 
the first time in [24] and can also be found in many other places like, e.g. 

[SUES]. 

Having discussed the Hadamard form to a large extent, let us state the 
already anticipated equivalence result obtained by Radzikowski in |103j . See 
also [107] for a slightly different proof, which closes a gap in the proof of [103] . 

Theorem 2.4.3. Let UJ 2 be the two-point function of a state on A{M). uj 2 
fulfils the Hadamard condition of Definition 12.4-21 if and only if it is of global 
Hadamard form. 

By the result of |104j . that a state which is locally of Hadamard form is al¬ 
ready of global Hadamard form, we can safely replace ‘global’ by ‘local’ in the 
above theorem. Moreover, from the above discussion it should be clear that the 
two-point functions of two Hadamard states differ by a smooth and symmetric 
biscalar. 

In past works on (algebraic) quantum field theory in curved spacetimes, 
one has often considered only on quasifree Hadamard states. For non-quasifree 
states, a more general microlocal spectrum condition has been proposed in m 
which requires certain wave front set properties of the higher order n-point func¬ 
tions of a non-quasifree state. However, as shown in [10811109] . the Hadamard 
condition of the two-point function of a non-quasifree state alone already de¬ 
termines the singularity structure of all higher order n-point functions by the 
CCR. It is therefore sufficient to specify the singularity structure of 0 J 2 also in 
the case of non-quasifree states. Note however, that certain technical results 
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on the structure of Hadamard states have up to now only been proven for the 
quasifree case m- 

We close the general discussion of Hadamard states by providing examples 
and non-examples. 

Examples of Hadamard states 

1. Given a Hadamard state oj on A{M), any ‘finite excitation of w’ is again 
Hadamard, i.e. for all A G A{M), loa defined for all B G A{M) by 
uja{B) = uj{A*BA)/uj{A*A) is Hadamard [109j . 

2. All vacuum states and KMS (thermal equilibrium) states on ultrastatic 
spacetimes (i.e. spacetimes with a metric ds'^ = —dt^ + hijdx^dx^, with 
hij not depending on time) are Hadamard states [SSJ 1106) . 

3. Based on the previous statement, it has been proven in in |55] that Hadamard 
states exist on any globally hyperbolic spacetime by means of a spacetime 
deformation argument. 

4. The Bunch-Davies state on de Sitter spacetime is a Hadamard state 
[21 [Sni El] . It has been shown in [321 El] that this result can be generalised 
to asymptotically de Sitter spacetimes, where distinguished Hadamard 
states can be constructed by means of a holographic argument; these states 
are generalisations of the Bunch-Davies state in the sense that the afore¬ 
mentioned holographic construction yields the Bunch-Davies state in de 
Sitter spacetime. 

5. Similar holographic arguments have been used in [34l [281 E21 ES] to con¬ 
struct distinguished Hadamard states on asymptotically flat spacetimes, 
to rigorously construct the Unruh state in Schwarzschild spacetimes and 
to prove that it is Hadamard in [22] , to construct asymptotic vacuum and 
thermal equilibrium states in certain classes of Friedmann-Robertson- 
Walker spacetimes in - see m and the next chapter - and to construct 
Hadamard states in bounded regions of any globally hyperbolic spacetime 
in [32]. 

6. A interesting class of Hadamard states in general Friedmann-Robertson- 
Walker are the states of low energy constructed in |95j . These states 
minimise the energy density integrated in time with a compactly supported 
weight function and thus loosely speaking minimise the energy in the time 
interval specified by the support of the weight function. This construction 
has been generalised to encompass almost equilibrium states in [84] and 
to expanding spacetimes with less symmetry in |114) . We shall review the 
construction of these states in the next chapter. 

7. A construction of Hadamard states which is loosely related to states of 
low energy has been given in m- There the authors consider for a given 
spacetime (M, g) with the spacetime {N, g) where is a finite-time slab 
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of M. Given a smooth timelike-compact function f on M which is iden¬ 
tically 1 on N, one considers A = ifEf which can be shown to be a 
bounded and self-adjoint operator on L?{N,dgx). The positive part A+ 
of A constructed with standard functional calculus defines a two-point 
function of a quasifree state w on A{N) via W 2 (/, g) = (/, A'^g) which can 
be shown to be Hadamard (at least on classes of spacetimes) [TS] . Taking 
for / the characteristic function of N gives the Sorkin-Johnston states 
proposed in [1] which are in general not Hadamard |45j . 

8 . Hadamard states which possess an approximate local thermal interpreta¬ 
tion have been constructed in |111] . see |119] for a review. 

9. Given a Hadamard state to on the algebra A{M) and a smooth solution 

4' of the field equation P'l' = 0, one can construct a coherent state by 
redefining the quantum held (j){x) as (j){x) >->■ (j){x) + The thus 

induced coherent state has the two-point function y) = y) + 

'k(a:)'k(j/), which is Hadamard since is smooth. 

10. A construction of Hadamard states via pseudodifferential calculus was 
developed in [59] . 

Non-examples of Hadamard states 

1. The so-called a-vacua in de Sitter spacetime [2] violate the Hadamard 
condition as shown in |20j . 

2 . As already mentioned the Sorkin-Johnston states proposed in [T] are in 
general not Hadamard [^ . 

3. A class of states related to Hadamard states, but in general not Hadamard, 
is constituted by adiabatic states. These have been introduced in [96] and 
put on rigorous grounds by [86] . Effectively, they are states which approx¬ 
imate ground states if the curvature of the background spacetime is only 
slowly varying. In [77], the concept of adiabatic states has been gener¬ 
alised to arbitrary curved spacetimes. There, it has also been displayed in 
a quantitative way how adiabatic states are related to Hadamard states. 
Namely, an adiabatic state of a specihc order n has a certain Sobolev wave 
front set (in contrast to the C°° wave front set introduced above) and 
hence, loosely speaking, it differs from a Hadamard state by a biscalar of 
finite regularity G". In this sense, Hadamard states are adiabatic states 
of ‘infinite order’. We will review the concept of adiabatic states in the 
next chapter 

Finally, let us remark that one can define the Hadamard form also in space- 
times with dimensions differing from 4, see for instance [TOTII^. Moreover, 
the proof of the equivalence of the concrete Hadamard form and the microlocal 
Hadamard condition also holds in arbitrary spacetime dimensions, as shown in 
[I 03 . A recent detailed exposition of Hadamard states may be found in [81] . 
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2.5 Locality and General Covariance 

And important aspect of QFT on curved spacetimes is the backreaction of quan¬ 
tum fields in curved spacetimes, i.e. the effect of quantum matter-energy on 
the curvature of spacetime. This of course necessitates the ability to define 
quantum field theory on a curved spacetime without knowing the curved space- 
time beforehand. It is therefore advisable to employ only generic properties of 
spacetimes in the construction of quantum fields. This entails that we have 
to formulate a quantum field theory in a local way, i.e. only employing local 
properties of the underlying curved manifold. In addition, we would like to take 
into account the diffeomorphism-invariance of General Relativity and therefore 
construct covariant quantum fields. This concept of a locally covariant quantum 
field theory goes back to many works, of which the first one could mention is 
[5H] . followed by many others such as |123l Chapter 4.6] and [1181170] . Building 
on these works, the authors of m have given the first complete definition of a 
locally covariant quantum field theory. 

As shown in ED, giving such a definition in precise mathematical terms re¬ 
quires the language of category theory, a branch of mathematics which basically 
aims to unify all mathematical structures into one coherent picture. A category 
is essentially a class £ of objects denoted by o6j(£), with the property that, for 
each two objects A, B m.€. there is (at least) one morphism or arrow : A ^ B 
relating A and B. The collection of all such arrows is denoted by hom(r{A, B). 
Morphisms relating a chain of three objects are required to be associative with 
respect to compositions, and one demands that each object has an identity 
morphism id^ : A —^ A which leaves all morphisms (j) : A ^ B starting from A 
invariant upon composition, i.e. (j) o idA = 4>- An often cited simple example of 
a category is the category of sets 6et. The objects of 6et are sets, while the 
morphism are maps between sets, the identity morphism of an object just being 
the identity map of a set. Given two categories £i and £ 2 , a, functor F : £1 —>■ £2 
is a map between two categories which maps objects to objects and morphisms 
to morphisms such that identity morphisms in £1 are mapped to identity mor¬ 
phism in £2 and the composition of morphisms is preserved under the mapping. 
This paragraph was only a very brief introduction to category theory and we 
refer the reader to the standard monograph [HD and to the introduction in [II3L 
Section 1.7] for further details. A locally covariant quantum field theory ac¬ 
cording to ED should be a functor from a category of spacetimes to a category 
of suitable algebras. The first step in understanding such a construction if of 
course the definition of a suitable category of spacetimes. 

We have already explained in the previous sections of this chapter that four¬ 
dimensional, oriented and time-oriented, globally hyperbolic spacetimes are the 
physically sensible class of spacetimes among all curved Lorentzian manifolds. 
It is therefore natural to take them as the objects of a potential category of 
spacetimes. Regarding the morphisms, one could think of various possibilities 
to select them among all possible maps between the spacetimes under consid¬ 
eration. However, to be able to emphasise the local nature of a quantum field 
theory, we shall take embeddings between spacetimes. This will allow us to 
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require locality by asking that a quantum field theory on a ‘small’ spacetime 
can be easily embedded into a larger spacetime without ‘knowing anything’ 
about the remainder of the larger spacetime. Moreover, a sensible quantum 
field theory will depend on the orientation and time-orientation and the causal 
structure of the underlying manifold, we should therefore only consider embed¬ 
dings that preserve these structures. To this avail, the authors in m have 
chosen isometric embeddings with causally convex range (see Section 12.11 re¬ 
garding an explanation of these notions), but since the causal structure of a 
spacetime is left invariant by conformal transformations, one could also choose 
conformal embeddings, as done in [98]. We will nevertheless follow the choice of 
[21] . since it will be sufficient for our purposes. Let us now subsume the above 
considerations in a definition. 

Definition 2.5.1. The category of spacetimes IBlan is the category having as its 
class of objects obj{'0Jlan) the globally hyperbolic, four-dimensional, oriented and 
time-oriented spacetimes {M,g). Given two spacetimes {Mi,gi) and (M 2 , 32 ) 
in obj{ifflan), the considered morphisms hormjyian ((Mi, 51 ), (M 2 , 32 )) are isomet¬ 
ric embeddings x ■ (Mi,gi) ^ {M^^g^) preserving the orientation and time- 
orientation and having causally convex range x(Mi). Moreover, the identity 
morphism id(^M,g) of a spacetime in obj(fOlan) is just the identity map of M and 
the composition of morphisms is defined as the usual composition of embed¬ 
dings. 

The just defined category is sufficient to discuss locally covariant Bosonic 
quantum field theories. However, for Fermionic quantum field theories, one 
needs a category which incorporates spin structures as defined in [11811108] . At 
this point we briefly remark that our usage of the words ‘Boson’ and ‘Fermion’ 
for integer and half-integer spin fields respectively is allowed on account of the 
spin-statistics theorem in curved spacetimes proved in [118j . see also [41] for a 
more recent and general work. 

To introduce the notion of a locally covariant quantum field theory and 
the related concept of a locally covariant quantum field, we need a few cate¬ 
gories in addition to the one introduced above. By TSllg we denote the cate¬ 
gory of unital topological *-algebras, where for two Ai, A 2 in obj(T2l[g), the 
considered morphisms hom.xaig(^ 1 , ^ 2 ) are continuous, unit-preserving, injec¬ 
tive *-homomorphisms. In addition, we introduce the category Test of test 
function spaces ro(M) of objects {M,g) in 9Jlan, where here the morphisms 
homxe 5 t(ro(Mi),ro(M 2 )) are push-forwards y, of the isometric embeddings 
X : Ml ^ M 2 . In fact, by V we shall denote the functor between 91lan and 
Test which assigns to a spacetime {M,g) in IDlan its test function space ro(M) 
and to a morphism in 911an its push-forward. For reasons of nomenclature, we 
consider T2llg and Test as subcategories of the category Top of all topological 
spaces with morphisms given by continuous maps. Let us now state the first 
promised definition. 

Definition 2.5.2. A locally covariant quantum field theory is a (covariant) func¬ 
tor A between the two categories 91tan and T2llg. Namely, let us denote by 
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the mapping A{x) of a morphism x in to a morphism in T 2 llg and by 

A{M, g) the mapping of an object in SOfon to an object in T2l[g, see the follow¬ 
ing diagram. 

^ (-^2,92) 

■"1 ^ 

A{Mi,gi) -^ A{M2,g2) 

Then, the following relations hold for all morphisms Xij G horngjianCC-Tfi, 5 i), {Mj,gj)) 

0x23 ° 0x12 = 0x230X12 J ^idM ~ ■ 

A locally covariant quantum field theory is called causal if in all cases where 
Xi S h.omman{{Mi, gi), {M, g)) are such that the sets xi{AIi) and X 2 {^ 2 ) are 
spacelike separated in {M,g), 

[axiiA{Mi,gi)),a^^{A{M2,g2))] = { 0 } 

in the sense that all elements in the two considered algebras are mutually com¬ 
muting. 

Finally, one says that a locally covariant quantum field theory fulfils the 
time-slice axiom, if the situation that y G homOTon((Mi, gi), (M 2 , 52 )) is such 
that x(Mi,gi) contains a Cauchy surface of (M 2 , 52 ) entails 

Ox iA{Mi,gi)) = A{M2,g2) ■ 

The authors of [ 2 T] also give the definition of a state space of a locally 
covariant quantum field theory and this turns out to be dual to a functor, by 
the duality relation between states and algebras. One therefore chooses the 
notation of covariant functor for a functor in the strict sense, and calls such a 
mentioned dual object a contravariant functor. We stress once more that the 
term ‘local’ refers to the size of spacetime regions. A locally covariant quantum 
field theory is such that it can be constructed on arbitrarily small (causally 
convex) spacetime regions without having any information on the remainder 
of the spacetime. In more detail, this means that the algebraic relations of 
observables in such small region are already fully determined by the information 
on this region alone. This follows by application of the above definition to the 
special case that (Mi, 51 ) is a causally convex subset of (M 2 , 52 )- 

As shown in [21] , the quantum field theory given by assigning the Borchers- 
Uhlmann algebra A{M) of the free Klein-Gordon field to a spacetime {M,g) 
is a locally covariant quantum field theory fulfilling causality and the time- 
slice axiom. This follows from the fact that the construction of A{M) only 
employs compactly supported test functions and the causal propagator E. The 
latter is uniquely given on any globally hyperbolic spacetime, particularly, the 
causal propagator on a causally convex subset (Mi, 51 ) of a globally hyperbolic 
spacetime (M 2 , 32 ) coincides with the restriction of the same propagator on 
(M 2 , 52 ) to (Ml, 51 ). Finally, causality follows by the causal support properties 
of the causal propagator, and the time-slice axiom follows by Lemma [2. 3. II 
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Let us now discuss the notion of a locally covariant quantum field. These 
fields are particular observables in a locally covariant quantum field theory which 
transform covariantly, i.e. loosely speaking, as a tensor and are in addition 
constructed only out of local geometric data. In categorical terms, this means 
that they are natural transformations between the functors D and A. We refer 
to m for the notion of a natural transformation, however, its meaning in our 
context should be clear from the following definition. 

Definition 2.5.3. A locally covariant quantum field $ is a natural transforma¬ 
tion between the functors V and A. Namely, for every object {M,g) in 9Jlan 
there exists a morphism $(M.g) : ^o{M, g) — ^ A{M, g) in Top such that, for each 
morphism y S homOTon((Mi, gi), (Mi, 51 )), the following diagram commutes. 


ro(Mi, 5 i) A(Mi,gi) 



ro(Af2,5'2) -^ A{M2,g2) 

*(M2,32) 


Particularly, this entails that 


o ^(Mi,gi) - ^(Ms.gs) O X* ■ 


It is easy to see that the Klein-Gordon field (j){f) is locally covariant. Namely, 
the remarks on the local covariance of the quantum field theory given by A{M) 
after Definition l2.5.2l entail that an isometric embedding y : (Mi, 51 ) ^ (M 2 , 52 ) 
transforms (f>{f) as 

ttx ('/'(/)) = <(’(X*/), 


or, formally. 


Q^x = </>(x(a;)) ■ 


Hence, local covariance of the Klein-Gordon field entails that it transforms as 
a ‘scalar’. While the locality and covariance of the Klein-Gordon field itself are 
somehow automatic, one has to take care that all extended quantities, like Wick 
powers and time-ordered products thereof, maintain these good properties. The 
prevalent paradigm in algebraic quantum field theory is that all pointlike ob¬ 
servables should be theoretically modelled by locally covariant quantum fields. 

A comprehensive review of further aspects and results of locally covariant 
quantum field theory may be found in [4^. 


2.6 The Quantum Stress—Energy Tensor and the 
Semiclassical Einstein equation 

The aim of this section is to discuss the semiclassical Einstein equation and the 
quantum-stress-energy tensor : : which is the observable whose expectation 

value enters this equation. As argued in the previous section, all pointlike 
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observables such as the quantum stress-energy tensor should be locally covariant 
fields in the sense of Definition [2331 Rather than discussing local covariance for 
non-linear pointlike observables only at the example of : it is instructive 

to review the construction of general local and covariant Wick polynomials. 

2.6.1 Local and Covariant Wick Polynomials 

The first construction of local and covariant general Wick polynomials was given 
in m based on ideas already implemented for the stress-energy tensor in |123] . 

Here we would like to review a variant of the construction of m in the spirit 
of the functional approach to perturbative QFT on curved spacetimes, termed 
perturbative algebraic quantum field theory, cf. nziisiiisa. Essentially, this 
point of view on local Wick polynomials was already taken in m- We review 
here only the case of the neutral Klein-Gordon field, however, the functional 
approach is applicable to general field theories 11051 HU [52] ■ 

In the functional approach to algebraic QFT on curved spacetimes, one 
considers observables as functionals on the classical field configurations. Upon 
quantization, these functionals are endowed with a particular non-commutative 
product which encodes the commutation relations of quantum observables. We 
have already taken this point of view in the discussion of the Borchers-Uhhnann 
algebra A{M) as the result of quantizing a classical symplectic space constructed 
in Section [2. 2 1 In particular, the smeared quantum scalar field 4>{f) consid¬ 
ered as the quantization of the linear functional Ff : F(M) —?> C, Ff{(j)) = (/, (f) 
with / £ Fg (M). The new aspect in the approach we shall review now is to con¬ 
sider a much larger class of functionals on T{M). To this avail, we view T{M) 
as the space of off-shell configurations of the scalar field, whereas Sol C F(M) is 
the space of on-shell configurations. For the purpose of perturbation theory it 
is more convenient to perform all constructions off-shell first and to go on-shell 
only in the end, and we shall follow this route as well, even though perturbative 
constructions are not dealt with in this monograph. 

To this end, we call a functional F : T{M) —>■ C smooth if the n-th functional 
derivatives 



exist for all n and all ipi,... ,tpn € ro(M) and if £ r'(M"), i.e. 

is a distribution. By definition F^^'^ [cf) is symmetric and we consider only 
polynomial functionals, i.e. {f) vanishes for n > TV and some TV. We define 
the support of a functional as 


supp F = {x £ M\'i neighbourhoods U oi x 3(f), if) £ TfM), supp if <ZU, 


such that Ffcf) + if)) ^ P{4’)} ■ 


(2.17) 
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which coincides with the union of the supports of over all (j) € r(M). 

The relevant space of functionals which encompasses all observables of the free 
neutral scalar field is the space of microcausal functionals 

= {-f" • r(M) —)> C I F smooth, compactly supported, 

WF n (yl u F") = 0 }, (2.18) 

where V+j- is a subset of the cotangent space formed by the elements whose 
covectors are contained in the future/past light cones and denotes is 

closure. contains two subspaces of importance. On the one hand, it contains 
the space T\oc of local functionals consisting of sums of functionals of the form 

FW = J dgx V[(j)]g,^...g,^{x) 

M 

where G T{T*M^) is a tensor such that is a (tensor) product of 

covariant derivatives of (p at the point x with a total order of n and / G Fg(TM") 
is a test tensor. The prime example of a local functional is a smeared field 
monomial 

= J dgX Pixffix) cp\f) , / G rC(M). (2.19) 

M 

One the other hand, contains the space Ereg of regular functionals, i.e. 
all microcausal functionals whose functional derivatives are smooth such that 
G for all p and all n. A prime example of a regular functional 

is a functional of the form 


F{4>) = \{{fj,4>) , GrC(M). (2.20) 

Linear functionals are the only functionals which are both local and regular. 

Given a bidistribution H G r(p(M^) which a) satisfies the Hadamard condi¬ 
tion Definition l2.4.2l b) has the antisymmetric part H{x, y) — H{y, x) = iE(x, y) 
defined by the causal propagator E and a real symmetric part 

Hsym{x,y) = i {H{x,y) +H{y,x)) , 

and c) is a bisolution of the Klein-Gordon equation P„.H{x, y) = PyH{x, y) = 0, 
we define a product indicated by -ku on via 

00 

{FkHCm = ^ - (</>),(</>)), ( 2 . 21 ) 

n—0 

i.e. by the sum of all possible mutual contractions of F and G by means of 
H. This is just an elegant way to implement Wick’s theorem by an algebraic 
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product as we shall explain in the following and may be understood in terms 
of deformation quantization, cf. |53| for a review of this aspect. Owing to the 
Hadamard property of H and the microlocal properties of microcausal function¬ 
als, the *//-product is well-defined on by m Theorem 8.2.13] and one can 
show that F -kff G £ F^c for all F,G £ J'^c, i-e. is closed under 

It is not necessary to require that H is positive and thus the two-point func¬ 
tion UJ2 of a Hadamard state oj on A{M). However, independent of whether or 
not we require H to be positive, its real and symmetric part is not uniquely fixed 
by the conditions we imposed. Given a different H' satisfying these conditions, 
it follows that d = H' — H = H!.y^ — i?sym is real, symmetric and smooth and 
that the product *h’ is related to -kn by 

F kH' G = ad {a-d{F) k^ a-d{H)) , (2.22) 


where aa '. F^c —t J'qic is the ‘contraction exponential operator’ 



(2.23) 


The previous discussion implies that, given a H satisfying the above-mentioned 
conditions, we can define a meaningful off-shell algebra W%{M) = {F^cAh)^ 
and a corresponding on-shell algebra WniM) = bV)?f(M)/I, where I is the ideal 
generated by Jyjc.soi C Fp,c, the microcausal functionals which vanish on on-shell 
configurations (f £ Sol C r(M). In fact, we have Wh{M) = (.F^c/-^/ic,soi,*ff) 
and in the following we shall indicate an equivalence class [T’j £ F^c/Ffj,c,soi by 
F for simplicity. 

As a further implication of the previous exposition we have that Wjj (M) and 
constructed with a different H of the required type are isomorphic via 
(kd ■ y^H[YI) —>• Wh'{M) with d = H' — H. In this sense, we can consider 
Wh{M) as a particular representation of an abstract algebra >V(M) which is 
independent of H, cf. [Hj. Wh{M) is in fact a ^-algebra and ad ■ WniM) — >■ 
Wh'{AI) is a ^-isomorphism if we define the involution (^-operation) on 
via complex conjugation by 

f*(</>) = fM 


which implies 


{FkH GY =G* kH F* 


by the conditions imposed on H. WniM) may be endowed with a topology in¬ 
duced by the so-called Hdrmander topology [BUS], and one can show that all 
continuous states on Wh{M) are induced by Hadamard states on the Borchers- 
Uhlmann algebra A{M), cf. [55] (in combination with |109j b In fact, we shall 
now explain why Wh{M) may be considered as the ‘algebra of Wick polynomi¬ 
als’ and in particular as an extension of A{M). 

To this avail, we first consider two linear functionals = Fij.((f>), i = 

1,2, cf. (I2.19|) . i.e. the classical field smeared with /i ,/2 £ Tq{M). The 
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definition of the ^//-product (|2.21l) implies 

mh), = Hfl) </-(/ 2 ) - </>(/ 2 ) Hfl) = . 

This indicates that the product -kn encodes the correct commutation relations 
among quantum observables. If we consider instead the quadratic local func¬ 
tionals = F 2 j-{(j)), i = 1,2, cf. (I2.19|) . i.e. the pointwise square of the 

classical field smeared with /i, /2 £ Tg (M), we find 

0"(/2)=0"(/l)</>"(/2) 

+ 4 f dgxdgyd{x,y) (l){x)(l){y)H{x,y)fi{x)f2iy) + 2H^{fi,f2), 

M2 

which we may formally write as 

(l)‘^{x) -kH (l>'^{y) = (l)‘^ix)(l)‘^{y) + mx)(l){y)H{x,y) +2H‘^{x,y). 

This expression may be compared by the expression obtain via Wick’s theorem 
:(l)^{x):H--4>^{y)--H=-(f^{x)(j)^{y):H +4 :0(x)0(y) -.h H{x,y) V 2H'^{x,y)l, 

if we define ■.-'.h to be the Wick-ordering w.r.t. the symmetric part i?sym of H, 

e.g. 

■.(j)[x)(j){y)-.H= (t>{x)(j){y) - Hsy^{x,y)l, :(j)^{x):H= hm :(j){x)(j){y) :h ■ 

x^y 

Consequently, local functionals F G J^oc considered as elements of F £ 
correspond to Wick polynomials Wick-ordered with respect to Hsy^, formally 
one may write :F:h= ct-HsyiniE) with defined as in (12.231) . In fact, the 

algebra F £ Wh{M) contains also time-ordered products of Wick polynomials 
[laizQiizi] and the perturbative construction of QFT based on Wh{M) implies 
that ‘tadpoles’ are already removed. 

We have anticipated that Wh{M) is an extension of the Borchers-Uhlmann 
algebra A{M), which contains only products of the quantum field (j) at different 
points. To see this, we consider the algebra ^'(M) = a-H^y^ ((•^reg/-^reg.soi, *//)), 
where Freg,soi — ^reg ^ ^c,sob (J'reg/J'reg,soi,*ff) IS a subalgebra of Wff(M) = 
(•^/jc/-^/jc,soi, ) because -k^ closes on regular functionals and we note that 

i® well-defined on J^eg although Hsym is not smooth. One can check 
that A'{M) is in fact the algebra A'{M) = {J-reg/Eyeg,soiAE) with the product 
■kE given by (j2.21l) with H replaced by iE/2, and that A'{M) is isomorphic to 
the Borchers-Uhlmann algebra A{M) defined in Definition l2.3.1l 

We have mentioned that Hadamard states oj on A{M) induce meaningful 
states on WniM) and that in fact all reasonable states on Wh{M) are of this 
form. To explain this in more detail, we consider a Gaussian Hadamard state on 
u) on A{M) with two-point function C02 and the algebra yVuj 2 {M) constructed by 
means of uj2- Given this, we can define a Gaussian Hadamard state on 

by 

u:{F) = = 0), V F £ 
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which corresponds to the fact that Wick polynomials Wick-ordered w.r.t. to UJ 2 
have vanishing expectation values in the state oj, e.g. u}{:(jP{x) = 0. Note 

that this definition implies in particular that 

UJ {(t>{fl)-ku:2 (t>{f2)) = W2(/l,/2), 

i.e. UJ 2 is, as required by consistency, the two-point correlation function of w also 
in the functional picture. If we prefer to consider the extended algebra Wh{M) 
constructed with a different i/, then we can define the state u) on Wh{M) by 
a pull-back with respect to the isomorphism ad ■ Wh{YI) 
d = (jJ 2 — H and ad as in (I2.23I1 . In other words, lo o ad dehnes a state on 
Wh{M), and this definition corresponds to e.g. 

= \m\{uj 2 {x,y) - H{x,y )), 

x^y 


i.e. to a point-splitting renormalisation of the expectation value of the observ¬ 
able (j)'^{x). 

We recall that observable quantities should be local and covariant fields as 
discussed in the previous section. However, not all local elements of the algebra 
Wh{M) satisfy this property, i.e. not all local functionals F S Jjoc considered 
as elements of yV//(M) correspond to local and covariant Wick polynomials. 
In particular the functional ^^(/) = i^ 2 ,/(^) — : -h, cf. (I2.19|l . does not 

correspond to a local and covariant Wick-square because H is by assumption 
a bisolution of the Klein-Gordon equation and thus : -h does not only 

depend on the geometry, i.e. the metric and its derivatives, in the localisation 
region of the test function /, but also on the geometry of the spacetime {M,g) 
outside of the support of / m- This is related to the observation that quite 
generally local and covariant Hadamard states do not exist [U] . Notwithstand¬ 
ing, local and covariant elements of Wh{M) do exist and, following [701 17^ - 
they can identified by means of a map Wh ■ iT\oc —t -^loc C Wh{M) which 
should satisfy a number of conditions: 

1. Wh commutes with functional derivatives, i.e. (ITf/(F))(^) = Wh{F^^'^) 
and with the involution Wh{F)* = Wh(F*). 

2. Wh satisfies the Leibniz rule. 

3. Wh is local and covariant. 

4. Wh scales almost homogeneously with respect to constant rescalings m i—>■ 
\m and g i—^ ^~^g of the mass m and metric g. 

5. Wh depends smoothly or analytically on the metric g, the mass m and the 
coupling ^ to the scalar curvature present in the Klein-Gordon equation 

(I2a . 

We refer to m for a detailed discussion of these conditions, and only sketch 
their meaning and physical motivation. To this avail, we consider the smeared 
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local polynomials = Fkj{(j)) defined in (12.191) and identify Wh{(I>^[x)) 

by : (jA{x) : omitting the smearing function for simplicity. The first of the 
above axioms then imply that :,(j){y)\ = imE{x,y) :(j)‘^~^{x) i.e. the 

Wick-ordered fields satisfy standard commutation relations. The Leibniz rule 
further demands that ■.(j)^{x) :=: i.e. Wick-ordering commutes 

with covariant derivatives, and the locality and covariance condition demands 
that : (jA{x) : is a local and covariant field in the sense of Definition 12.5.31 
Finally, the scaling condition requires that under constant rescalings m i—>■ \m 
and g i—?> X~^g, :(jA(x ): scales (in four spacetime dimensions) as ■.(jA{x) :i—?> : 

(jA{x): +0(log A), which among other things implies that: (jA{x): has the correct 
‘mass dimension’, and the smoothness / analyticity requirement implies that e.g. 
■.(jV{x ): may not contain a term like e.g. exp(^“^)m“'^i?“^(a:)(i?^y 
which would be allowed by the previous conditions. 

It has been demonstrated in [70l [72] (see also [91]) that a prescription of 
defining local and covariant Wick polynomials exists, but that this prescription 
is not unique. In fact, if we consider (in a geodesically convex neighbourhood) a 
H of the form (12.4.31) with w = 0, i.e. a purely geometric Hadamard parametrix, 
then Wick-ordering w.r.t. to this H, e.g. : (j)^{x) :=: (j)^[x) :h= a-{ x)) 
satisfies all conditions reviewed above. However, this prescription is not the 
only possibility, but one can consider e.g. 

: (j)'^ (x): = : (x) :h +aR{x) + jSrri^ 

with arbitrary real and dimensionless constants a and j3 which are analytic 
functions of These constants parametrise the renormalisation freedom of 
Wick polynomials, or, in the context of perturbation theory, the renormalisation 
freedom inherent in removing tadpoles. Note that a change of scale A in (I2.4.3|) 
can be subsumed in this renormalisation freedom as a particular one-parameter 
family. 

The coefficients parametrising the renormalisation freedom of local Wick 
poly-nomials may not be fixed within QFT on curved spacetimes, but have to 
be determined in a more general framework or by comparison with experiments. 
We will comment further on this point when discussing the renormalisation free¬ 
dom of the stress-energy tensor in the context of cosmology in the next chapter. 
Note that, by local covariance, these coefficients are universal and may be fixed 
once and for all in all globally hyperbolic spacetimes (of the same dimension). 
Admittedly, in view of the above presentation of locality and general covariance, 
one might think that this holds only for spacetimes with isometric subregions 
(or spacetimes with conformally related subregions on account of |98j 1. How¬ 
ever, given two spacetimes (Mi, 51 ) and (M 2 , 52 ) with not necessarily isometric 
subregions, one can employ the deformation argument of [54] to deform, say, 
(Ml, 51 ) such that it contains a subregion isometric to a subregion of {M 2 ,g 2 )- 
As the renormalisation freedom is parametrised by constants multiplying cur¬ 
vature terms or dimensionful constants which maintain their form under such 
a deformation, one can require that the mentioned constants are the same on 
(Mi,(/i) and (M 2 , 32 ) in a meaningful way. 
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2.6.2 The Semiclassical Einstein Equation and Wald’s Ax¬ 
ioms 

The central equation in describing the influence of quantum fields on the back¬ 
ground spacetime - i.e. their backreaction - is the semiclassical Einstein equa¬ 
tion. It reads 


G^,{x)=%ttGuj{-.T^,{x):), (2.24) 

where the left hand side is given by the standard Einstein tensor — 

\R9tiui G denotes Newton’s gravitational constant, and we have replaced the 
stress-energy tensor of classical matter by the expectation value of a suitable 
Wick polynomial '.T^vix ): representing the quantum stress-energy tensor eval¬ 
uated in a state lo. Considerable work has been invested in analysing how such 
equation can be derived via a suitable semiclassical limit from some potential 
quantum theory of gravity. We refer the reader to [IHl Section II.B] for a review 
of several arguments and only briefly mention that a possibility to derive (I2.24[) 
is constituted by starting from the sum of the Einstein-Hilbert aetion <S'eh( 5) 
and the matter action S'matter(<7, ‘k), 

S{g, $) = SEnig) - S'matter(g, ^') , (2.25) 

SEH(.g) = j dgX R= j dx^J\deig\ R 

M M 

formally expanding a quantum metric and a quantum matter field around a 
classical (background) vacuum solution of Einstein’s equation, and computing 
the equation of motion for the expected metric while keeping only ‘tree-level’ 
{h^) contributions of the quantum metric and ‘loop-level’ (fi^) contributions of 
the quantum matter field. In this work, we shall not contemplate on whether 
and in which situation the above mentioned ‘partial one-loop approximation’ is 
sensible, but we shall take the following pragmatic point of view: (I2.24|l seems 
to be the simplest possibility to couple the background curvature to the stress- 
energy of a quantum field in a non-trivial way. We shall therefore consider (I2.24p 
as it stands and only discuss for which quantum states and Wick polynomial 
definitions of : Tg^{x) : it is a self-consistent equation. A rigorous proof that 
solutions of the semiclassical Einstein equation actually exist in the restricted 
case of cosmological spacetimes has been given in 1101) . 

We first observe that in (12.241) one equates a ‘sharp’ classical quantity on 
the left hand side with a ‘probabilistic’ quantum quantity on the right hand 
side. The semiclassical Einstein equation can therefore only be sensible if the 
fluctuations of the stress-energy tensor : (x): in the considered state w are 

small. In this respect, we already know that we should consider w to be a 
Hadamard state and '.Tg^ix ): to be a Wick polynomial Wick-ordered by means 
of a Hadamard bidistribution. Namely, the discussion in the previous sections 
tells us that this setup at least assures finite fluctuations of : Tg^{x) : as the 
pointwise products appearing in the computation of such fluctuations are well- 
defined distributions once their Hadamard property is assumed. In fact, this 
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observation has been the main motivation to consider Hadamard states in the 
first place |12011123] . However, it seems one can a priori not obtain more than 
these qualitative observations, and that quantitative statements on the actual 
size of the fluctuations can only be made a posteriori once a solution of (12.241) 
is found. An extended framework where the left hand side of the semiclassical 
Einstein equation is also interpreted stochastically is discussed in [lOOj . 

Having agreed to consider only Hadamard states and Wick polynomials con¬ 
structed by the procedures outlined in the last section, two questions remain. 
Which Hadamard state and which Wick polynomial should one choose to com¬ 
pute the right hand side of (12.2411 ? The first question can ultimately only be 
answered by actually solving the semiclassical Einstein equation. Observe that 
this actually poses a non-trivial problem as the formulation of the semiclassi¬ 
cal Einstein equation in principle requires to specify a map which assigns to a 
metric g a Hadamard state Wg, whereas we know that a covariant assignment of 
Hadamard states to spacetimes does not exist [44]. This problem can be par¬ 
tially overcome by defining such a map only on a particular subset of globally 
hyperbolic spacetimes as we shall see in Section [3.2.21 The question of which 
Wick polynomials should be taken as the definition of a quantum stress-energy 
tensor is also non-trivial, as Wick-ordering turns out to be ambiguous in curved 
spacetimes, see [Toiizg and the last section. We have already pointed out at 
several occasions that one should define the Wick polynomial : T^y(a;) : in a 
local, and, hence, state-independent way. In the context of the semiclassical 
Einstein equation the reason for this is the simple observation that one would 
like to solve (12.241) without knowing the spacetime which results from this pro¬ 
cedure beforehand, but a state solves the equation of motion and, hence, already 
‘knows’ the full spacetime, thus being a highly non-local object. On account 
of the above considerations, we can therefore answer the question for the cor¬ 
rect Wick polynomial representing : (cc): without having to choose a specific 

Hadamard state ui beforehand. The following review of the quantum stress- 
energy tensor will be limited to the case of the free neutral Klein-Gordon field. 
An analysis of the case of Dirac fields from the perspective of algebraic QFT on 
curved spacetimes may be found in |26| . whereas the case of interacting scalar 
helds in treated in m- 

To this avail, let us consider the stress-energy tensor of classical matter helds. 
Given a classical action S'matter, the related (Hilbert) stress-energy tensor can 
be computed as fmiis] 


T — 

fii/ — 


6S„ 




(2.26) 


For the Klein-Gordon action S{g,(j)) = ^ {(j), P(j)) we hnd 


Tp.v =(1 - 


(2.27) 
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A straightforward computation shows that the classical stress-energy tensor 
is covariantly conserved on-shell, i.e. 

= 0. 

Moreover, a computation of its trace yields 

= (6C - 1) (</>□</- + . 

Particularly, we see that in the conformally invariant situation, that is, m = 0 
and ■? = gj the classical stress-energy tensor has vanishing trace on-shell. In 
fact, one can show that this is a general result, namely, the trace of a classical 
stress-energy tensor is vanishing on-shell if and only if the respective field is 
conformally invariant |491 Theorem 5.1]. 

In view of the discussion of local and covariant Wick polynomials in the 
previous section, it seems natural to define the quantum stress-energy tensor 
: T^v{x) : just as the Wick polynomial : T^v{x) :h Wick-ordered with respect 
to a purely geometric Hadamard parametrix H = of the form (12.4.31) with 
w = 0. Given a Hadamard state uj whose two-point function is (locally) of the 
form UJ 2 = H + the expectation value of : T^^(x) :h in this state may be 
computed as 


{■.T^^{x) :h) = lim 


D'ilS'{.x,y)w{x,y) 

87r2 


87r2 


where the bidifferential operator may be obtained from the classical stress- 
energy tensor as 




1 

2 S(l){x)6(l){y) 


(1 - 205];'- 2CV^V, + 

+ |2en. + (^2^ - 0 g;'vPVp> - , 


and where we recall the Synge bracket notation for coinciding point limits of 
bitensors, cf. Section [211 Here, unprimed indices denote covariant derivatives 
at X, primed indices denote covariant derivatives at y and g]) is the bitensor 
of parallel transport, cf. Section 12.41 However, as we shall discuss in a bit 
more detail in the following, this particular definition of a quantum stress- 
energy tensor is not satisfactory because it does not yield a covariantly conserved 
quantity which is a necessary condition for the semiclassical Einstein equation 
to be well-defined. 

The first treatment of the quantum stress-energy tensor from an algebraic 
point of view was the analysis of Wald in [120] . At the time [120] appeared, work¬ 
ers in the field had computed the expectation value of the quantum stress-energy 
tensor by different renormalisation methods like adiabatic subtraction, dimen¬ 
sional regularisation, ((-function regularisation and DeWitt-Schwinger point¬ 
splitting regularisation (see [T3] and also [5H1|S2]) and differing results had been 
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found. From the rather modern point of view we have reviewed in the previ¬ 
ous sections, it is quite natural and unavoidable that renormalisation in curved 
spacetimes is ambiguous. However, the axioms for the (expectation value of) 
the quantum stress-energy tensor introduced in |120) helped to clarify the case 
at that time and to understand that in principle all employed renormalisation 
schemes were correct in physical terms. Consequently, the apparent differences 
between them could be understood on clear conceptual grounds. These axioms 
(in the updated form presented in |123j l are: 

1. Given two (not necessarily Hadamard) states uj and to' such that the dif¬ 
ference of their two-point functions UI 2 — 0 J 2 is smooth, uj{:T^l,{x) :) — w'(: 
T^uix ):) is equal to 

{^2{x,y) - uj'2{x,y))\ . 

2. uj{:Tf^^{x ):) is locally covariant in the following sense: let 

X-.{M,g)^{M',g') 

be dehned as in Section [23] and let denote the associated continuous, 
unit-preserving, injective *-morphisms between the relevant enlarged (ab¬ 
stract) algebras W{M,g) and W[M',g'). If two states w on W{M,g) and 
oj' on VV(M', g') respectively are related via uj = uj' o then 

uj' {:T^>i,'{x'):) = x*a; (:T^,.(a;):) , 

where x» denotes the push-forward of x in the sense of covariant tensors. 

3. Covariant conservation holds, i.e. 

V^a;(:T^,(a:):) = 0. 

4. In Minkowski spacetime M, and in the relevant Minkowski vacuum state 
Wm 

wm {■■T^y{x)-.) = 0. 

5. uj{:Tfj^L,{x ):) does not contain derivatives of the metric of order higher than 

2 . 

Some of these axioms are just special cases of the axioms for local Wick 
polynomials reviewed in Section 12.6. II In fact, the first of these axioms is just 
a variant of the requirement that local Wick polynomials have standard com¬ 
mutation relations. In the case of : : this implies that two valid definitions 

of this observable can only differ by a term proportional to the identity. The 
second axiom is just the locality and covariance of Wick polynomials here formu¬ 
lated on the level of expectation values. The condition that the quantum-stress 
energy tensor has vanishing expectation value on Minkowski spacetime in the 
corresponding vacuum state is not compatible with the requirement that local 
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Wick polynomials depend smoothly on the mass m, cf. [TO Theorem 2.1 (e)] 
and can be omitted because it is not essential. The last axiom is motivated by 
the wish to ensure that the solution theory of the semiclassical Einstein equation 
does not depart ‘too much’ from the one of the classical Einstein equation. In 
particular one would like to have that all solutions of the semiclassical Einstein 
equation behave well in the classical limit h —?> 0. Wald himself had realised, 
however, that this axiom could not be satisfied for massless theories without 
introducing an artificial length scale into the theory; therefore, the axiom has 
been discarded. 

Using these axioms as well as a variant of the scaling requirement for local 
Wick polynomials, Wald could prove that a uniqueness result for uj{-.Tf^y :) (and 
thus for : T^i,: itself by the first axiom) can be obtained, namely that two valid 
definitions : T ^^,: and : T^i ,of the quantum stress-energy tensor can only differ 
by a term of the form 

■ . . T^^if. — 0^1777. T o;2777 G T o.'^I^u E 0 . 4 . Jiiif T cK^if , (2.28) 

where ai are real and dimensionless constants and the last three tensors appear¬ 
ing above are the conserved local curvature tensors 

M 

J dgX RapR ^ 

M 

+ DR) + - DR^, + 2R^pR^f, , (2.29) 

J ^9^ Ra/SjSR 

M 

'^RatiR°‘i, — YARfj^y -|- 2R.fj^i,. 

This uniqueness result follows by using the first of Wald’s axioms in order to 
observe that : T^i/: is a c-number. From the locality and covariance 
axiom in combination with conservation it follows that this c-number must be a 
local and conserved curvature tensor whereas the scaling condition implies that 
it has the correct mass dimension. As m pointed out, these requirements do 
not fix : : to be of the above form, but demanding in addition 

that : Tg^i, : depends in a smooth or analytic way on m and g as in [zniiTi] is 
sufficient to rule out the additional terms mentioned in m so that (12.281) indeed 
classifies the full renormalisation freedom compatible with the axioms of local 
Wick polynomials introduced in [701171] and the conservation of 

In fact, we will see in the next section that changing the scale A in the 
regularising Hadamard bidistribution amounts to changing : : exactly by 


dgX R — ~g^v ( 'zR + 2ni? ) -|- 2R.g^ii -f 2RR^ 


dfiiy — 


J — 


I 


vTdet^ 

I 

vTdet^ 




= 


vTdet^ 
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a tensor of this form and, furthermore, the attempt to renormalise perturba¬ 
tive Einstein-Hilbert quantum gravity at one loop order automatically yields a 
renormalisation freedom in form of such a tensor as well [ni0. Moreover, using 
the Gauss-Bonnet-Chern theorem in four dimensions, which states that 



M 


is a topological invariant and, therefore, has a vanishing functional derivative 
with respect to the metric laine], one can restrict the freedom even further 
by removing from the list of allowed local curvature tensors as = 
— IfjLu- Finally, the above tensors all have a trace proportional to Di? and 
thus the linear combination I^^, — is traceless. 

2.6.3 A Conserved Quantum Stress—Energy Tensor 

After some dispute about computational mistakes (see the discussion in [121) 1 it 
had soon be realised that the quantum stress-tensor : T^i, -.h defined by Wick¬ 
ordering the canonical expression by means of a purely geometric Hadamard 
bidistribution H is not conserved although it satisfies the other conditions for 
local Wick polynomials mentioned at the end of Section 12.6.11 The reason 
for this is the fact that, in contrast to the two-point function of a Hadamard 
state, a purely geometric Hadamard bidistribution fails to satisfy the equation of 
motion. Consequently ://=: :h is (in four spacetime dimensions) 

the covariant divergence of a non-vanishing and non-conserved local curvature 
term ■T^y:H= ^ 0. The obvious solution to this problem has been 

to compute and then to define a conserved stress-tensor by 


(2.30) 



i.e. by just subtracting this conservation anomaly. It was found that this con¬ 
served stress tensor still has a trace anomaly, i.e. the trace of any stress-energy 
tensor which satisfies the Wald axioms is non-vanishing in the conformally in¬ 
variant case TO = 0, ^ = g [121j . 

Two proposals have been made in order to motivate the ad-hoc subtrac¬ 
tion (12.301) on conceptual grounds. In [21] it has been suggested that the 
Wick-ordering prescription should be kept, but that the classical expression 
T^j/ entering the definition of : T^^, \h should be modified in such a way that 
it coincides with the canonical classical stress-energy tensor on-shell but gives 
a conserved observable upon quantization. In m instead it has been argued 
that the Wick-ordering prescription should be modified without changing the 
classical expression. This point of view has the advantage that it fits into the 
general framework of defining local and covariant Wick polynomials introduced 

^In fact, at least in the case of scalar fields, the combination of the local curvature tensors 
appearing as the finite renormalisation freedom in ESI is, up to a term which seems to be an 
artifact of the dimensional regularisation employed in that paper, the same that one gets via 
changing the scale in the regularising Hadamard bidistribution. 
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in [701 [72]. In particular it is possible to alter uniformly the definition of all 
local Wick polynomials induced by Wick-ordering with respect to a purely geo¬ 
metric Hadamard parametrix H in such a way that a) all axioms for local Wick 
polynomials mentioned at the end of Section 12.6.11 are still satisfied and b) the 
canonical stress-energy tensor Wick-ordered with respect to this prescription 
is conserved In the aforementioned reference it has also been argued that 
this point of view has the further advantage that it is applicable even to pertur- 
batively constructed interacting models. Notwithstanding, we shall review the 
approach of [91] in the following for ease of presentation. 

To this avail, we modify the classical stress-energy tensor by setting 

where is the canonical expression and c is a suitable constant to be fixed 
later. We then define : : by 


■ fiv ■ — ■ fj,v -H J 


where : • :h indicates Wick-ordering (i.e. point-splitting regularisation) with 
respect to a purely geometric H = of the form (12.4.31) with w = Q. Following 
the arguments of the previous section, the expectation value of : T^i, : in a 
Hadamard state w whose two-point function is (locally) of the form u :2 = H+-^ 
may be computed as 




87r2 ’ 


(2.31) 


where 


= DIT + (2.32) 

= (1 - 205"'- 2eV^V, -b 

+ 0 gN^^P' ~ + ^9^.i^Px ■ 

The following result can now be shown m- 

Theorem 2.6.1. Let be defined as in (I2.31|) with c = |. 

a) uj{:Tfj,^(x):) is covariantly conserved, i.e. 

V'^w(:T^4x):) =0. 


b) The trace of uj(;.T^,j{x)\) equals 

g^'^uj i-.T^Ax) 0 = ^ bi] - i (^ (^ - ^) ° + 

= 2^ -1) + 6(1 - 50ni? + - y ) 
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1 (rrA (6^ — 1) m^R 

'4^ V~8~ 24 


1 


+ ™ It"']’ 


which, for m = 0 and f = ^, constitutes the trace anomaly of the quantum 
stress-energy tensor. 

c) The conservation and trace anomaly are independent of the chosen scale 
A in the Hadamard parametrix h. Namely, a change 


results in 


U}{:Tf_i„{x):) uj{:Tf^^{x):)' = uj{:Tf^^{x):) + STfj,,^ , 


where 


^ 21ogA/A' . 

- s_2 


Stt^ 




21ogA/A' — l)G^i/ m 


Stt^ 


Stt^ 

1 


12 


'''' I /T Q 7 ^ 


(2.33) 


(6^-1)% 

144 


is a conserved tensor which has vanishing trace for m = 0 and C = g • 

Proof. 1. Leaving c undetermined and employing Synge’s rule (cf. Lemma 
[2XT]) . we compute 


87r2v'‘cc(:r^,(x):)=V'‘ w] = )D 


in' 


(^gf V„'Px + cigl Vu'Px + VvPx)^ w 

Let us now recall that Px{h pw) =0 and, hence, PxW = —Pxh. Inserting 
this and the identities found in Lemma 12.4.21 we obtain 

(:T^.(x):) = - [(-<?/'V.'P. + + V,P,)) h 

= (6c-2) [vi].^ . 

This proves the conservation for c = ^. 

2. It is instructive to leave c undetermined also in this case. Employing 
Synge’s rule and the results of Lemma [2.4.21 we find 

{■.T^,,{x) :) = w] 


= (4c - l)[Pxw\ - (3(--^)n + TO )[«;] 


= (4c - l)6[ui] - 1^3 1^- - ^ ) □ + TO )[«;]. 
Inserting c = -f yields the wanted result. 
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3. The proof ensues without explicitly computing 6 T^i, in terms of the stated 
conserved tensors from the following observation. Namely, a change of 
scale as considered transforms ic by a adding a term 21ogA/A'u. Hence, 
our computations for proving the first two statements entail 


2 log A/A' 




~9v ^f'Px + ^v'Px + VjyPi 


= -(4c - i)[P.ci - (3 (i -«) □ + M. 

The former term vanishes because PxV = 0 as discussed in Section 12.41 
and the same holds for the latter term if we insert ? = g and m = 0. 

□ 


The proof the second statement clearly shows that there is a possibility to as¬ 
sure vanishing trace in the conformally invariant case, but this possibility is not 
compatible with conservation. Moreover, we have stated the last result in ex¬ 
plicit terms in order to show how a change of scale in the Hadamard parametrix 
is compatible with the renormalisation freedom of the quantum stress-energy 
tensor. Finally, we stress that the term added to the canonical stress-energy 
tensor is compatible with local covariance because the corresponding change of 
the quantum stress-energy tensor is proportional to gij,i^[vi], i.e. a local curva¬ 
ture tensor. 

We would also like to point out that the above explicit form of the trace 
anomaly has also been known before Hadamard point-splitting had been devel¬ 
oped. Particularly, the same result had been obtained by means of so-called 
De Witt-Schwinger point-splitting in [23] . This renormalisation prescription is a 
priori not rigorously defined on Lorentzian spacetimes and the Hadamard point¬ 
splitting computation in |121) had therefore been the first rigorous derivation 
of the trace anomaly of the stress-energy tensor. However, DeWitt-Schwinger 
point-splitting can be reformulated on rigorous grounds, cf. [62) . 


2.7 Further Reading 

The review of algebraic quantum field theory on curved spacetimes in this chap¬ 
ter has covered aspects of this framework which are relevant for the applications 
discussed in the following chapter. Recent reviews which deal with aspects and 
constructions not covered in the present chapter, or provide further details, are 
[SU [73] and [ini IZl nil ns]) which are part of (TB]. A historical account of 
quantum field theory in curved spacetimes may be found in [124] . 
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Chapter 3 


Cosmological Applications 


Abstract. In this chapter we discuss two cosmological applications of algebraic 
quantum field theory in curved spacetimes. In the Standard Model of Cosmology 
- the ACDM-model - the matter-energy content of the universe on large scales is 
modelled by a classical stress-energy tensor of perfect fluid form. Motivated by 
the fact that this matter-energy is considered to have a microscopic description 
in terms of a quantum field theory, we demonstrate as a first application how 
the classical perfect fluid stress-energy tensor in the the ACDM-model may be 
derived within quantum field theory on curved spacetimes by showing that there 
exist quantum states on cosmological spacetimes in which the expectation value 
of the quantum stress-energy tensor is qualitatively and quantitatively of the 
form assumed in the ACDM-model up to corrections which may have interesting 
phenomenological implications. In the simplest models of Inflation, it is assumed 
that a classical scalar field on a cosmological spacetime coupled to the metric 
via the Einstein equations drives an exponential phase of expansion in the early 
universe. As a second application, the standard approach to the quantization 
of the perturbations of this coupled system, which makes heavy use of the 
symmetries of cosmological spacetimes, is re-examined by comparing it with 
a more fundamental approach which consists of quantizing the perturbations 
of a scalar field and the metric field in a gauge-invariant manner on general 
backgrounds and then considering the symmetric cosmological backgrounds as 
a special case. 


3.1 A Brief Introduction to Cosmology 

According to the well-known cosmological principle, our universe is homogeneous 
and isotropic. This postulate implies that, on large scales, the cosmos looks ‘the 
same’ everywhere and in all directions, see [621 Chapter 5] for a precise definition 
and a discussion of these issues. A remarkable confirmation of the isotropy of our 
universe is the fact that the temperature of the Cosmic Microwave Background 
(CMB) is isotropic up to relative fluctuations of the order 10“® [1]. Based 
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on the cosmological principle, we shall regard Friedmann-Lemaitre-Robertson- 
Walker (FLRW) spacetimes as the curved manifolds describing our universe on 
large scales, where here ‘large’ means scales of the size of galaxy superclusters, 
i.e. around 10® light-years or lO^^m, which is just about a thousandth of the 
diameter of the observable universe! The underlying manifold of such spacetimes 
is It where R denotes an open interval in R and is a three-dimensional 

manifold of constant curvature k, and their metric is given by the line element 


ds^ = —dv 


\t) 


dr^ 


1 — Kr^ 


Here, a(t) is a strictly positive smooth function called the scale factor, t £ R 
denotes cosmological time, and 0,ip are coordinates on the 2-sphere the 
canonical line element of which is denoted by If we recall the discussion in 
Section I^Tl we immediately realise that (or rather {t} x E„ for all t £ R) is a 
Cauchy surface, and hence all FLRW spacetimes are globally hyperbolic. In the 
following, we shall restrict attention to k = 0 as this is the situation strongly 
favoured by experimental data, the CMB measurements in particular [2]. In 
this case, one speaks of flat FLRW spacetimes and E^ = K®, while r denotes 
the Euclidean distance in R®, which in the cosmological context is called the 
comoving distance. 

Obviously, a possible time function on a FLRW spacetime is given by the 
cosmological time itself, but there are further possibilities, which are often more 
convenient: the conformal time t], the scale factor a itself, and the redshift z. 
These time variables are related by 


dt = adr] 


da 


dz 

{l + z)H^ 


(3.1) 




r] = 



H=^ 


ao 1 
z =-1 


where H denotes the Hubble rate and to and oq are conventional, as is the 
interval Ir^ ^ rj determined by R and the above integral expression for r]. In 
cosmology, oq is interpreted as the scale factor of today and usually chosen as 
oo = 1 so that the present redshift is z = 0. The conformal time is always a well- 
defined time variable, whereas using a and z as time variables is only meaningful 
if the Hubble rate H has a definite sign. Observations indicate H > 0, i.e. 
an expanding universe, and we will make this assumption when discussing the 
cosmological evolution. In fact, as the name suggests, z is a useful time variable 
in cosmology because it is actually a direct observable that can be measured by 
comparing the observed spectra of distant luminous objects with known spectra 
and assigning the ‘time-label’ z to the time when the light observed today was 
emitted from these objects. From the theoretical point of view, using z as a 
time variable has the advantage that - as we shall discuss in the following - 
the (semiclassical) Einstein equation in the simplest cases becomes an ordinary 
differential equation for H{z), whereas by using t or rj as variables one has to 
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deal with differential equations where both H and a appear as functions of t or 
rj. The latter variable is called conformal time because the metric of flat FLRW 
spacetimes in the Cartesian coordinate system ( 77 , x) S x reads 

= 0 ( 77 )^ {—drf' + dx^') , 

which indicates that flat FLRW spacetimes are conformally flat. The conformal 
time is a convenient time variable because one can eliminate terms which are 
first order in dt in the Klein-Gordon equation on FLRW by passing to 77 as 
we shall see. In the following, we shall use ' and ' to indicate derivatives with 
respect to t and 77 respectively. 

The functional behaviour of the scale factor a describes the ‘history’ of our 
universe, which, according to General Relativity, is completely determined by 
the specification of the matter-energy content of our universe in terms of the 
stress-energy tensor and its coupling to gravity via the Einstein equation 

= SttGT^, . (3.2) 

In order to obtain flat FLRW spacetimes as solutions of this equation, must 
have the have the form of a stress-energy tensor of a perfect fluid, namely, in 
comoving coordinates (t,x), 

= diag(p,p,p,p). (3.3) 

In (|3.3I) . the energy density g and the pressure p of matter-energy are related 
by the equation of state 


P = p{p) =w{q)q, (3.4) 

where the last form of the equation of state is convenient because in cosmology 
w{q) is often a constant. 

In fact, according to the Standard Model of Gosmology - the AGDM-model 
- our universe contains matter, radiation, and Dark Energy, whose combined 
energy density determines the expansion of the universe. In the AGDM-model, 
these three kinds of matter-energy are modelled macroscopically by perfect flu¬ 
ids with equations of state p = wp, w = 0, — 1 for matter, radiation and Dark 

Energy (assuming that the latter is just due to a cosmological constant), respec¬ 
tively. In the context of cosmology, the terms ‘matter’ and ‘radiation’ subsume 
all matter-energy with the respective macroscopic equation of state such that 
e.g. ‘radiation’ does not comprise only electromagnetic radiation, but also the 
three left-handed neutrinos present in Standard Model of Particle Physics (SM) 
and possibly so-called Dark Radiation, and ‘matter’ includes both the baryonic 
matter which is well-understood in the SM and Dark Matter. Here, Dark Mat¬ 
ter and Dark Radiation both quantify contributions to the macroscopic matter 
and radiation energy densities which exceed the ones expected from our knowl¬ 
edge of the SM and are believed to originate either from particles respectively 
fields not present in the SM or from geometric effects such as modifications of 
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General Relativity. Observations indicate that the current matter-energy con¬ 
tent of the universe is composed of roughly 30% matter and 70% Dark Energy, 
while the relative contribution of radiation is orrly of order 10“"^, see e.g. [2] 
for the exact numbers from the 2013 data release of the Planck satellite. The 
measured w of Dark Energy is in good agreement with a constant w = — 1 [2], 
where in this context Dark Energy is mostly just taken to be all matter-energy 
which does not behave macroscopically like matter or radiation. 

In flat FLRW spacetimes, the Einstein equation reduces to the Friedmann 
equations 




8ttG 


3 


9, 


d 

a 


IttG 

~3“ 


{g + 3p). 


(3.5) 


The Einstein equation is only consistent if the stress-energy tensor is co- 
variantly conserved = 0. In (not necessarily flat) FLRW spacetimes, the 

covariant conservation of the stress-energy tensor implies 


g + 3 H{q + p) = 0 , 


(3.6) 


and this equation can be obtained directly from the Friedmann equations since 
= 0 has been implicitly assumed in their derivation. Alternatively, given 
conservation of , it is sufficient to solve only the first Friedmann equation in 

( 1 ^ . 

Under the assumption that interactions between the individual matter- 
energy components are negligible on large scales - which is well-motivated for 
redshifts z < 10®, cf. Section fl.2.11 -. the stress-energy tensor of each matter- 
energy component is conserved on its own and (1^ implies for the energy 
densities of radiation, matter and a cosmological constant A respectively 


Prad _ ^rad 

go 


Qmat _ ^mat 

^0 


— = Da. 
£*0 


(3.7) 


Here, by convention, go is the energy density of today and thus the constants 
Drad — 10“^, Dmat — 0.3, Da — 0.7 quantify the relative contributions of ra¬ 
diation, matter and the cosmological constant to the present energy density. 
Clearly, matter, radiation, and the cosmological constant have very different 
scaling behaviours with respect to a. The first Friedmann equation implies 
that, if ^3 > 0 for all times and a > 0 at one instant of time, then a will 
be strictly increasing for all times. Consequently, if we consider the present 
matter-energy content described above and assume that Dark Energy is a cos¬ 
mological constant, then our universe must have had two phases of evolution 
preceding the present era dominated by Dark Energy: a phase where radiation 
has determined the behaviour of a followed by a matter-dominated era. This 
motivates examining the solutions of the Friedmann equations separately for 
each matter-energy component and one finds 


arad(0 oc 


amat(i) (X {t-to)N 


aA{t) oc evT*. 


(3.8) 
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The outcome of the preceding discussion is that, under the mentioned as¬ 
sumptions, our universe must have inevitably faced a Big Bang at some point of 
time in the past, i.e. there has been a to > —oo with a(to) = 0 and we shall set 
to = 0 in the following. Note that the occurrence of a Big Bang follows already 
from the second Friedmann equation and the assumptions a > 0, p > 0, since 
then g + 3p > 0 and therefore d < 0 if we take into account the sum of all three 
matter-energy constituents in the radiation-dominated era at early times. 

The Big Bang scenario is known to have (at least) one hitch, usually termed 
horizon problem. We refer to e.g. [MlIlTlIlI] for a quantitative discussion of 
this issue and only consider its qualitative aspects here. To wit, the isotropy of 
the temperature of the CMB radiation entails that the so-called last scattering 
surface, i.e. the region from where the CMB photons we see today have been 
emitted, must have lied in the forward lightcone of some event responsible for 
the thermal equilibrium of such region. The size of the last scattering surface is 
the radius rem of our past lightcone at the time tem of CMB photon emission, 
namely, the speed of light times the conformal time difference pit now) — vitem)- 
The isotropy of the CMB therefore entails that the following inequality must 


hold 


em 



0 


and one can compute that this is not the case in the standard model of cosmol¬ 
ogy; this is precisely the horizon problem. A prominent possibility to solve the 
horizon problem is Inflation, see for instance dn m mi [ST]. In this scenario, 
one usually assumes that, in the very early universe, there has been an addi¬ 
tional matter-energy component mimicking a large cosmological constant and 
thus leading to phase of exponential expansion. Inserting this assumption into 
(13.91) leads to a large negative p{0) and therefore allows for (13.91) to be fulfilled. 
In the simplest models of Inflation the matter-energy component responsible 
for the exponential expansion is a classical scalar field with a suitable potential. 


3.2 The Cosmological Expansion in QFT on Curved 
Spacetimes 


Following the program outlined in Section 11.2.11 we shall now explain how the 
post-BBN cosmological evolution can be understood within QFT on curved 
spacetimes. In particular, we shall argue that there exist phenomenologically 
well-motivated Hadamard states for free quantum fields on FLRW spacetimes 
which approximately solve the semiclassical Einstein equation to a good degree 
in such a way that the energy density in these states and on the spacetime 
provided by the solution of the semiclassical equation is qualitatively and quan¬ 
titatively of the form assumed in the ACDM-model up to small corrections 
whose possible interpretation we shall discuss. As argued in Section 11.2.11 we 
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shall make the simplified assumption that matter and radiation are microscop¬ 
ically modelled by a pair of conformally coupled free neutral scalar fields which 
are massive and massless respectively. 

3.2.1 The Renormalisation Freedom of the Quantum Stress- 
Energy Tensor in the Context of Cosmology 

We start our analysis by analysing the renormalisation freedom of the quantum 
stress-energy tensor in the cosmological context. To this avail, we recall the 
discussion in Section [TH We consider a Hadamard state oj whose two-point 
function uj 2 is (locally) of the form uj 2 {x,y) = H{x,y) + where H is 

the purely geometric, singular and state-independent part and w encodes the 
smooth and state-dependent part in (12.4.31) . The expectation value of '■T^^{x): 
in the state w reads 

u) ;) = ^ oiiui T 02 m -t- -t- , (3.10) 

where is the bidifferential operator 


= (1 - 2^91^- 2yy, + 



the bracket [•] denotes the coinciding point limit, I^ and are conserved local 
curvature tensors (j2.29|) and ai are real dimensionless constants which are an¬ 
alytic functions of the coupling to the scalar curvature One should think of 
the Qfi as encoding the combined renormalisation freedom of all quantum matter 
fields. In our model, only the massive scalar field ‘contributes to’ ai and a 2 , 
whereas 03 and 04 encode the combined freedom of the massive and massless 
conformally coupled scalar fields. 

The parameters Ui are free parameters of the theory which are independent 
of the spacetime (M, g) and can in principle be fixed by experiment, just like the 
mass m. The general physical interpretation of the occurrence of these a priori 
undetermined parameters is as follows. In usual particle physics experiments we 
always measure the difference of the expectation value of : T^^ : in two states, 
e.g. the vacuum and a many-particle state. However, gravity is sensitive to 
the absolute value of w (:T^iy:), thus the unambiguous specification of w (:T^iy:) 
would correspond to a specification of a ‘zero point’ in the absolute energy scale, 
but this is impossible within quantum field theory in curved spacetime. 

ai and 02 can be interpreted as a renormalisation of the cosmological con¬ 
stant and a renormalisation of Newton’s constant, respectively. In the following 
we will take the point of view that 02 is not a free parameter because Newton’s 
constant has been measured already. In order to do this, we have to fix a value 
for the length scale A in the Hadamard parametrix H (12.4.31) as the Definition 
of the smooth part w and thus the expression [D^^w] depend on A, cf. The¬ 
orem 12.6.11 and in particular (12.3311 . We do this by confining A to be a scale 
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in the range in which the strength of gravity has been measured. Because of 
the smallness of the Planck length, the actual value of A in this range does not 
matter as changing A in this interval gives a negligible contribution to w :). 
Moreover, in the case of conformal coupling which we shall assume most 

of the time, a 2 is independent of A as one can infer from (12.331) . One could 
also take a more conservative point of view and consider a 2 to be a free pa¬ 
rameter, in this case comparison with cosmological data, e.g. from Big Bang 
Nucleosynthesis, would presumably constrain a 2 to be very small once A is in 
the discussed range. Omitting the freedom parametrised by q; 2 , we are left with 
three free parameters in the definition of uj {'.T^u :). One of them corresponds to 
the cosmological constant which is already a free parameter in classical General 
Relativity, whereas the other two parameters do not appear there and thus will 
by themselves lead to an extension of the ACDM-model. 


At this point, we would like to highlight the point of view on the so-called 
cosmological constant problem taken in this work, as well as in most works on 
algebraic QFT on curved spacetimes and e.g. the review [^. It is often said that 
QFT predicts a value for the cosmological constant which is way too large in 
comparison to the one measured. This conclusion is reached by computing one 
or several contributions to the vacuum energy density in Minkowski spacetime 
Avac and finding them all to be too large, such that, at best, a fine-tuned 
subtraction in terms of a negative bare cosmological constant Abare is necessary 
in order to obtain the small value Avac + Abare we observe. In this work, we 
assume as already mentioned the point of view that it is not possible to provide 
an absolute definition of energy density within QFT on curved spacetimes, and 
thus neither Avac nor Abare have any physical meaning by themselves; only 
Avac + Abare Is physical and measurable and any cancellation which happens in 
this sum is purely mathematical. The fact that the magnitude of Avac depends 
on the way it is computed, e.g. the loop or perturbation order, cf. e.g. [53], 
is considered to be unnatural following the usual intuition from QFT on flat 
spacetime. However, it seems more convincing to us to accept that Avac and 
Abare have no relevance on their own, which does not lead to any contradiction 
between theory and observations, rather than the opposite. 


In the recent work m it is argued that a partial and unambiguous relevance 
can be attributed to Avac by demanding Abare to be analytic in all coupling con¬ 
stants, as required by the axioms for local and covariant Wick polynomials 
reviewed in Section [2.6.1l taking this point of view, one could give the contribu¬ 
tion to Avac which is non-analytic in these constants an unambiguous meaning. 
Indeed the authors of [27j compute a non-perturbative and hence non-analytic 
contribution to Avac which turns out to be small. In the view of this, one could 
reformulate our statement in the above paragraph and say that contributions 
to Avac and Abare which are analytic in masses and coupling constants have no 
physical relevance on their own. 
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3.2.2 States of Interest in Cosmological Spacetimes 

We shall now discuss Hadamard states of interest in flat FLRW spacetimes. As 
argued in Section [1.2.11 in the context of describing the cosmological expansion 
within QFT on curved spacetimes we are interested in Hadamard states which 
may be phenomenologically interpreted as generalised thermal excitations of 
generalised vacuum states. 

The spatial translation and rotation invariance of flat FLRW spacetimes 
allows to give a Fourier decomposition of solutions of the Klein-Gordon equation 
in terms of modes. To this end, we recall that the metric g of a flat FLRW 
spacetime is conformally related to the Minkowski metric go via g = a^go- This 
implies that the Klein-Gordon operator on FLRW spacetimes can be written as 

P = -□ + ^ a , 


where V denotes the gradient with respect to the (comoving) spatial coordinates. 
We thus see that scalar field in a FLRW spacetime is conformally equivalent to 
a scalar field with time-varying mass in Minkowski spacetimes. This leads us 
to define a mode solution of P(j) = 0 as 




(27r)ia(77) ’ 


(3.11) 


where k = \k\ and the temporal mode Xk{v) is a solution of the ordinary differ¬ 
ential equation 


d; + k- 


+ a{r]f 


^ “ 6 ) Riv) + aivym- 


Xkiv) = 0, 


(3.12) 


which depends on k only via k. For each fc, the solution space of this equation 
is two-dimensional and can be parametrised without loss of generality by any 
solution Xkiv) and its complex conjugate Xkiv) which satisfies the normalisation 
condition 


Xk{v)d^Xk{r]) - Xk{v)d,^Xk{v) = i ■ (3.13) 

Note that this condition is well-defined because the Wronskian on the right 
hand side of (j3.13l) is constant in g on account of (13.1211 and furthermore purely 
imaginary. Given a family of temporal modes Xk which satisfy (13.1311 . one can 
define a new family by a Bogoliubov transformation 

Tk = Kk)Xk + k-{k)xJ^, |A(A:)p - |/i(A:)p = 1 (3.14) 

where the Bogoliubov coefficients A(fc) and /i(fc) have to satisfy the normalisation 
condition above in order for (13.1311 to hold. 

The normalisation condition (13.1311 is chosen in such a way that a pure and 
Gaussian isotropic and homogeneous state for the Klein-Gordon field on a flat 
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FLRW spacetime (satisfying mild regularity conditions) is determined by a two- 
point correlation function of the form [39) 

^2{x,y) = o 3 / ^ W / dk Xkir]a:)Xk{Vy)e"''^'"~'^^e~^'^. 

£ 4,0 8TT-^a[rix)a{r]y) J 

R3 

In particular, we recall that, by Theorem l2.2.2l the causal propagator E{x, y) = 
i~^{uj 2 {x,y) —uj 2 [y,x)) satisfies 

V;^F;(a:, j/)|exe = 5^{x,y) 


for any Cauchy surface S of the spacetime (M, g) with future pointing normal 
vector held iV, where S-s{x,y) is the d-distribution with respect to the measure 
dS on S induced by the covariant volume measure dgX . Observing that {rj} x 
is a Cauchy surface with N = a{g)~^dy and dT, = a{g)^dx, 113.131) ensures that 
the causal propagator satishes 

-^—rdg^E{x,y)\n,=r,y = —^S{x-y). 
airioo) a{rixr 


Choosing a pure, Gaussian, homogeneous and isotropic state uj of the quan¬ 
tized free Klein-Gordon held on a spatially hat FLRW spacetime amounts to 
choosing a solution of (13.121) and (13.131) for each k. In order for w to be a 
Hadamard state, the temporal modes Xk have to satisfy certain conditions in 
the limit of large k which are difficult to formulate precisely. Heuristically, a 
necessary but not sufficient condition is that the dominant part of Xk for large fc, 
when the mass and curvature terms in (13.121) are dominated by is 
i.e. a positive frequency solution. A Bogoliubov transformation (13.141) of a 
Hadamard state dehned by a family of modes Xk is Hadamard if and only if 
fi{k) is rapidly decreasing in k [F7l [63] . 

A particular class of states often discussed in the literature are the adiabatic 
states introduced in [44] . They are specihed by modes of the form 


Xkiri) 


^/20{k,r]) 


exp —z 


f})dfj 


(3.15) 


where f 2 (fc, t]) solves a non-linear differential equation in t] obtained by inserting 
this ansatz into (I3.12|) and hnding 

^{k, vf = fiO{k, g)", H(fc, 77 )', H(fc, g),a{g)) 


for a suitable function /. While this ansatz in principle holds for any state, 
the adiabatic states are specihed by solving the differential equation for Lt{k,g) 
iteratively as 

On+iik, gY = fi^nik, g)'', H„(fc, g)', H„(fc, g),a(g)) 


starting from Llo{k, g) = ^k'^ + + (C ~ g) dla^. Truncating this iteration 

after n steps dehnes the adiabatic states of order n. Note that, while the result¬ 
ing modes satisfy the normalisation condition (13.131) exactly, they satisfy (I3.12p 
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only up to terms which vanish in the limit of constant a or of infinite k and/or 
TO. Thus they constitute only approximate states. This can be cured by using 
the adiabatic modes of order n only for the specification of the initial condi¬ 
tions for exact solutions of (j3.12L see [39]. Regarding the regularity properties 
of such defined ‘proper’ adiabatic states, it has been shown in |2^ (for space- 
times with compact spatial sections) that they are in general not as UV-regular 
as Hadamard states, but that they approach the UV-regularity of Hadamard 
states in a certain microlocal sense in the limit of large n. Consequently, while 
Hadamard states have sufficient regularity in order to compute products and ex¬ 
pectation values of Wick polynomials with arbitrarily many derivatives, one has 
to consider adiabatic states of higher and higher order when dealing with Wick 
polynomials with an increasing number of derivatives, which is conceptually un¬ 
satisfactory. In the following we shall often use the ‘improper’ adiabatic modes 
of order 0, Xo,k{il) = exp(—f Llo{k, fi)dfi))/y/2Llo{k, rf). Adiabatic states have 
also been constructed for Dirac fields, see 111133, and general curved spacetimes 

mus]. 

A further class of states of interest in cosmology, and in fact our candidates 
for generalised vacuum states, are the states of low energy (SLE) introduced in 
[42], motivated by results of [20]. These states are defined by minimising the 
energy density per mode gk 

Qk = (iXfcP + (6f - 1) aH (IxfeP)' + {k^ + - (6^ - 1) H'^a^) |xfcp) 

integrated in (cosmological) time with a sampling function / and thus loosely 
speaking minimise the energy in the time interval where the sampling function 
is supported. The minimisation is performed by choosing arbitrary basis modes 
Xk and then determining the Bogoliubov coefficients A(fc), /i(/c) with respect to 
these modes, such that (for ^ S [0, i]) the resulting modes of the state of low 
energy are 

Xf,k = Mk)xk + d{k)Xk 

with 


X{k) = + 1 = ^|A(fc)P-l 

V v^nky - |c2(fc)|^ ^ 

Ci(fc) = ^J dtf{t)^ (\x'k\^ + (6^ - 1) aH {\xk\^y 
C 2 (fc) = ^ y dtf(t)-^ (x'k^ + (6C - 1) aH (xiy 
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[25] only discusses the case of minimal coupling, i.e. ^ = 0 and proves that the 
corresponding SLE satisfy the Hadamard condition for sampling functions / 
which are smooth and of compact support in time. However, we shall use these 
states for the case of conformal coupling C = and, although we do not prove 
that they satisfy the Hadamard condition, we shall find them to be at least 
regular enough for computing the energy density. Moreover, it is not difficult 
to see that the SLE construction yields the conformal vacuum 

and thus a Hadamard state [25], for all sampling functions / in the massless case. 
This demonstrates both that the SLE construction for ^ = i yields Hadamard 
states at least in special cases and that states of low energy deserve to be 
considered as generalised vacuum states on curved spacetimes. In [52] it is 
conjectured that states of low energy in fact only exist for ^ G [0, g] and that 
they are Hadamard in all these cases. The SLE construction has recently been 
generalised to spacetimes with less symmetry in [58] . 

A conceptual advantage of states of low energy is the fact that they can 
be consistently defined an all FLRW spacetimes at once just by specifying the 
sampling function / once and for all (with respect to e.g. cosmological time 
and a fixed origin of the time axis). Thus, they solve the conceptual problem 
mentioned in Section 12.6.21 namely the necessity to specify a state in a way 
which does not depend on the spacetime in order for the semiclassical Einstein 
equation to be a priori well-defined. Moreover, in m it has been proven that, 
on spacetimes which are asymptotically de Sitter towards t] ^ —oo, every state 
of low energy converges to the Bunch-Davies state (the unique maximally sym¬ 
metric Hadamard state) upon sending the support of /(t(??)) in ij to negative 
infinity. This is a rigorous variant of the statement that every state on de Sitter 
spacetimes converges to the Bunch-Davies state for positive asymptotic times. 

We now proceed to construct the anticipated generalised thermal states on 
the basis of states of low energy. To this avail, we recall a result of [14]: given a 
pure, isotropic and homogeneous state, i.e. a set of modes Xk, one can construct 
generalised thermal states with a two-point correlation function of the form 


u}2{x,y) 


8TT^a{r]x)a{riy 




B3 


[ XkiVx)Xk{Vy) 

I 1 - 


XkiVx)XkiVv) 

g/3feo — 1 


(3.16) 


with _ 

fco = \Jk‘^ + nPafp. 

It has been shown in [14] that for the case of conformal coupling, special FLRW 
spacetimes and particular generalised vacuum modes Xk on these spacetimes, 
these states satisfy certain generalised thermodynamic laws and the Hadamard 
condition, and one can show that they satisfy the Hadamard condition on gen¬ 
eral FLRW spacetimes if the pure state specified by Xk is already a Hadamard 
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state by using results of gTlIM]- Essentially, this follows from the fact that the 
difference of the momentum space integrand in (13.161) and the corresponding 
integrand of the pure state induced by Xk is rapidly decreasing for large k. 

We shall assume in the following that the quantum fields in our model are in 
a generalised thermal state of the form as above, with generalised vacuum modes 
Xk specified by a state of low energy with suitable sampling function /. If m > 0, 
the phenomenological interpretation of these states is that they are the quantum 
state of a massive field which has been in thermal equilibrium in the hot early 
universe and has departed from this equilibrium at the ‘freeze-out time’ a = ap 
and ‘freeze-out temperature’ Tp = ap/p. In the massless conformally coupled 
case, these states are just conformal rescalings of the thermal equilibrium (KMS) 
state with temperature 1//3 in Minkowski spacetime. 

The generalised thermal states we use here have also been constructed and 
analysed for Dirac fields, see [14]. Moreover, we would like to mention that 
several definitions of generalised thermal states on curved spacetimes have been 
proposed so far, including almost equilibrium states |3^ and local thermal equi¬ 
librium states [n ED] 159] . A comparison of these different proposals in the 
context of cosmological applications would certainly be interesting, but is be¬ 
yond the scope of this work. 


3.2.3 Setup and (Computational) Strategy for Approxi¬ 
mately Solving the Semiclassical Einstein Equation 


We now approach the first main result of Section 13.21 a demonstration that 
the energy density in the ACDM-model can be reproduced from first principles 
within quantum field theory in curved spacetime. To this avail, we consider the 
following setup: we model radiation by a conformally coupled massless scalar 
quantum field and matter by a conformally coupled massive scalar quantum 
field. We choose the conformal coupling also for the massive scalar field because 
this considerably simplifies analytical computations and we also found numerical 
computations to be more stable with this value of non-minimal coupling to the 
curvature. Moreover, both quantum fields are assumed to be in generalised 
thermal equilibrium states as introduced in the previous section, where the state 
and field parameters P (possibly different values for the two quantum fields), m 
and Of , as well as the sampling functions / determining the generalised vacuum 
states of the two fields, are considered to be undetermined for the time being. 
Let us stress once more that there is no principal obstruction for formulating 
this model with more realistic quantum fields of higher spin, we just consider 
scalar quantum fields for simplicity and ease of presentation. 

An exact computation of the energy density of the two quantum fields in 
the generalised thermal states would require to solve the coupled system - the 
so-called backreaction problem - consisting of the quantum fields propagating on 
a FLRW spacetime, which in turn is a solution of the semiclassical Friedmann 
equation 


= 


8ttG 


+ g^) , 


3 


(3.17) 
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where = a;™(: = uj°{:Tqq:) are the energy densities of the two quan¬ 

tum helds in the respective generalised thermal states and the 00-component 
of the stress-energy tensor is here taken with respect to cosmological time t. 
An exact solution of the backreaction problem is quite involved, as it requires 
solving simultaneously the mode equation (13.1211 for all k and the semiclassical 
Friedmann equation. Notwithstanding, there have been quantitative numerical 
treatments of the backreaction problem, see e.g. [1131 [3 [5], as well as nu¬ 
merous qualitative treatments including |19j . where the backreaction problem 
in FLRW spacetimes is set up in full generality from the point of view of the 
algebraic approach to QFT on curved spacetimes, [13], where the same point 
of view is considered and the coupled system is solved exactly for conformally 
coupled massless scalar quantum fields and approximately for massive ones, and 
|46l|49], where exact solutions of the backreaction problem are shown to exist. 

However, in this work we follow a simplified strategy in order to avoid solving 
the full backreaction problem, which is justified in view of our aim. We assume 
that the two quantum fields in our model are propagating on a FLRW spacetime 
which is an exact solution of the Friedmann equation in the ACDM-model, i.e. 


1/2 

7/2 


PACDM 

QO 


= Oa 






-,4 ’ 


(3.18) 


where Hq ~ 10“^^eV denotes the Hubble rate of today, the so-called Hubble 
constant, go — 10“^^eV^ is the energy density of today and Ha, Hrad and 
Hmat denote respectively the present-day fraction of the total energy density 
contributed by the cosmological constant, matter and radiation, cf. Section [3.1l 
For definiteness we consider the sample values Hmat = 0.30, Hrad = 10“"^, Ha = 
1 — Hmat ~ Hrad, rather than currently measured values from e.g. the Planck 
collaboration |5], because the exact values are not essential for our results. Given 
this background spacetime, we aim to prove that the field and state parameters 
of our model, as well as the SLE sampling functions, can be adjusted in such a 
way that the energy density of the quantum fields in our model matches the one 
in the ACDM-model up to negligible corrections for all redshifts z € [0,10®], 
i.e. 




, Hy,2 , H,. 

— "I r ”1 T = 


PACDM 


go a-^ o'* go 

Once we succeed to obtain this result, we have clearly solved the full coupled sys¬ 
tem in an approximative sense to a good accuracy in particular. We recall that 
the restriction to post-Big Bang Nucleosynthesis (BBN) redshifts z € [0,10®] 
is motivated by the fact that interactions between matter fields are assumed to 
have negligible large-scale effects in this period of the cosmological history, cf. 
Section 11 .2.11 

In order to compute the quantum energy density p®-l-p™, we start from (13.101) 
which parametrises the freedom in defining the energy density as an observable 
and gives a possible ‘model definition’. The renormalisation freedom for the 
energy density is readily computed via goo = — 1, Gqo = and 


Joo = \l.. = m^-l2HH-mHHV 


( 3 . 19 ) 
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In order to compute the energy density for each quantum field following from 
(|3.10[) . one has to first subtract the geometric Hadamard parametrix from the 
twO”point function of the given state and then to apply a suitable bidifferential 
operator followed by taking the coinciding point limit. As the states we consider 
here are given as integrals over spatial momenta, it seems advisable to try to re¬ 
write the geometric Hadamard parametrix also in this form, in order to perform 
a mode-by-mode subtraction and a momentum space integral afterwards. This 
is indeed possible, as elaborated in [m mills [501 US |63]. The details are quite 
involved, thus we omit them and present directly the result. To this avail we 
follow m, where results of [50] are used. In [15] only the minimally coupled 
case ^ = 0 is discussed, but it is not difficult to generalise the results there to 
arbitrary 

Doing this, we find the following result for the total energy density of the 
massless and massive conformally coupled scalar fields in the generalised thermal 
states. 


g° + g^ 

go 


- gvac 


- gvac 




^^gth 






go 




Jqo 
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(3.20) 


^ “ 3607r2 


Ha 


SirGai 


5 = 


87rG'a2 


e 


St:GHI 

3 


(3Q!3 + 0:4) . 


Here Ha, 6 and e parametrise the freedom in the definition of the energy den¬ 
sity as per (13.101) . The number of free parameters in this equation has been 
reduced to three, because I^v and are proportional in Robertson-Walker 
spacetimes, cf. (I3.19p . As already discussed in Section 13.2.11 we omit the free¬ 
dom parametrised by 5 in the following, as it renormalises the Newton constant 
and we consider this to be already given as an external input. For now we will 
also neglect the contribution parametrised by e, as it turns out to be negligible 
for 0 < e <C 1 ; we will analyse the influence of this new term, which does not 
appear in the ACDM-model, separately in the next section. Thus, for the re¬ 
mainder of this section, Ha parametrises the residual freedom in the definition 
of the quantum energy density. 

The term proportional to 7 , which is also not present in the ACDM-model, 
is the contribution of the trace anomaly, cf. Section 12.6.31 This term is fixed 
by the field content, i.e. by the number and spins of the fields in the model and 
always proportional to barring contributions proportional to Jqo which we 
prefer to subsume in the parameter e. We have given here the value of 7 for 
two scalar fields, see Table 1 on page 179 of nni for the values in case of higher 
spin. As 7 ~ and H < Hqz^ in the ACDM-model for large redshifts, 

this term can be safely neglected for z < 10 ®; we will comment on the relevance 
of this term for z ^ 10® in Section [3.2.51 

Finally, the remaining terms in (13.2011 denote the genuinely quantum state- 
dependent contributions to the energy densities of the two quantum fields. We 
have split these contributions into parts which are already present for infinite 
inverse temperature parameter /3 in the generalised thermal states, and thus 
could be considered as contributions due to the generalised vacuum states (£>^aci 
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£'gvac)j into the remaining terms, which could be interpreted as purely ther¬ 
mal contributions fi'gth)- Note that fi’gvac not uniquely defined 

in this way, but only up to the general renormalisation freedom of the quantum 
energy density, i.e. one could “shuffle parts of” LIa, S and e into e.g. and 
vice versa, without changing any physical interpretation of the total energy den¬ 
sity. We have have fixed the renormalisation freedom of fi^ac (and thus £»gvac for 
TO = 0) in such a way that it has a particularly simple form in the conformally 
coupled case, cf. (I3.23|) . With this in mind, the state-dependent contributions 
read as follows, where the massless case is simply obtained by inserting to = 0, 
and we give here the result for arbitrary coupling ^ for completeness. 



0 


(3.21) 


-f (fc^ -f TO^a^ - (6^ - 1) \xk\'^) 

k rrV' — — 1) 

2a^ Aa^k 



+ 36ij2 - 72HH 
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+ (A:2 + TO^a^ - (6^ - 1) ifV) Ixfcp) 


In the conformally coupled case C = | simplifies considerably to 


OO 



(3.23) 
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where xo,k are the adiabatic inodes of order 0, cf. the previous section. This 
implies that the so-called Hadamard point-splitting regularisation of the energy 
density coincides with the so-called adiabatic regularisation of order zero up to 
the trace anomaly term and terms which can be subsumed under the regulari¬ 
sation freedom. 

3.2.4 Computation of the Energy Density 

In the following we shall analyse the individual state-dependent terms in the 
energy density. 

Computation of 

Following our general strategy, we first aim to show that in states of low en¬ 
ergy defined by a sampling function of sufficiently large support in time on the 
FLRW spacetime specified by (13.181) . is for all z S [0,10®] negligible in 
comparison to the total energy density in the ACDM-model. Results in this 
direction have been reported in |15j for the simplified situation of a de Sitter 
spacetime background (corresponding to flmat = f^rad = 0). Here we gener¬ 
alise these results to ACDM-backgrounds. One can easily see that = 0 in 
the case of m = 0. For masses in the range of the Hubble constant m ~ Hq 
and states of low energy we have performed numerical computations and found 
fi'^ac/fi'ACDM < 10“^^®, see Figures Ism 13.21 To achieve this result, we have 
rewritten all expressions in terms of the redshift z as a time variable and solved 
the equation (13.131) with initial conditions at z = 0 given by the value and deriva¬ 
tive of the adiabatic modes of order zero Xk,o there. Note that a state of low 
energy does not depend on the choice of a mode basis, but the choice we made 
seemed to be numerically favoured. To fix the state of low energy, we chose a 
sampling function which was a symmetric bump function in z supported in the 
interval z € (10“®, 10“® +10“"^) for definiteness. In order to make the numerical 
computations feasible, we chose a logarithmic sampling of k with 10® sampling 
points, where the boundaries of the sampling region have been chosen such that 
the integrand of cf. (13.231) . was vanishing in fc-space to a large numerical 
accuracy outside of the sampling region for all z G [0,10®]. We have computed 
the mode coefficients Ci(fc) in the mode basis chosen at each sampling point by 
a numerical integration in z and finally the energy density by means of a sum 
over the sampling points in fc-space. Thus we have approximated the integral 
in (13.231) by a Riemann sum with logarithmic sampling. As our main aim here 
is to demonstrate that g^ac/fi’ACDM 1 in general, we have not performed an 
extensive analysis of the dependence of g^^^. on the width of the sampling func¬ 
tion, but we have observed that the maximum amplitude of p^ac/P acdm seems 
to be monotonically growing with shrinking width of the sampling function, in 
accordance with the computations of |15j in de Sitter spacetime. 

Unfortunately, we have not been been able to compute fii^ac/ffACDM for 
m > 10®iJo in the way outlined above because for large masses the modes os¬ 
cillate heavily, and thus it costs a lot of computer power to solve the mode 
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z 


Figure 3.1: A£i™ 3 ^j,/eACDM for z < 1 for various values of m (rescaled for ease 
of presentation). The dotted line corresponds to m = IOOHq and A = 10“^, 
the dashed line to m = 10 Hq and A = 1 and the solid line to m = Hq and 
A = 10^. One sees nicely how the energy density is minimal in the support of 
the sampling function at around z = 10“^. 


equation for such a large z-interval we are interested in and to the numerical 
accuracy which is necessary to obtain reliable results for the coefficients of the 
state of low energy and g^ac/Q acdm, see also [52l Section 8.4] for related con¬ 
siderations. However, realistic field masses in the GeV regime are rather of the 
order of 10'^^ Hq. In the numerical computations outlined above we have ob¬ 
served that g^ac/^lACDM seemed to grow quadratically with m, see Figure 15^ 
but looking at the results of m in de Sitter spacetime one could maybe ex¬ 
pect that g^ac/^lACDM decreases for large masses. Moreover, even if a potential 
quadratic growth of fi'^ac/fi'ACDM with m would still imply ^?^ac/^?ACDM "C 1 
for realistic masses and given ^?^ac/^?ACDM 10“^^° for m = Hq, it would be 
better to have a more firm understanding of the large mass regime. 

In view of the numerical problems for large masses we had to resort to an 
approximation in order to be able to compute g^ac/2 acbm- In fact, we have 
taken the adiabatic modes of order zero as basis modes for computing the state 
of low energy. Of course these modes are not exact solutions of the mode 
equations, but the failure of these modes to satisfy the exact mode equation is 
decreasing with increasing mass and thus one can expect that the error in all 
quantities derived from these modes rather than exact modes is also decreasing 
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Figure 3.2: p^ac/PACDM for z > 1 for various values of m. The upper line 
corresponds to to = IOOHq, the middle lines to to = IOHq and the lower lines to 
TO = iJo; solid lines (dashed lines) indicate results obtained with exact modes 
(zeroth order adiabatic modes). p^ac/fi'ACDM becomes constant for large z 
because there both energy densities scale like a~‘^, c.f. (13.251) and the related 
discussion. 


with increasing mass. We have checked numerically that the energy density 
computed with adiabatic modes rather than exact modes matched the ‘exact’ 
result quite well already for masses in the regime to ~ Hq, see Figure [XU For 
more details regarding error estimates for adiabatic modes we refer the reader 
to [43]. 

Inserting the adiabatic modes xo,k we obtain the following expressions for 
the coefficients Ci{k) of the states of low energy. 


ci(fc) = J dzf{z) 


mrH 


16r2o(fc)®(l + z) 


+ 


C2(fc) = J dzf{z) 


nrH 


m 


— i- 


16r2o(^)^(l H"-2:) 4r2o(^) 


fIo(fc)(l + zf 
2H 


exp ( -2i r ^^dz' 


We now perform a further approximation. We take as a sampling function 
a Gaussian with mean zq and variance tr <C 1 


/W = 


1 




exp - 


[z - Zof 

2 ct 2 
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and take the zeroth term of the Taylor expansion of both the expressions in 
the curly brackets in the integrands of Ci{k) and of the integrand appearing in 
the exponent of the exponential in C 2 {k). Without performing a detailed error 
analysis we note that this is justified for tr <C 1 because the higher coefficients of 
the associated Taylor series differ from the lowest coefficient roughly by factors 
of either dzH/H\ 2 =zo or H{zo)/m, both of which are either smaller than or of 
order one under the assumption of large masses and a ACDM-background. We 
can now compute the z-integrals, which corresponds to considering the Fourier- 
transform of / in the case of C 2 . Using Ho/m <C 1 (and thus H{zo)/m <C 1), 
we can estimate the resulting coefficients as follows 


Cl (A:) > 1 


|c 2 (fc)| < exp - 


HizofJ^'^^y H{zo?{l + zo)\ 

For H{zo){I + zo)/{ma) <C 1, |c 2 (fc)| <C 1 and we can approximate the Bogoli- 
ubov coefhcients \{k) and /r(fc) as /r(fc) ~ ^2ci^k) ’ ^ estimate 

eS-ac as 


le^acl J dkk^iu^ + ^^\M) i\Xo,k\^ + ^o\Xo,k\^) J dkk'^fj.\\\no 
0 0 


1 H{zoVm 
<-I --exp 


H{zo)Hl + zoy 


such that, barring our approximations, we indeed get a result which shows that 
the energy density decreases - exponentially - for large masses. Note that for 
not too small a the bound we found is in general small compared to pacdm 
even if we forget about the exponential because Horn is much smaller than the 
the square of the Planck mass, i.e. 1/G. We also see that the bound grows with 
growing zq, i.e. if we ‘prepare’ the state of low energy further in the past, and 
that it diverges if the width of the sampling function goes to zero; this is in 
accord with the results of m in de Sitter spacetime. Note that we could have 
chosen any rapidly decreasing or even compactly supported sampling function 
in order to obtain a bound which is rapidly decreasing in m/H q, thus one could 
say that the result does not depend on the shape of the sampling function as 
long as its width is not too small. Finally, one could of course directly take the 
point of view that for large masses the adiabatic modes xo,k define ‘good states’ 
themselves and conclude that in these states — 0 on account of (I3.23|) . 


Computation of P^herm 

We now proceed to analyse the thermal parts of the state-dependent contribu¬ 
tions to the total energy density. Inserting ^ = i in (I3.22p . we find 

OO 

= Jdkk‘^ ^^k!_ 1 (IXfcP + (fc^ +mV) Ixfcp) (3.24) 

0 












100 


CHAPTER 3. COSMOLOGICAL APPLICATIONS 


with fco = . Before performing actual computations, we would 

like to mention a general result about the scaling behaviour of the energy 
density with respect to a [48]. To wit, using the equation of motion (I3.13|) 
and the assumption that H > 0, one can compute the derivative of Qk = 
IXfeP+(fc^ + IXfcP with respect to a and obtain the following inequalities 


Qkja = 1) ^ Qkja) ^ Qkja = 1) 

k2 _|_ ^2 q 4 “ 


(3.25) 


From these one can already deduce that has a scaling behaviour with respect 
to a which lies between and and approaches in the limit of vanishing 
a, in fact this still holds if we replace the Bose-Einstein factors in the generalised 
thermal states by arbitrary functions of k. Moreover (I3.25|) also implies that 
can not scale with a power of a lower than —4 for small a on ACDM 
backgrounds, cf. (13.231) . 

Proceeding with actual computations we find that in the massless case 
can be computed exactly and analytically and the result is 

As in the massless case the state of low energy is the conformal vacuum and the 
associated generalised thermal state is the conformal temperature state with 
temperature parameter /3 = 1/T, this result in fact holds for fields of all spin, 
i.e. the generalised thermal energy density in this case is always the one in 
Minkowski spacetime rescaled by a~^. Thus a computation with e.g. photons 
or massless neutrinos yields the same result (13.261) up to numerical factors due 
to the number of degrees of freedom and the difference between Bosons and 
Fermions. 

In the massive case it is not possible to compute analytically and exactly, 
but we have to resort to approximations once more. We recall that the massive 
scalar field in our model should represent baryonic matter and Dark Matter in a 
simplified way. Thus we take typical values of /3, ap and m from Chapter 5.2 in 
[34] computed by means of effective Boltzmann equations. A popular candidate 
for Dark Matter is a weakly interacting massive particle (WIMP), e.g. a heavy 
neutrino, for which |34j computes 


xp =/lapm 15 + 3 log(m/GeV), (3.27) 

ap ~ 10“^^(m/GeV)“^. 

We shall take these numbers as sample values although working with a scalar 
field, because for large masses m » i/o, the ‘thermal energy densities’ g^i^ 
in generalised thermal states for free fields of spin 0 and ^ can be shown to 
approximately coincide up to constant numerical factors on the basis of the 
results of [2] and [131 Section IV. 5]. 

Considering m > IGeV, we can compute g^i^ approximatively as follows. 
We recall from the computation of ^^at for large masses Hq one can 
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consider the adiabatic modes of order zero xo.fc as approximative basis modes 
for the computation of the state of low energy and that with respect to this basis 
one finds for the coefficients of the state of low energy A ~ 1, ^ ~ 0, thus we 
can insert those modes in (j3.26p instead of the modes of the state of low energy. 
Using Hq once more, we have |xo feP + + m?a^) Ixo.feP — + m?a^ 

and using xp >15 we can approximate the Bose-Einstein factor in (I3.26P as 
l/(g^fco _ 1 ) ~ Finally, we can rewrite the integral in (13.261) in terms 

of the variable y = k/^apm) and compute, using a/ap >> 1 for the redshift 
interval z G [0,10®] we are interested in, 


apiri 


“ 27r2 flS 


OO 

0 


This already gives the desired result oc a~^. The remaining integral can 
be computed numerically, however, for xp >> 1 only y <C 1 contribute to the 
integral and one can approximate -^/y^ + 1 ~ 1 +y'^ / 2 and compute 


dgth - 


(27r)3/2 /33 q3 


which for a = ap = 1 (unsurprisingly) coincides with the thermal energy density 
for massive scalar fields in Minkowski spacetime. 


The Total Energy Density 

Collecting the results of this section, we find for the total energy density of our 
model 


gO_^gm^ 1 m I 1 TT^ 1 

go ^ (27r)^/2 pfa^go ^ ’HDP'^go oA ’ 

where we wrote /3i, (32 in order to emphasise that the generalised thermal states 
for the massive and massless conformally coupled scalar fields can have differ¬ 
ent temperature parameters /3. We recall that the thermal contribution of the 
massless scalar field has been computed exactly, while the one of the massive 
scalar field is an approximative result. The above result shows that we indeed 
succeeded in modelling radiation by a massless scalar field and matter by a mas¬ 
sive scalar field in suitable generalised thermal states. Obviously, we can choose 
the free parameters m, At, xp in such a way that the prefactors of the matter 
and radiation terms have their correct ACDM-values Umat and Urad, e.g. for 
the former we could choose the sample values (13.271) with m ~ lOOGeV, and for 
the latter 1//3 ~ IK, i.e. the temperature of the CMB. Finally, we model Dark 
Energy simply by a cosmological constant, which in our context appears as a 
parametrisation of the freedom in defining energy density as an observable. 

Since our description of the the standard cosmological model within quantum 
field theory on curved spacetime reproduces the energy density of the original 
ACDM-model up to negligible corrections, it can obviously be matched to the 
observational data as good as this model. 


















102 


CHAPTER 3. COSMOLOGICAL APPLICATIONS 


3.2.5 Deviations from the Standard Model and their Phe¬ 
nomenological Conseqnences: eJoo and Dark Radia¬ 
tion 


Our analysis in the previous sections implies that there exist quantum states 
in which the total energy density in quantum field theory on curved spacetimes 
differs from the one in the ACDM-model only by the higher derivative term e Jqo 
and terms which are generally negligible or become important only at redshifts 
z 10®. The prefactor e of Jqo is not determined by the theory but a free 
parameter so far and it seems advisable to study its impact on the cosmological 
expansion. Indeed, as the second main result of Section 13.21 we demonstrate 
in the following that such a term can provide a natural explanation of Dark 
Radiation. The fact that the contribution eJoo to the energy density can look 
like radiation for large z has already been observed in e.g. [33], but to our 
knowledge the relevance of this for the phenomenon of Dark Radiation has not 
been discussed so far. 

At this point we would like to mention that in some other works, see e.g. m 
and references therein, the parameter e is not considered to be a free parameter, 
but it is rather taken as determined by the field content of the QFT model, just 
like the parameter 7 in (13.201) . This is motivated by the fact that most common 
computational schemes for regularising the quantum stress-energy tensor yield 
the same result for e, which is thus taken to be the correct value. In this work 
we follow the point of view of [601ES] and thus start from the premise that there 
is no physical principle within QFT on curved spacetimes which determines the 
value of e a priori; thus in particular we do not attribute any physical ‘meaning’ 
to a specific computational scheme, even if most common schemes give the same 
result. 

To start our analysis, we briefly review the notion of Dark Radiation and 
the related observations. The fraction Dj-ad of the radiation energy density in 
the ACDM-model is computed as 


^rad — 




(3.28) 


where r2.y ~ 5 x 10“® is the fraction due to electromagnetic radiation, which can 
be computed by inserting into (13.261) the CMB temperature Tcmb — 2.725K, 
dividing by today’s energy density go = SHq/^SttG) ~ 1.33 x 10“^^eV (and 
multiplying by two for the two degrees of freedom of the photon). Moreover, Ngg 
is the number of neutrino families and the factor 7/8(4/11)'^/^ = 0.2271 takes 
into account that neutrinos are Fermions and ‘colder’ than the CMB photons, 
because they have decoupled from the hot early bath in the universe earlier than 
electrons and positrons and have thus not been ‘heated up’ by the decoupling 
of the latter in contract to the photons. The standard value for Nqs is not 3 
as one would expect, but rather TVeff = 3.046 because the value 7/8(4/11)^^/^ 
in (13.281) is computed assuming e.g. instantaneous decoupling of the neutrinos 
and corrections have to be taken into account in a more detailed analysis |40] : 
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it is customary to take these corrections into account by considering N = 3.046 
as the standard value of the ‘neutrino family number’ rather than changing the 
factor 7/8(4/11)'^/^ in this formula, hence the nomenclature Neff. Consequently, 
it is convenient to parametrise any contribution to flrad which is not due to 
electromagnetic radiation and the three neutrino families in the standard model 
of particle physics by AiVegf = Neff — 3.046. 

One of the two main observational inputs to determine flrad and thus Ngff is 
the primordial fraction of light elements in the early universe as resulting from 
the so-called Big Bang Nucleosynthesis (BBN), which has occurred at around 
z ~ 10® and thus in the radiation-dominated era. This is the case because 
the nucleosynthesis processes which happened at that time depend sensitively 
on the expansion rate H ~ Ho^Llrad/o,^, see e.g. [331 IMl [IZ]- The other 
main observational source for the determination of Neff is the cosmic microwave 
background radiation (CMB). This radiation was emitted at about z ~ 1100, 
but the CMB power spectrum is sensitive to the expansion before this point, e.g. 
to the redshift Zeq at which the energy densities of matter and radiation were 
equal, see [2l Section 6.3] and the references therein for details; for standard 
values, Zeq — 3000. 

The observations to date do not give a conclusive value for Neff, and the 
value inferred from observational data depends on the data sets chosen. The 
Planck collaboration [2] reports e.g. values of Neff = 3.36^0^4 at 95% confidence 
level from combined CMB power spectrum data sets, Neff = 3.52+Q 45 at 95% 
confidence level from combining these data sets with direct measurements of 
the Hubble constant Hq and of the power spectrum of the three-dimensional 
distribution of galaxies (so-called baryon acoustic oscillation, BAO), and Neff = 
3.41 ± 0.30 at 68% conhdence level from combining the CMB power spectrum 
data sets with BBN data. Yet, one can infer from these values that there is a 
mild, but not very significant, preference for ANeff > 0. Thus there has been an 
increasing interest in models which can explain a potential excess in radiation 
and thus ANgff > 0. Most of these models assume additional particles/fields, 
e.g. a fourth, sterile, neutrino, whereas other consider geometric effects from 
e.g. modifications of General Relativity. Moreover, in most models ANgff is 
constant and thus affects BBN and CMB physics alike, while in others ANgff is 
generated only after BBN and thus affects only CMB physics. 

In the following we shall propose a new and alternative explanation for Dark 
Radiation which follows naturally from our analysis of the ACDM-model in 
quantum field theory on curved spacetimes and has the interesting characteristic 
that it generates a value of ANgff which increases with z and thus affects BBN 
physics more than CMB physics. 

Following the motivation outlined at the beginning of this section, we solve 
the equation 


WN ~ 


0 


rad 


yoo 

-^0 


(3.29) 


which can be rewritten as a second order ordinary differential equation for H 
in z, numerically with ACDM-initial conditions H{z = 0) = Hq, dzH{z = 0) = 
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Figure 3.3: ANeff{z) depending on e for z = 10® (BBN, solid line) and z = 3x 10® 
(CMB, dashed line). For e < 0 and e positive and large enongh, the values at 
the two redshifts coincide because the maximum value of A7Veff(z) is reached 
already for z < 3 x 10® in these cases. 


77o(3flniat + 4flrad)/2. As before, we consider for definiteness flmat = 0.30, 
r^A = 1 — ~ because the exact values of these parameters are not 

essential for our analysis. Looking at the characteristics of the solution to this 
ordinary differential equation, it turns out that a non-zero e generates a time- 
varying ANeS > 0. In more detail, we define for the solution H of (j3.29|) 


AAeff(z) 


~ f2niat(l + 2:)® — r2rad(l + 

-^0 _ 

0.2271(1 + z)4 


and sample this observable at the redshift z = 10® associated to BBN physics 
and at the redshift z = 3000 associated to CMB physics. We collect our results 
in Figure 1331 

As can be inferred from this hgure, A7Veff(z) is monotonically increasing in 
e, where positive and negative values of e result in very different behaviours. For 
positive values of e one hnds that AiVeff(z) vanishes in the limit of vanishing 
e, as one would expect. On the other hand, it turns out that for negative 
values of e, ANes{z) diverges as e approaches zero. While this seems to be 
puzzling at first sight, it fits well with previous qualitative analyses of the effect 
of the higher derivative term Jqo- In fact, it is known that the inclusion of 
this higher derivative term can lead to unstable solutions of the semiclassical 
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Figure 3.4: A7Veff(^) for e = 2 x 10" 


Einstein equations, where for e < 0 (e > 0) the class of solutions we consider 
here, effectively fixed by the ACDM initial conditions, turns out to be unstable 
(stable), see e.g. [U [221 |2S1133 1331133] and the discussion at the end of this 
section. Thus, the divergence of ANesiz) as e approaches zero from below can 
be just interpreted as a sign of this instability. 

In [STl [T3] , e = 0 has been chosen on conceptual grounds in order to discard 
unstable solutions altogether. However, as we see here, a non-zero e can have 
interesting phenomenological implications. After all, taking quantum field the¬ 
ory on curved spacetimes seriously, e is a free parameter of the theory, which 
we can only fix in a more fundamental theory or by observations. Indeed, 
we see in Figure 13.31 that e < 0, corresponding to an unstable solution of the 
semiclassical Einstein equation, is already ruled out by observations because it 
generally leads to ANes{z) ^ 1 which is certainly not compatible with value of 
ANes — 0.5 — 1.0 inferred from observations as mentioned above. On the other 
hand we see that in order to not exceed AA"eg = 1 at both BBN and CMB we 
have to choose 0 < e < 2 x 10 “^®, thus, without performing a detailed fit of 
BBN and CMB data, we can say that the values for ANes reported e.g. by the 
Planck collaboration in [2] give an upper bound of about 2 x 10“^^ for e. We 
plot ANes{z) for for this value of e and redshifts 0 < z < 10® in Figure l3Hl As 
already anticipated in Figure [3731 one can nicely see how ANeg(z) is monoton- 
ically growing in z, with ANgg^z = 0 ) = 0 as fixed by our initial conditions. 
Moreover one can see clearly that if one wants to meet the bounds on ANgS at 
the BBN redshift, the excess in the effective number of neutrinos at the CMB is 
negligible, which is the characteristic signature of this potential explanation for 
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Dark Radiation. We have not considered the influence of the initial conditions 
for (13.291) on AA^eff, but we expect that for the initial conditions compatible with 
low-z observational data such as Supernova type la data and Baryon Acoustic 
Oscillation data, ANes will not differ considerably from the form we found as 
these data will not allow for large deviations in the initial conditions from the 
ACDM ones we chose. 

As a further, rather pedagogical remark, we would like to comment on the 
fact that, for large absolute values of e, ANgg does not depend on the sign 
of e, as can be seen from Figure. 13.31 This phenomenon can be understood 
as follows. Naturally, for large absolute values of e, the other terms in (I3.29|) 
become negligible and one effectively solves for Jqq = 0. The solution of this 
ordinary differential equation with initial conditions H{z = 0) = c, dzH{z = 
0 ) = d is 


H{z) 


22/3 


and thus, inserting the ACDM initial conditions H{z = 0) = Hq, dzH{z = 0) = 
7Jo(3Dmat + 4Drad)/2 we, find for large z, Nee{z) ~ 10“^ as in Figure [?31 

The above analysis may be interpreted as providing an upper bound for 
the parameter e in eJoo, which quantifies higher order derivative corrections to 
classical General Relativity, based on BBN and CMB data. Two further ways 
to obtain bounds on e are known to us. To this avail, let us recall that the 
tensors and in (13.101) can be obtained as variational derivatives with 
respect to the metric of the Lagrangean densities y/^R? and 
Consequently, an upper bound on e can be obtained by studying modifications 
of the Newtonian potential arising from such higher-derivative modifications of 
the gravitational action. To wit, the Lagrangean 


L = 


i— 

ViGttG 


+ C\R? + C2RauR^‘' 


leads to the Newtonian potential of a point mass m 


—mG 


4 

~ 3® 


(3.30) 


TOi = 


TO 2 = 


i/167rGc2 


v'327rG(-3ci - C2) 

Using recent data m from torsion-balance experiments to test the gravitational 
inverse-square law at ^ 10“^m and assuming that the two Yukawa corrections 
don’t cancel each other at this length scale, one obtains —ci,C 2 < 10®^ [T^ . 
To compare this with our results, we recall that in our treatment these higher 
curvature terms appear on the right hand side of the semiclassical Einstein 
equation and that we have computed in units of Hq, thus we have 


(—3ci — C2)8ttGHq 


~ (— 3 ci — C2) X 10 
















3.2. THE COSMOLOGICAL EXPANSION IN QFT ON CURVED SPACETIMES107 


which would imply e < 10“®° and thus a stronger bound than the one we 
inferred from cosmological observations. Of course such a low value of e leads 
to AA^eff 1 at both BBN and CMB and thus the amount Dark Radiation 
generated by this term would be negligible. 

Notwithstanding, there are still several aspects of our analysis which are 
of interest. First of all, our bound on e is completely independent from the 
one inferred from laboratory experiments and can thus be considered as an 
additional confirmation of those results. Moreover, it is still possible that the 
Yukawa corrections in (I3.30|) cancel each other on the length scales relevant for 
the experiments described in m, such that e could be as large as our upper 
bound, which in this case would give a real bound on one and hence both Yukawa 
corrections. Finally, the bounds inferred from |31] and from our analysis stem 
from phenomena on completely different length scales. As a rough estimate 
we note that the diameter of our observable universe, which today is about 
6/Ho ~ lO^^m, was at e.g. z = 10° still 10^® m and thus much larger than the 
submillimeter scales relevant for the torsion-balance experiments. Thus it could 
be that effects we have not considered so far, e.g. state-dependent effects which 
are due to the small-scale structure of the quantum states we have fixed only 
on cosmological scales so far, affect the comparison between the two different 
sources of input for the determination of e. 

Finally, a further bound on e can be obtained by analysing the effects of 
higher derivative contributions to the gravitational Lagrangean in the context 
of Inflation. In fact, the first inflationary model proposed by Starobinsky in 
[54] was based on the observation that the semiclassical Friedmann equation 
for massless (and in the scalar case conformally coupled) quantum fields in the 
conformal vacuum state, namely. 


H^ 


= 7 




_Joo 

^0 


Oa : 


(3.31) 


has solutions which correspond to a de Sitter phase a{t) oc exp(iJ_|_t) with 
~ Ho/^// which for e > 0 ‘decays’ to a de Sitter phase a{t) oc exp(iJ_t) 
with H- ~ Ho^/LIa. For various reasons, this model was discarded as a viable 
inflationary model, see e.g. the discussion in Section 9.4 of [35], and Starobinsky 
instead proposed an explanation for Inflation based only on (I3.3ip with 7 = 0 
in [55] : thus Inflation in this model occurs due to an i?°-modification of the 
gravitational action. This model is now widely known as the Starobinsky model 
and still in very good agreement with CMB data |3]; interestingly, it is mostly 
attributed to Ref. |S3j rather than Ref. [5S] , although the two models are quite 
different. 

Neglecting the trace anomaly term jH'^/Hq in (I3.3ip can be justified if 
e 7 and if dzH/H is not too small, but note that 7 > 0 is in contrast to 
e > 0 an o priori prediction of QFT on curved spacetimes and that the trace 
anomaly term will always dominate for z large enough in almost all solutions 
of (]3.31l) . However, this occurs generally where H of the order of the Planck 
mass Mp = l/yfC ~ I0^°GeV (recall that 7 ~ in the case of a scalar 
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field and thus H+ ~ ~ Mp in the original Starobinsky model), and 

one can argue that the semiclassical Einstein equation should not be trusted in 
this regime, see e.g. m- Notwithstanding, taking (13.311) seriously and solving 
it with ACDM-initial conditions, i.e. requiring that today we are in the ‘lower 
de Sitter branch’ ruled by H-, already leads to the lower bound e > 0 because 
with e = 0 and these initial conditions, H becomes imaginary for large z, cf. 

ISI1II2- 

A more concrete lower bound on e can be obtained by confronting the 
Starobinsky model with CMB data. In this model, the parameter here called e 
(not to be confused with the usual slow-roll parameter e in Inflation) is quan¬ 
tified in terms of the ‘scalaron mass’ [54l [55] 


M = 


Ho 


whose magnitude determines the amplitude of the primordial fluctuations, which 
turns out to be linear in M [35] . Comparison with CMB data yields [50] 

M ~ 10"^Mp ^ e ~ 10"“^. (3.32) 


Thus, if Inflation occurred due to the e Jqo contribution to the energy density, 
then e is too small for generating a considerable amount of Dark Radiation. 
However, if Inflation has a different origin, then we can interpret the value in 
(|3.32|) as an an upper bound on M and thus obtain the lower bound e > 

After discussing potential phenomenological implications of the energy den¬ 
sity contribution eJoo which is not present in the ACDM-model, we note that 
potential further corrections to the ACDM-model can come from specihcs of 
the quantum state we have neglected in our analysis. We have chosen the quan¬ 
tum states in our discussion such that their characteristic energy density was 
entirely of thermal nature, but we have seen that also pure, non-thermal states 
can have contributions to the energy density which scale like a“^, cf. Figure 
o It could be that there exist states which are compatible with observations 
and have sizable energy-density contributions of this kind; these states would 
then provide a further alternative explanation for Dark Radiation which does 
not call for the introduction of new fields and the associated particles. 


3.3 A Birds—Eye View of Perturbations in Infla¬ 
tion 

A prominent application of quantum field theory in curved spacetime is the 
analysis of perturbations in Inflation. In the simplest models of Inflation, the 
characteristic exponential expansion is driven by the energy density of a classical 
scalar field ip. The perturbation (j) of this held, combined in a gauge-invariant 
way with the perturbations 7 of the classical metric g, are considered as quantum 
helds propagating on the classical background M^p = {M,g,p), where M = 
{M,g) is a spatially hat FLRW spacetime and p depends only on time in the 
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homogeneous and isotropic FLRW coordinates. Rephrased in more abstract 
terms, one may say that the analysis of perturbations in inflation consists in 
quantizing the field theory of the tuple 0 = {g, (p)* perturbatively around a 
background which satisfies the classical Einstein equation and is of FLRW-type. 
This perturbative quantum field theory is usually truncated at linear order. 

The standard textbook treatment of the quantum theory of perturbations in 
Inflation, see e.g. m, consists in using the FLRW-symmetry of the background 
in order to split the metric perturbation 7 into components which transform 
as scalars, vectors and tensors under the isometry group of FLRW-backgrounds, 
the Euclidean group E^. Subsequently, gauge-invariant linear combinations of 
these components and the scalar field perturbation cj) are identified and quan¬ 
tized in a canonical fashion. As the above-mentioned splitting is non-local and 
depends heavily on the FLRW-symmetry, it is a priori not clear whether it cap¬ 
tures all local observables of the theory. A systematic analysis of this issue from 
the point of view of algebraic quantum field theory has been performed in |24j . 
In this section we shall review the main steps and results of this analysis. For 
simplicity we consider the special and commonly assumed case where the scalar 
field is minimally coupled to the metric, the case of general coupling is treated 
in |24j . 


3.3.1 The Linearised Einstein—Klein—Gordon System on 
General Curved Spacetimes 

We first consider the quantization of the linearised Einstein-Klein-Gordon sys¬ 
tem on arbitrary backgrounds M^p = (M, g, p), such that M is four-dimensional, 
M = (M, g) is globally hyperbolic, and © = (g, p)* satisfies the coupled Einstein- 
Klein-Gordon equations. Later we point our attention towards such back¬ 
grounds Alp which are in addition of FLRW-type and compare the general 
quantization procedure with the usual approach to the quantization of pertur¬ 
bations in Inflation. 

The linearised Einstein-Klein-Gordon system is a Bosonic gauge theory with 
linearised diffeomorphisms as local gauge transformations. Thus, the strategy 
to quantize this system is to show that it fits in the axiomatic framework of 
general linear gauge theories on curved spacetimes reviewed in Section [2.2.2l In 
particular, we have to show that this model satisfies the axioms of Definiton 

To this avail, we introduce the vector bundles over AI V = T*AI © 
(M X M), where \J denotes the symmetric tensor product, and W = TAI. The 
space of smooth sections of a vector bundle such as V will be denoted by r(V) 
as before. We recall that important subspaces of r(V) are ro(V) and rsc(V) 
the space of smooth sections of compact and space-like compact support, re¬ 
spectively. Both the background fields 0 = [g,pY and their perturbations 
F = are elements of r(V), whereas gauge-transformations (linearised 

diffeomorphisms) will be parametrised by c S r(yV). We introduce on such 
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sections symmetric and non-degenerate bilinear forms by 

(ri,r 2 )v= J dgX (5“^5‘'‘^7l.ac72,M + , 

M 


(3.33) 


(^1 N 2 )w = J dgX • 


M 


(3.34) 


These bilinear forms are well-defined for pairs of sections with compact over¬ 
lapping support. 

The starting point of the analysis is the Einstein-Hilbert-Klein-Gordon ac¬ 
tion for © = {g,(pY G r(V) 



2 



where V ((/?) is an arbitrary smooth potential. In this section, we shall use units 
in which SttG = 1, G being Newtons gravitational constant; this renders the 
Klein-Gordon field, and thus 0 = {g,(pY, dimensionless. The Euler-Lagrange 
equations of 5'(©) are the Einstein-Klein-Gordon-equations 

r(v) 3 EL(e) = (^ I? ) i ^ 

In order to obtain the linearised theory, we split 0 into a background 0 
(which, slightly abusing notation, we denote by the same symbol) and a pertur¬ 
bation r and formally expand (omitting boundary terms) 

5(0 + T) = 5(0) - (EL(0), r)^ - i (r, PT)^ + G(r3). 

Here P is the second order partial differential operator 


p : r(v) ^ r(v) p = 


(3.36) 


(fb7)a/3 = ^ (-V^V^7a/3 + 2V^V(„7,g)^ - 5a/3V'"V'"7^i, - VaV^7/ 
+5a/3V^V^7/ -k gagR^'^lg.^ - Rlag) + Q(V^</?)(V''(^) -k 7 ^/ 3 -^ 
-i<7a/37M.(VV)VV - 

- (V(a(/?) V^)|(() 

P 37 = ((VV)V'^ - + (V'^VV) + 7m. 


{P2Y)o.g = |5a/3 (^(VV)V^ + \dg,V 


Tdlv)cl>. 
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P is formally self-adjoint with respect to (•,-)v even if the background is off- 
shell, i.e. even if EL(0) ^ 0, cf. the discussion in Section 2.2.1]. 

S'(© + r) is invariant under diffeomorphisms of M, in particular under those 
generated by an arbitrary but fixed compactly supported vector field c G ro(W). 
Given such a diffeomorphism, © + F transforms as 


© + r © + r + £,© + £,r + , 


where denotes the Lie derivative w.r.t. <^. To first order in and F, the 
diffeomorphism-invariance of S'(© -I- F) reads 


P£^© = £^EL(©) 


where the term on the right hand side arises from the £,;F contribution of the 
transformed ©-f F. These observations imply the following: we may consistently 
truncate the diffeomorphism-invariant field theory for ©-I-F at joint linear order 
in F and c if and only if we assume that the background © is on-shell, i.e. 
EL(©) = C0; for this reason we shall only consider on-shell backgrounds in the 
following. In this case, one may think of the linearised Einstein-Klein-Gordon 
theory as originating from the quadratic action 


5(2)(r) = _i(F,PF)^ , 


which is for all G Fo(W) invariant under the affine transformation 

F F-f £^©. 

Defining K : F(>V) —>■ F(V) by iFc = £^©, we may express this gauge-invariance 
as 


PoK = 0. 


We recall that this automatically implies that the equation of motion PF = 0 
for the perturbations does not have a well-posed Cauchy problem as non-trivial 
solutions with compact support exist. 

In order to cast the linearised Einstein-Klein-Gordon theory into a form 
which satisfies Dehnition 12.2.51 such that the results of Section 12.2.21 may be 
readily applied, we introduce a field redefinition which generalises the trace- 
reversal well-known from linearised gravity (the appearing numerical factors 
are introduced in order to homogenise the normalisation of the 7-7 and 
components of the principal symbol of P). 



^Strictly speaking Zi.;EL(©) = 0 is satisfied even if Eq = —□(/? + d^V = c with c constant 
but non-zero. However, one may absorb c by redefining 
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Using this field redefinition and observing that it is invertible, we may now 
define 

P = K = ^oK = 0=(^^“/5^=r. 

These definitions are tailored in such a way that the second order action for T 
may now be re-written as 

5(2) (D = -l(r,pr)v = -i(0,P0)v = ^^"^(0) • 

Moreover, S'(^^(0) is invariant under the affine transformation 0 i—>■ 0 + Kc; 
for all c G ro(>V), (•,-)y is symmetric and non-degenerate and P is formally 
self-adjoint with respect to the redefined bilinear form. 

We now provide the expressions for P, K and (•, •)y. As we always assume 
that the background metric and scalar field satisfy the full Einstein-Klein- 
Gordon equations, we have used these to simplify P. 

P : r(V) ^ r(V) P = (3.37) 

(PoOU = + 2V^V(„0^)^ - - 2ve^p 

{P2<l>)^p = <(> 

Pad = 4 + (V“^) O^p 

Pi<^= (-V„V“ + 92y)<^. 

K : r(>V) ^ r(V) Fc = (^2 (3.38) 

V ) 

(0i)02)v = y dgX (^AOi^pO^^ — 29i^^92i,^ + 4’i(t>2^ (3.39) 

M 

In order to satisfy the remaining conditions of Definition 12.2.51 we need to 
check whether R = Ail o AT is Cauchy-hyperbolic and whether there exists a 
‘gauge-fixing operator’ T : r(>V) —>■ r(V) such that the ‘gauge-fixing equation 
of motion operator’ P = P + T o is Green-hyperbolic and Q = o T is 
Cauchy-hyperbolic as well. It is not difficult to check that these conditions are 
satisfied by choosing T = 2 AT, cf. [24]. In particular, with this choice of P, P is 
of the form 

P : r(v) ^ r(V) P = 



(3.40) 
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(Po0)«/3 = + 2(V'^(^)(V(„<^)0;3 )m 

iP24’)a/3 = {VaVpfp)4> 

p^e = 4 (^(v“v'5(^) - e^p 

Pi(l> = (-V„V“ + d^V + (t>. 

Consequently, we have the following result. 

Proposition 3.3.1. The linearised Einstein-Klein-Gordon system defined hy 
{M^,V,W,P,K), where 

• Mip = (M,g,ip) with {M,g) a four-dimensional globally hyperbolic space- 
time and {g, (p) a solution of the Einstein-Klein-Gordon equations (I3.35P 

• V = T*M (B {M X R) and W = TM real vector bundles over M 

• the spaces of smooth sections r(V) and r(>V) of V and W over M are 
endowed with the bilinear forms (•,-)y (13.391) and (•,-)w (|3.34p 

• P is the differential operator P : r(V) —>■ r(V) defined in (13.371) and K is 
the differential operator K : r(W) —>■ r(V) defined in (13.381) 

satisfies the axioms of linear gauge theories in Definition lil.il. 51 

1. (•,-)y and (•,-)w symmetric and non-degenerate. 

2. P is formally selfadjoint with respect to (•, •)y and satisfies P o K = 0. 

3. The differential operator R : r(>V) —>■ r(>V), R = o K, with : 

r(V) —>■ r(>V), •)w = (•, K-)^ is a multiple of a normally hyperbolic 

operator and thus has a well-posed Gauchy problem. 

4 . There exists a differential operator T : r(>V) —>■ r(V), e.g. T = 2K, such 
that P : r(V) ^ r(V), P = P + ToK^ (l3^ and Q : r(>V) ^ r(>V), 
Q = K^oT are normally hyperbolic and have a well-posed Cauchy problem. 

We may now follow the construction reviewed in Section 12.2.21 in order to 
construct a pre-symplectic space of gauge-invariant observables which may then 
be canonically quantized as outlined in Section 12.31 To this end, we consider 
the spaces 


Sol = {0 e r(v) IP0 = 0}, Soisc = Sol n rsc(v), 
g = K [r(w)] , ^sc = e n r,e(v), Ssc.o = k [r,e(w)] . 

The space of gauge-equivalence classes of solutions Sol/C/ may be interpreted as 
the space of pure states in the classical linearised Einstein-Klein-Gordon field 
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theory. The space of (regular) linear functionals on Sol/t/ is parametrised by 
the labelling space 

S = Kero (F^) /P [ro(V)] , Kero (Ft) = {h € ro(V) | K^h = 0}. 

Indeed, the dual pairing 

Sd/e X F 9 ([0], [/]) ^ ([0], [/]) = (0, f)^ 

is by construction independent of the representatives and thus well-defined. £ 
thus conveniently parametrises linear and gauge-invariant local on-shell observ¬ 
ables in the classical field theory. In order to endow £ with a pre-symplectic 
structure, we consider the causal propagator of the normally hyperbolic 
gauge-fixed equation of motion operator P = P + T o and define a bilinear 
form r on F by 

TmAf2]) = {h,E^h)_. 

Proposition 12.2.41 and Proposition 12.2.21 then imply: 

1 . T is independent of the representatives and thus well-defined. 

2 . T is antisymmetric. 

3. T is independent of the gauge-fixing operator T appearing in the definition 
of P. 

4. (F, r) is a well-defined pre-symplectic space. 

An ‘equal-time version’ of the covariant pre-symplectic space (F, r) can be 
constructed by setting 

(01,02)^ = , ®2)y , 

where 0)*' is the future part of 0i. Theorem l2.2.5l then implies that (Solsc/I/sc, 0 ) (•> ■)^) 
is a well-defined pre-symplectic space and that the causal propagator of P de¬ 
scends to an isomorphism 

E^ : (F,r) -)> (Sobc/Ssc.o, (•, •)^) • 

Moreover, given the current j 

7 : r(v) X r(v) ^ t*m (3.41) 

(01,02) ^ 7a(01, 02) = - 40 i^"V„ 02^. + + 20i/V„02/- 

- </>iV„(A2 + Hi - ‘1 o 2’ 

whose covariant divergence satisfies 

V“J„(01,02) = ((01, F02))^ - ((F01, 02 ))^ , 


(3.42) 
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where ((•,-))y denotes the integrand of (t)^, Theorem 12.2.51 further implies 
that 

(ei,02)s^ = J dE 7V“7„(ei,02) 

E 

for an arbitrary but fixed Cauchy surface E with future pointing unit normal 
vector N. 

The held redehnition ” was helpful for uncovering important structural prop¬ 
erties of the linearised Einstein-Klein-Gordon system {M^, V, W, P, K), and to 
show that it hts into the axioms listed in Dehnition 12.2.51 The original sys¬ 
tem {M^,V,P,K) does not directly satisfy these axioms because e.g. the 
operator o K does not have a well-posed Cauchy problem. Yet, the held 
redehnition is merely a computational trick and not of physical signihcance. 
To see this explicitly, we consider the spaces related to the original system 
(M^,V,W,P,X): 

Sol = {T e r(v) I pr = 0}, Soisc = Soi n r,c(v), 
g = K [r(>v)], = e n r,,(v), ^sc.o = k [rsc(w)]. 

£ = Kero {K^) /P [ro(V)] , Kero {K^) = {h € ro(V) | K^h = 0}. 

It is not difficult to see that £ — £ and as before, we may now observe that £ 
parametrises linear functionals on Sol/5. Due to 

(r,h)v = (r, h)^ 

the physical interpretation of these linear functionals in terms of local observ¬ 
ables is manifestly independent of the held redehnition. Finally, it follows from 
the previous discussion that 

P o~ = P + T o o ~ 

is not normally hyperbolic, but still Green-hyperbolic. The causal propagator 
ofPo“is“~^o E^ and the induced bilinear form on £ 

T{[hi], [h 2 ]) = {hi, o Ep] /i 2 )^ 

manifestly equals r. Similarly, (Solsc/5sc,o, ("i •)soi) with 

(0i,02)soi = (-P0m02)y , 

is a pre-symplectic space which is isomorphic to (Solsc/5sc,0) (•> •)s^) via the 
held redehnition map. Gonsequently, we have the isomorphisms 

(Solsc/5sc,0j ("j ■)sjq) — — (Solsc/5sc,0j ("i ■)soi) ■ 

These considerations in classical held theory clearly carry over to the quantum 
theory constructed as in Section 12.31 and the algebras of quantum observables 
induced by the pre-symplectic spaces {£, r) and {£, t) coincide. 
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For the analysis of perturbations in Inflation it is useful to observe that, for 
an arbitrary but fixed Cauchy surface S with future pointing unit normal vector 

N, 

(ri,r2)soi = J dT iv“j„(ri,r2) 

s 

with the current j 

j : r(V) X r(V) T*M (3.43) 

ja(ri,r2) = - - ^7l/V/i72cf 

- - (plVa(l)2 + ^7l/ (Va7>) </>2 + 72a/x 

- ‘1 O 2’ 


whose covariant divergence satisfies 

r2) = ((Fi,pr2))v - ((pFi,F2))v , (3.44) 

where ((•,•))-,; denotes the integrand of This follows from the identity 

j(ri,F2) = j(Fi,r2). 


3.3.2 The Standard Approach to Perturbations in Infla¬ 
tion and Comparison of Approaches 

After discussing the gauge-invariant quantization of the linearised Einstein- 
Klein-Gordon system on general on-shell backgrounds = {M,g,(f), we 
consider the special case of FLRW-type on-shell backgrounds with flat spa¬ 
tial sections, which is the field-theoretic model of perturbations in Inflation. 
We recall that these backgrounds are characterised by 


g = a{r])‘^{-dr]‘^ + dx'^), 


(p{g,x) = (fig). 


Two reoccurring important quantities are 


n = -=aH, 
a 


z = 


a(p' 

IT 


and we shall indicate spatial indices in the homogeneous and isotropic (i.e. 
comoving) FLRW-coordinates by Latin letters i,j,k, _ We shall use the con¬ 

vention that these indices will be raised and lowered by means of the Euclidean 
metric Sij rather than by means of the induced metric a(g)‘^5ij. For the dis¬ 
cussion of perturbations in Inflation it is convenient to work with the original 
formulation ,P,K) of the linearised Einstein-Klein-Gordon system 

rather than its redefined version {M^, V, W, P, K). 

In the following discussion, a special role is played by sections which vanish 
at spatial infinity (with all derivatives). 


roo(V) = {F G r(V) I i9ii • • ■ 9i„r(?7, x) vanishes for |a;| —>■ oo for all n G Nq} 
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roo(W) = {<? € r(>V) I x) vanishes for \x\ —>■ oo for all n G No} 

Namely, one can uniquely decompose F = ( 70 / 3 , (j))*' € FooCF) as 

^ f )2 f 2A {d,B - Vf \ 

7a/3 a(7^) y Q^B-V -2 {d,djC + ) 

A,B,C,D V,W d%=d^Wi = Q 

TeC^ (M, , T/ = 0 , d^Tij = 0 . 

The components B, C, D, Wi are solutions of certain Poisson equations, e.g. 

AS = , A = 

and the uniqueness of the above decomposition results from the unique solvabil¬ 
ity of such equations under the assumption that the solutions vanish at infinity. 
For the same reason, such a decomposition for general F G F(V) can only be 
unique up to harmonic functions. In fact, this non-uniqueness is a non-trivial 
obstacle for proving that a decomposition which is smooth in rj exists in general, 
though we presume that this is the case. Notwithstanding, we shall only need 
the existence and uniqueness of the decomposition for Foo(V) in the following. 

Owing to their transformation properties under the Euclidean group the 
components A, B,C, D,(j) of a section in F G Foo(V) are called ‘scalar’, Vi,Wj 
are called ‘vector’ and Tij are called ‘tensor’ and similarly for the components 
of c G Foo(FV). Following this nomenclature, we say that F = {jap,4’Y G r(V) 

... is of scalar type if jap can be decomposed as above with Vi = Wi= T^ = 

0 . 

... is of vector type if (/> = 0 and 7 ^^ can be decomposed as above with 
A = B = C = D = Tij = 0. 

... is of tensor type ii (j) = 0 and 70^3 can be decomposed as above with 
A = B = C = D = Vi = Wi=0. 

Based on this we define the following section spaces. 

= {F G F(V) I F is of scalar/vector/tensor type} 

= roo/o(v) n F^/'"/^(v) 

The existence and uniqueness of the decomposition of sections vanishing at 
spatial infinity may be cast in the following form 

roo(V)=F^(V)©F^(V)©F^(V) 

and one may check that the decomposition is orthogonal w.r.t. (•, •)v and (•, •)w. 
However, the decomposition is non-local in space, and thus one has 

F^(V)©Fo^(V)©f;^(V)CFo(V)! 
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Similar decomposition results hold for G roo(W) and we extend the above 
nomenclature to this case in the obvious way. 

By the existence and uniqueness of the decomposition, the individual compo¬ 
nents induce well-defined functionals on roo(V) and roo(W), e.g. A : roo(V) —>■ 
C^ (M, R). This existence and uniqueness further implies that there exist pro¬ 
jectors 




which are formally selfadjoint w.r.t. (•, •)v and (•, •)w- It is not difficult to check 
that the gauge transformation operator K and the equation of motion operator 
P commute with these projectors, i.e. 



Thus the equations of motions and gauge transformations decouple for sections 
which vanish at spatial infinity and we may consider subspaces 


Sol^/^/^ C Soloo and C G^ 


which are defined in the obvious way. 

The decomposed equations of motion may be expressed in terms of gauge- 
invariant linear combinations of the decomposition components, i.e. in terms of 
functionals on Soloo/f/oo- 


^^J = A-idr,+n){B + C') <I> = D-n{B + C') 


X = (t^-+ E') f^ = --q> + x 


a 


X^ = W'- K Tj 


and <i> are the so-called Bardeen potentials whereas /i is the Mukhanov-Sasaki 
variable and is of particular physical significance because it is related to the 
perturbation of the scalar curvature of the spatial slices. In terms of these 
gauge-invariant quantities, the equations of motion PT = 0 for T G roo(V) read 
(see [23] for details): 


scalar: 




df] azp -I- Ao , 



where 770 is arbitrary and Aq is the unique solution of 
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vector: 


tensor: 


AXi = 0, {dn + 2'H)X, = 0 
P^T,, = ^ ((d^ + 2n)d^ - A) T,, = 0 


We see that /i is a conformally coupled scalar field with a particular time- 
dependent mass and that $ and x ^'^e non-local functionals of /i, whereas 
no non-trivial vector solutions vanishing at spatial infinity exist. Moreover, 
the tensor components Ty satisfy a normally hyperbolic equation. The last 
statement may have been deduced directly from the previous analysis by ob¬ 
serving that the field redefinition “ acts trivially on r^(V) and that r^(V) lies 
in the kernel of . Thus P restricted to r^(V) coincides with the normally 
hyperbolic P = P o~~^ -\-T o . 

The standard treatment of the quantum theory of perturbations in Inflation 
may be rephrased in the present context as follows. One considers the symplectic 
spaces 




X 9 ([/l], [h]) H9 r'^([/i], [/2]) = 

E^ causal propagator of P^ 


and 



M 


P^ and are both formally self-adjoint and one may show that is 

non-degenerate on C^{M,T) [24]. Thus and are antisymmetric and 
non-degenerate by standard results reviewed in Section 12.2.11 

In the standard treatment of perturbations in Inflation one effectively as¬ 
sumes that all gauge-invariant (polynomial) local observables in the quantum 
theory are spanned by the local observables obtained from the canonical quan¬ 
tization of {E^,t^) and On the other hand, one may take the point 

of view that the construction outlined in the previous discussion of the lin¬ 
earised Einstein-Klein-Gordon system, i.e. the canonical quantization of the 
pre-symplectic space (£, r) should yield all (polynomial) local gauge-invariant 
observables. In order to compare these two approaches, we define 


Kero/^/"^ (a:^) = Kero (X^’) O r'^/^/^(V) 

gS/V/T jp 


/ 


One may then prove the following results [24] . 
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Theorem 3.3.1. The following relations hold. 

1 gS/V/T ^ g 

2 . = { 0 }. 

3. and are isomorphic. 

4- and (^£^^are isomorphic. 

5. £^ ®£'^ C£. 

6. £^ 0 is separating on Soloo/Soo = 0 Sol^/5^. 

7. T is non-degenerate on £. 

The proofs of most of these statements are unfortunately quite cumbersome. 
In particular, one needs to analyse how the solution theory of hyperbolic op¬ 
erators and the results of Theorem 12.2.51 intertwine with the decomposition of 
perturbations into scalar, vector and tensor contributions. Once this is under¬ 
stood, the third and fourth statement e.g. follow essentially from the fact that 
the equal-time versions of the compared symplectic spaces are isomorphic be¬ 
cause the corresponding symplectic forms are given by Cauchy surface-integrals 
of (normal components of) ‘conserved currents’ which can be seen to match by 
explicit computations. 

Even if the proof of the above theorem is rather lengthy, the statements have 
a straightforward physical interpretation, which is valid both for the classical 
and the quantum theory of perturbations in Inflation as should be clear from the 
discussion of the quantization of pre-symplectic spaces in Section 12.31 The first 
statement implies that there are local observables which may be meaningfully 
classified as ‘scalar’ and ‘tensor’. This may seem surprising in view of the fact 
that the decomposition of configurations in scalar/vector/tensor components is 
a priori non-local. However the second statement entails that there are indeed 
no non-trivial local ‘vector’ observables. From the third and fourth statement 
one can infer that the standard treatment of perturbations in Inflation captures 
the same local scalar and tensor observables that one obtains from the general 
gauge-invariant quantization of the linearised Einstein-Klein-Gordon system. 
An interesting result found in [18] implies that the scalar field /i is in fact the 
unique field with a normally hyperbolic equation of motion whose associated 
symplectic space is equivalent to (£’'^,r). 

However, statement 5. implies that not all local observables of the linearised 
Einstein-Klein-Gordon system are spanned by local observables of scalar and 
tensor type. In this sense, the standard approach to the quantization of pertur¬ 
bations in inflation ‘misses’ some local observables. However, the sixth state¬ 
ment entails that the observables captured in the standard approach are still 
sufficient to measure configurations of the perturbations which vanish at spatial 
infinity. These configurations are considered to be ‘small’ in a certain sense. 
Presumably this statement can be generalised by proving that local observables 
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of scalar and tensor type separate quantum states whose correlation functions 
vanish at spatial infinity in each argument. 

Finally, the last statement is somewhat independent of the others and may 
be interpreted such as to say that the quantum theory of the linearised Einstein- 
Klein-Gordon system on FLRW backgrounds does not contain non-trivial clas¬ 
sical observables. 
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